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0.1 Foreword 



Classical Electrodynamics is probably the most fascinating and complete 
part of the Classical Field Theory. Intuition, free thought, perspicuity and 
research skill of many years finally brought about the synthesis of experi- 
ment and theory, of physics and mathematics, which we have been calling 
for short the Maxwell equations for the duration of a century and a half. 
From the beginning of the second half of the 19th century till its end these 
equations turned from abstract theory into daily practice, as they are today. 
Their profound study during the first half of the 20th century brought for- 
ward a new theoretical concept in physics known as relativism. Brave and 
unprejudiced workers enriched and widened the synthesis achieved through 
the Maxwell equations, and created a new synthesis known briefly as quan- 
tum electrodynamics. Every significant scientific breakthrough is based on 
two things: respect for the workers and their work and respect for the truth. 
" May everyone be respected as a personality, and nobody as an idol" one of 
the old workers used to say. We may paraphrase that saying: "may every 
scientific truth be respected, but no one be turned into dogma" . 

In this book I tried to follow the values this creed teaches, as far as my 
humble abilities allow me to. Together with the analysis of the classical elec- 
trodynamics and the quantum concept of the structure of the electromagnetic 
field, the path followed brought me to the conclusion that a soliton-like solu- 
tion of appropriate non-linear equations characterized by an intrinsic period- 
ical process is the most adequate mathematical model of the basic structural 
unit of the field - the photon. The fact that neither the Maxwell equations 
nor the quantum electrodynamics offer the appropriate tools to find such 
solutions, unmistakably emphasizes the necessity to look for new equations. 
The leading physical ideas in this search were the dual ("electro-magnetic") 
nature of the field on the one hand and the local Energy-Momentum Con- 
servation Laws on the other. The realization that every such soliton-like 
solution determines in an invariant way its own scale factor, as well as the 
suitable interpretation of the famous formula of Planck for the relation be- 
tween the full energy E of the photon and the frequency u = 1/T of the bc- 
forementioned periodical process, which I prefer to write down as /i = E.T, 
advanced Planck's constant h to the rank of an estimator, separating the re- 
alistic soliton-like models of the photon from among the rest. The resulting 
soliton-like solutions possess all integral qualities of the photon, as described 
by quantum electrodynamics, but also a structure, organically tied to an in- 
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trinsic periodical process, which in its turn generates an intrinsic mechanical 
momentum - spin. I consider this soliton- like oscillating non-linear wave 
much more clear and understandable than the ambiguous " particle-wave " 
duality. 

The dual 2-component nature of the field predetermined to a great extent 
the generalization of the equations in the case of an interaction with another 
continuous physical object, briefly called medium. The proposed physical 
interpretation of the classical Frobenius equations for complete integrability 
of a system of non-linear partial differential equations as a criterium for the 
absence of dissipation, turned out to be relevant and was effectively used. The 
fruitfulness of the new non-linear equations is clearly shown by the family of 
solutions, giving (3+l)-dimensional interpretation of all (l+l)-dimensional 
1-soliton and multisoliton solutions. I have to say that the use of difi^erential 
geometry proved extremely useful. 

This book is addressed to all who love theoretical physics and try to 
build up their own point of view while in the same time show due respect for 
others' opinions. The stress is laid on the conceptual and generic framework, 
while in many cases the intermediate calculations were omitted. I did not 
propose examples of complete description of actual systems, as this was not 
my purpose. A feature of this book is the lack of citations. In my opinion 
the reference list in the end will be enough for that purpose. 

1 would like to express my most sincere gratitude to all my friends and 
colleagues, with whom I discussed to one extent or another the issues and 
results herein presented. Each and every conversation was extremely valuable 
to me both as an actual analysis of the issue and as a stimulus for its further 
study. I would highly appreciate any remarks and opinions concerning this 
book from anyone interested in the subject matter being analyzed. 
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0.2 To the Reader 



Dear Reader, 

Opening the first pages of a new book, devoted to the part of science 
you hke and want to know as better as you can, the first and most natural 
question you ask is: what new shall I learn if I read this book? Having in 
view the huge quantity of monographs and other issues, devoted to every 
of the various and numerous branches of physics, the natural respect to you 
requires a short but sufficiently informative and fair response to this question. 
And this is the contents of these preliminary notes. In other words, I'll try 
to explain in short what is this extended electrodynamics and what aims to 
achieve this extension of the well known from the University electrodynamics 
of Faraday and Maxwell. 

First of all, let me specify the kind of extension I mean. In theoretical 
physics by means of mathematical relations and equations some of the really 
existing objects and processes are modelled and described. An equation 
separates those values of a given quantity, which are of some interest to that 
person, who has written this equation. For example, the elementary equation 
X — 1 — Q separates the value 1 of the variable x. If, for some reasons, we 
want to separate two values of the variable x, e.g. (1,-1), we, naturally, 
write down a new equation, which has two solutions: 1 and (—1). The new 
equation will be — 1 = {x + l)(x 1) = 0. In this way we extended 
the set of admissible values of the variable x, or we may say that we have 
extended the equation x — 1 = to the equation x'^ — 1 — 0. The idea 
is now quite clear: we write down new (algebraic or differential) equation, 
whose solutions comprise as a subclass all solutions of the old equation, but 
have also new ones. The important point is the new solutions to have the 
properties, desired by us. 

The basic equations of Classical electrodynamics are the well known equa- 
tions of Faraday-Maxwell. So when we talk about extended electrodynamics 
it is clear that just these equations we have to extend in the sense, mentioned 
above. In order to motivate such an extension we should honestly say two 
things: which properties of Maxwell equations we do not like sufficiently, and 
what new properties will be required from the new solutions of the extended 
equations. Let me briefly discuss this matter. 

As it is well known, at zero electric current Maxwell's equations require 
the components of the electric and magnetic vectors to satisfy the wave equa- 
tion (d'Alembert's equation). The solutions of this equation in the whole 3- 
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space have the following 2 properties: first, every solution can be represented 
by a sum (finite or infinite) of a subclass of solutions, known as plane waves; 
second, every solution, defined by some finite localized and smooth enough 
initial condition, "blows up" radially and goes to infinity with the velocity of 
hght. The plane waves are characterized by the condition, that there exists 
a system of canonical coordinates on TZ^ in which these solutions depend 
on one space variable only. This means, that at every moment every such 
solution occupies the whole 3-space, or an infinite its subregion. Therefore, 
they are infinite and their integral energy is infinitely large. We note that 
this infinity does not come from some singularity of the solutions. We set 
the qTicstion: do there exist in Nature electromagnetic fields with properties 
adequate to these exact solutions called 'plane waves' ? The natural answer 
is negative, since in order to create such an object it will be necessary to 
transform an infinite quantity of some kind of energy into electromagnetic 
energy, which will take infinite period of time. 

On the other hand, according to the second property, every finite solution 
is strongly time-unstable, so it could hardly be considered as a realistic model 
of a real object. As for the static case, the components of the electric and 
magnetic vectors satisfy the Laplace equation and, as it is well known the 
solutions of this equation in the whole 3-space are singular, or when they 
are finite and nonsingular, they are constant. We conclude, that the vacuum 
Maxwell equations can not describe finite and time-stable, i.e. soliton-like 
electromagnetic formations. And our purpose is to describe just such soliton- 
like electromagnetic formations. 

If we now recall the basic idealization of classical mechanics, the material 
point, we'll see that it is a full antipode of the plane wave: the material point 
has no structure, occupies zero volume and has finite energy. Are there such 
objects in Nature? According to my opinion the real objects are finite, i.e. 
at every moment of its existence every really existing object occupies finite, 
comparatively small volume and it has definite properties of constancy and 
stability. When they are subject to external perturbations they survive or 
transform into other objects, obeying definite conservation laws. Usually, 
a survival is connected with some change in the behavior of the object as a 
whole, but it also keeps some of its basic characteristics unchanged, otherwise, 
we could not say that the object is the same, i.e its identification after the 
perturbation should be possible. 

At the beginning of our century it has become clear that the electro- 
magnetic field has a discrete nature, i.e. it consists of many non-interacting 
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(or very weakly interacting) objects, carrying energy, momentum and in- 
trinsic angular momentum. Moreover, the integral energy of these objects 
originates intrinsically from a periodic process with frequency of v in corre- 
spondence with the Planck's formula E — hv. The experiment shows soon 
that these objects are finite, they move as a whole with the velocity of light 
along straight lines, carry momentum hvjc and intrinsic angular momentum 
h = ET = E/v. The problem to describe such (free) objects, called later 
photons, appeared. Because of the finite nature and time-stabihty of photons 
it is clear that the solutions of the wave equation can not give good enough 
mathematical models of them. On the other hand, the material point (or 
particle) can not also serve as a model, since there is no any sensible way to 
assign the characteristic frequency to a free particle. The frequency always 
comes from an outside elastic force, so it depends strongly on this force and 
is not a proper characteristic of the particle. The classical free particle moves 
always along straight lines with a constant velocity. There have been some 
unsuccessful efforts to build an extended model of the photon, but pressed by 
the quickly advancing experiment, the major part of the physicists assume 
the prescription schemes of quantum electrodynamics. A basic assumption 
in these schemes is that the photon is a point-like object, so its frequency 
(or spin) is an integral characteristic of the same kind as the proper mass 
and the electric charge. I do not share this point of view, since a free struc- 
tureless object can not carry the physical characteristic of intrinsic angular 
momentum and no periodic process can be associated with its existence. I 
think that a periodic process and an intrinsic angular momentum can be as- 
sociated with a free object only if it has structure, and these characteristics 
may have finite values only if the object is finite. The photon moves uni- 
formly as a whole, but it is not a point-like (or structureless) object, and its 
existence is strongly connected with some intrinsic periodic process with a 
constant frequency. This process occurs in the whole volume and in this way 
generates the intrinsic angular momentum. Since at this stage we do not see 
any other more-natural way to combine the known photon's features than 
the notion of soliton-like objects, we turn to appropriate (3+l)-soliton-like 
solutions of a definite system of nonlinear partial differential equations as 
possible mathematical models. 

As it was mentioned earlier. Maxwell's equations have no such solutions. 
On the other hand, as an working tool in macrosituations they have proved 
their wide applicability, so it does not seem reasonable to leave them off 
fully and to search for new equations. Such a step would reject also the 
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corresponding well proved energetic relations for finite volumes, obtained by 
these equations. Therefore an appropriate their extension seems natural and 
reasonable provided that the new incorporated solutions have the desired 
properties. 

In the non-vacuum case, i.e. in the presence of energy-momentum ex- 
change between the field and some medium, there are also problems coming 
mainly from the hypothesis that the Faraday induction law is valid for all 
media, which leads to the notion that the field is able to exchange energy- 
momentum only with electrically charged particles. No doubt, such media 
exist, but the assumption that this is true for all media, we do not accept. 
Further, the usual way to treat the field through introducing the polarization 
and magnetization vectors and to describe the media by means of corre- 
sponding permeability tensors and bound currents and so on, may be good 
operational skills, but they do not lay on approved physical principles and do 
not generate new ideas. In this case, as well as in the vacuum case, we choose 
a different way, namely we extend the equations in such a way, that to keep 
everything, that classical electrodynamics is able to do, and to incorporate 
new solutions with good enough properties. Let's say quite clear, that our 
purpose in this case is the same as in the vacuum case, namely, to describe 
soliton-like field configurations with sufficiently clear physical interpretation. 

In the both cases, the leading idea for extension of the equations is an 
analysis and appropriate formulation of the local energy-momentum conser- 
vation laws in relativistic terms. In result we obtain nonlinear equations. 
In the continuous medium case an additional and entirely new moment is 
considered, namely, a physical interpretation of the Frobenius integrabil- 
ity equations for Pfaff systems. Also, the classical concept of a continuous 
medium is extended and understood as an arbitrary continuous physical sys- 
tem, exchanging energy-momentum with the field. The encouraging results 
we obtain in the both cases give us reasons for hope that we are on the right 
way. 

Finally, everyone, who decides to read this book, will follow the complete 
version of the author's mental way and will get to the results obtained fol- 
lowing an easy and smooth way of reasoning. The real road, that 1 passed 
through, was very different from what is given here. There were many turns, 
unexpected traps and various positive and negative surprises, and all this 
was taking place, when the fashion in theoretical physics was quite differ- 
ent. The deep belief in the conservation laws, in their universal character, 
was the point of light, that was leading me through the jungle of unknown 
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possibilities and hard to evaluate hypotheses and helped me to withstand 
the falling on one after another fashionable topics. Writing this book I was 
doing my best to pay a maximum respect to every reader, to every positive 
opponent. And, in order not to appear any doubts about my great respect 
to the prominent men, who created Classical electrodynamics, I would like 
to give here the creed that was the basic stimulus during this period: the 
respect and esteem paid to the creators can not be honest and genuine if they 
are not in correspondence with the respect and esteem paid to the Truth. 

Now, De£ir reader, after you have got some idea of what could be learned 
from this book, and if you have already made up your mind to become well 
acquainted with the small harmonic world I tried to create, turn over this 
page and be my fair corrective up to the last line of the last page. 
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Chapter 1 



From Classical Electrodynamics 
to Extended Electrodynamics 



1.1 Basic Notions and Equations of 
Classical Electrodynamics 

1.1.1 Nonrelativistic approach 

In the mathematical model of the electromagnetic phenomena, called shortly 
CLASSICAL ELECTRODYNAMICS (CED), the Newton's notions of space 

and time are assumed, namely, a mathematical model of the real space is the 
real 3-dimensional mathematical space TZ^, considered as a 3-dimensional 
Euclidean manifold with the standard Euclidean metric g, and a mathemat- 
ical model of the real time is the 1-dimensional mathematical space TZ. The 
standard orthogonal coordinates {x, y, z) are usually used, so that the metric 
tensor Qik has canonical components, i.e. the diagonal elements are equal to 
1, and all other components are zero. Then the isomorphism between the 
tangent vectors and the 1-forms, defined by g and denoted by the same letter 
g, (or g~^), looks as follows: 



where Vk = gnV^ . This isomorphism is extended naturally (i.e. component- 
wise) for arbitrary tensors (or tensor fields). The metric g defines also a 
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volume element: 



ljq = y \detgik\dx A dy A dz = dx A dy A dz, 

a covariant derivative V, which in standard coordinates reduces to the usual 
derivative, and the Hodge *-operator * : A^CR,^) —>■ A^~p{TZ^) according to 
the rule 

aAp = g{*a, P)u;o, 

where a and f3 are p and (3 — p) forms respectively. The operator * is linear 
and in the canonical orthonormal basis {dx,dy,dz} from the above defining 
relation we obtain 

*1 = dx A dy A dz, *dx = dy A dz, *dy = —dx A dz, *dz = dx A dy, 

*dx Ady — dz, *dx Adz — —dy, *dy Adz — dx, *dx Ady Adz — 1. 

We note that the defining relation for the * operator is equivalent to the 
following: if a and (5 are p-forms then 

a A*I3 = g{a, P)wo. 

If d is the exteriour derivative, then the operator S = (—1)^ d* is called 
coderivative and we get {5a)ik... = —'^j(^\k ■ easily checked, that the 

operator 5 is the dual to d in the vector space of forms with compact support 
with respect to the inner product {a, P) = J a A and for forms with 
different degree is assumed {a,l3) — 0. 

From physical point of view CED assumes the following: 

1. The Electromagnetic (EM) field is a real object (or a set of real objects) 
possessing energy and momentum. 

2. The EM-field is able to interact (i.e. to exchange energy and momen- 
tum) with material particles, possessing the characteristic electric charge 
which takes positive, negative and zero values. Such charged particles are 
called field sources. 

3. The sources arc called 

Si) free - if under the action of the external field they can move throughout 
the whole space, and 

b) bounded - if their motions are bounded inside comparatively small space 
regions. 

In order to describe the EM-field CED introduces the following quantities: 
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I. A couple of vector fields {E, B), defined on TZ^, and a parametric de- 
pendence of {E, B) on the time t is admitted. 

II. A scalar function t), called charge density, such that the inte- 
gral / pujo, computed over the region V C TZ^, gives the entire electric charge 
in V. 

III. A vector field \{x,y,z,t), describing the mechanical motion of the 
charge-carriers, and the vector field j = pv is called electric current. 

Remark. Further a vector field V and the corresponding 1-form g{V) will 
be denoted by the same latter, since from the context it will be clear if the 
object is tangent or co-tangent. 

An important concept in CED is the so called flow of a vector fi,eld U 
through a given 2-dimensional surface S. By definition, it is the integral of 
the 2-form *U over S: Js*U. It is important if the surface is closed (most 
frequently homeomorphic to the 2-sphere S^), or it is not closed and has for 
a boundary a given 1-dimensional manifold. 

As a generalization of the experimental data in CED is assumed the 
following: 

1°. The flows of the electric E and magnetic B fields through a closed 
2-surface S"^, surrounding the 3-volume V are defined as follows: 

<p *E — A-K I pujo, f *B — 0. 
Js^ Jv Js^ 

2°. // the 2-surface S is not closed and its boundary is the closed curve I, 
then the following relations hold: 

^ [ *E^c [b-Att [ *i, ^Ub^-cIe. 
at Js Ji Js at Js Ji 



Prom these relations, making use of the notations g ^ * dg — rot, 
g~^Sg — div and the Stokes' formula 



Is^^ I I 



a 



we get the differential equations 



IdE Att 

- =rotB j, divB = 0, (1.1) 

c at c 
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1 dB 

-— = -rotE, divE = At: p. (1.2) 

c dt ' ^ ^ 

Because of the identity diVoTot = 0, from the first equation of (1.1) and the 
second equation of (1.2) it follows the continuity equation 

^ = -divi, (1.3) 

the sense of which becomes clear from its integral form 

d 



-I 

dt J\ 



V JSv 

This relation shows, that any change of the electric charge inside the region 
V is caused by some processes of charge transport through the boundary of 

V , i.e. charges do not vanish and do not arise. 

Using the scalar product g from the first equation of (1.1) and the first 
equation of (1.2) wc obtain 

\^{^' ^) + ^K^' ^) = - 9{B,rotE) - %(£;j). 

Since g does not depend ot time and 

g{E, rotB) - g{B, rotE) = -div{E x B) 

we get 

d E^ + B^ 

= —\E — divS, 

dt St: 

where the vector 

S = -^E X B 

is the Poynting vector. This last relation describes the local balance of the 
energy and momentum in the system EM-field and free current in an unit 
space-time volume. 

The equations (1,1), (1,2) are linear, so that if {Ei,Bi) and (£'2,^2) are 
two solutions, then any linear combination 

E^aEi + bE2, B = mBi + nB2 

with constant coefficients is again a solution. The following question arises 
naturally: do there exist constants (a, b, m, n) , such that the linear combina- 
tion 

E' = aE + bB, B' = mE + nB 



4 



is again a solution? The answer to this question is positive iS m = —b,n = a. 
The new solution will have energy density w' and momentum S' as follows: 

S'={a'' + b^)-^ExB. 

An 

Obviously, the new and the old solutions will have the same energy and 
momentum if + 6^ = 1. 

These observations show that besides the usual linearity, Maxwell's equa- 
tions admit also " cross" -hnearity, i.e. linear combinations of E and S of a 
definite kind define new solutions. Therefore, the difference between the elec- 
tric and magnetic fields becomes non-essential. The important point is that 
with every solution {E, B) of Maxwell's equations a 2-dimensional real vector 
space, spanned by the couple {E,B), is associated, and the linear transfor- 
mations, which transform solutions into solutions, are given by matrices of 
the kind 

a b 
—b a 

If these matrices are unimodular, i.e. if -|- 6^ = 1, then the initial and the 
transformed solutions have the same energy and momentum. It is well known 
that matrices of this kind do not change the canonical complex structure J 
in 71^: if the canonical basis of TZ^ is denoted by (ei, 62) then J is defined by 
J{ei) = 62, J(e2) = -ei. 

The above remarks suggest to consider E and B as two vector-components 
of an 7?.^ -valued 1-form Q: 

fl ^ E ^ ei + B ^ 62- 



So, the current j becomes 1-form J 
charge density becomes a function Q 
equations take the form 



= j (8) ei with values in V?, and the 
— p<Si ei with values in TZ"^. Maxwell's 



1 dfl An 

= J-*dJ{n), Sn = AnQ, (1.4) 

c at c 

where J(Q) = E ® J{ei) + B J{e2) — E®e2— B ®ei. Note that according 
to the sense of the concept of current in CED and because of the lack of 
magnetic charges, it is necessary to exist a basis of 7?.^, in which J and Q to 
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have components only along ci. Nevertheless, this point of view shows that 
even at this non-relativistic level the separation of the EM-field to "electric" 
and "magnetic" is not adequate to the real situation. 

1.1.2 Relativistic formulation 

We pass now to the relativistic formulation of CED. We begin with the note 
that the relativism, considered as a theoretical conception for understanding 
and modeling the natural entities and processes, arises as a result of the anal- 
ysis of the invariance properties of Maxwell's equations with respect to linear 
transformations of space-time coordinates {x,y,z,^ = ct). The assumption 
for the linear character of the space-time transformations comes on the one 
side from the linearity of the model space TZ"^ and Maxwell's equations, and 
on the other side, from the idea, that any straight- line uniform motion should 
not affect the character and the course of all physical processes. The latter 
formally means, that the parameters of the admissible space-time transfor- 
mations, interpreted as straight-line uniform motions, should be determined 
by the relative constant velocity between two frames. So the cartesian fram- 
ings in TZ"^ model the physical inertial frames. The basic conclusion of the 
analysis carried out during the end of last and the beginning of the current 
century consists in, that the invariance of Maxwell's equations requires a 
pseudo-eucledean space-time metric tensor rj and, additionally, any compo- 
nent of E and B depends linearly on all components of E and B when a 
space-time transformation is performed. (In connection with this we note, 
that the "cross" -linearity mentioned above admits "mixing" only for corre- 
sponding components of E and B: Ei with E2 with B2 and E?^ with ^3.) 
This undoubtedly shows, that the adequate mathematical object, describing 
the i?M-field, must have 6 independent components and its transformation 
properties should be determined entirely by the admissible space-time trans- 
formations. In view of the 4-dimensionality of space-time such objects are 
only the antisymmetric tensor fields of second order. Because of the possible 
identification of the contravariant and covariant tensor fields by means of the 
pseudo-metric tensor rj^j^y^ the obvious candidates for models of the E'Af -field 
are the differential 2-forms. So, starting with the 3-dimensional vector fields 
E and B and with the pseudoeucledean structure 77^1^ of the model space TZ^ 
we have to build a differential 2- form F G K^ilZ'^, rj). Further, the space-time 
{V}, rj), where rj^^ = — 1 for /i = z/ = 1, 2, 3; 7744 = 1, and rj^^, = for ^ z/ 
in standard coordinates {x^ , , , x^) = {x,y,z,^ = ct), will be denoted by 
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M and will be called Minkowski space. The pseudo-metric r^^j, defines in the 
same way a volume element cuo and the Hodge *-operator. We have 



cuq = \J\detrj^j_y\dx A dy A dz A d^, 

a A l3 = rj{*a, I3)ujq <^=^ a A*I3 = —r]{a, (3)ujq, 

**p = -{-ly^^-^hd, = -i-iy^^-^hp, *u;o = 1, *1 = -uo, 

*dx — dy A dz A d^ *dx A dy A dz = d^ 

*dy = —dx Adz Ad^ *dx A dy A d^ = dz 

*dz = dx A dy Ad^ *dx A dz A d^ = —dy 

*d^ = dx A dy A dz *dy A dz A d^ = dx 

*dx Ady = —dz A d^ *dy A dz = —dx A d^ 
*dx A dz = dy A d^ *dy Ad^ = —dx A dz 
*dx A d^ — dy A dz *dz Ad^ — dx A dy. 

An arbitrary 2-form F on M in standard coordinates looks as follows: 



1 

2' 

+Fudx Ad^ + F24dy Ad^ + Fudz A d^. 



F = -F^^dx^ A dx" = Fi2dx Ady + Fisdx A dz + Fasdy Adz + 



Then for *F we obtain 
1 

-F22,dx Ad^ + Fisdy Ad^- F^dz A d^. 



* F = --s^i.arF'^'^dx^ A dx" = F2,Adx Ady — F2idx Adz^- F^dy A dz 



The definition of the components F^jy by means of the components of E and 
B is made in the following way. Let i:TZ^ — > TZ'^ be the standard immersion 
(x, y, z) — > (x, y, z, 0). Then we define i*F and i*{*F) by 

i*F^*B, i*{*F)^*E, (1.5) 

where on the right-hand side of these equalities the Eucledean *-operator is 
used. Relations (1,5) define F uniquely, and we get 
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-^3; -^13 — —B2, F23 — Bi, Fi4 — El, F24 — E2, F34 — F3 



Recalling that in the static case with zero current Maxwell's equations reduce 
to dE = 0, d * F = 0, dB = 0, d* B = 0, and the well known relation 
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di* = i*d, we obtain for this static case i*dF — 0, i*d * F — 0. The map i* 
cancels all terms where stays. Removing this map i*, we get the equations 
dF = 0, d * F = 0, so we keep all terms with non-canceled, and having 
in view the above component interpretation of Ffj,„ we obtain exactly the 
left-hand sides of Maxwell's equations (1,2), (1,1) respectively. We introduce 
now the 4-current j'* by j'^ = pu'^, where u'^ is the 4- velocity of the charged 
particles. Maxwell's equations take the form 



dF = 0, d * F = 47r * J. 
These equations (1.6) may be written in various forms: 

S^F = 0, 6F = Anj 



(1.6) 



dF. 



+ 



dF,„ dF 



+ 



CFfl 



= 0, 



dx'^ dx^^ dx^ 
Of definite importance for the theory is the quantity 

in 11 

C ^ F ^F^/^S" - F F""" = 



-F^F^'^ - {*F)^,{*Fy 



since on the solutions of (1.6) the following relation holds: 



F,,{5Fy + {*F),4S*Fr 



(1.7) 



(1.8) 



This relation describes the local energy-momentum balance in the system 
EM-field and free current. The quantity F^yj'^ is the Lorentz force and it 
determines the energy-momentum, which the charge carriers get from the 
field in an unit 4-volume. In regions with zero current = we have 
VvQ^^ = 0, so we may create integral conserved quantities. Because of its 
importance for the theory we shall consider this point more in detail. 

In classical field theory we build integral conserved quantities, i.e. time- 
independent quantities, by means of a symmetric second rank tensor 
with zero divergence V^Q'^^ = by making use of isomctries, i.e. symmetries 
of the metric tensor, in the following way. By definition a symmetry of a 
tensor field g is a map f : M ^ M, which keeps this tensor field unchanged: 
f*g = g. When an one-parameter group of symmetries ft, defined by the 
generator X, (or the vector field X) is given, then X is called local symmetry 
of g and the following relation holds 



Lxg = limt^o 



ft9 - 9 
t 



0, 
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which means that the local symmetries of g are those vector fields X along 
the integral lines of which g stays unchanged. The expression on the right is 
called Lie-derivative of g along X. The local symmetries of the metric tensor 
are also called Killing vector fields. The equation Lxg — 0, where g is given, 
looks as follows 

where V is the corresponding symmetric Riemannean connection. If now 
Q^u is a conservative tensor field, i.e. V^Q'j^ — 0, and X is a local isometry, 
we obtain 

V^iQ'^X^) = (V.QpX'^ + Q'^'V^X, = Q'^'^V^X,. 

Because of the symmetry of Q, in the sum Q'^'^Vf^X^ only the symmetric part 
of V^Xi, may contribute, but this symmetric part is zero since X is a local 
isometry. In this way with every local isometry X of the metric the 1-form 
Q^ijX'^dx" is associated, and this 1-form has zero divergence. The last means 
that the 3-form *{Q^i,X'^dx'') is closed, so according to the Stokes theorem, 
the integral of this 3-form over TZ^ will not depend on time. Of course, these 
considerations make sense only for finite valued such 3-integrals, i.e. for finite 
field functions. 

In order to complete the energy-momentum balance picture we have to 
point out how the charged particles of the 4-current use the gained from the 
field energy-momentum to change their behaviour, i.e. we have to write down 
the equations of motion of these charge carriers. 

Assuming that only interaction between the particles and the field takes 
place, the energy-momentum tensor of the particles is defined by 

where Ho denotes the invariant mass density of the particles. So, the full 
local energy-momentum conservation law requires V^iQ]^ + K^) = 0. Since 
particles do not vanish and do not arise, which formally means that the mass 
continuity equation VuiHoU'^) — holds, we obtain 

HoC^u^'VyU^ = -F^^f. (1.9) 

This equation for the vector field u describes the mechanical evolution of the 

charge-carriers. We note that this is a compact form of a nonlinear system 
of partial differential equations for the components of m, while Maxwell's 
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equations are linear. This shows that the current j'^ = fm^ cannot he de- 
fined independently^ i.e. it strongly depends on the field F . Therefore, the 
variational procedure for verifying Maxwell's equations with "given and not 
subject to variation" current does not seem to be quite correct. Otherwise, we 
lose the full energy-momentum conservation law, which is hardly preferable. 

Remark. It is clear that the particles with mass distribution /io and charge 
distribution p, are considered as sources of the field F according to the usual 
interpretation of equations (1.6) amd (1.7), and, therefore, they cannot he 
considered as test particles, i.e. as not disturbing the field. If the field is not 
disturbed, i.e. if it does not exchange energy- momentum with the particles, 
then VjyQ^ = and the corresponding solution of the field equations, which 
is meant to define "Lorentz force", has to satisfy the current- free Maxwell 
equations dF = 0, = 0, as in the static Coulomb case. In this case the 
dynamical equations for the particles appear as an additional assumption, 
and talking about full energy-momentum conservation law is hardly sensible 
from the point of view of the field, nevertheless, it may get some sense from 
the point of view of the particles. 

We consider now the conformal invariance of pure field Maxwell equa- 
tions. In accordance with (1.6) if j = we get 

dF = 0, d*F = 0. (1.10) 

The only factor not permitting a full invariance of these equations is the 
*-operator, more exactly, its restriction on 2-forms. On a 2n-dimensional 
riemannian manifold the restriction of * on n-forms is always conformally 
invariant because 2 conformal metrics g and g = f^g,f{x) ^ 0,x & M, 
generate the same operator. In our case n = 2 and g = rj, so 

*F = -F^,*{dx^ A dx") = - - F^,frrr^e^ro.iisj\ detfjp^ \ dx" A dx^ = 

~F^J-V"v''"e,rapf^/\detr),,\dx" A dx^ = *F 

If we recall the expression (1.7) for the energy-momentum tensor, we'll see 
that the ^-operator is applied there also on 2-forms only, which shows, that 
it is conformally invariant too. Moreover, the expressions (1.7), (1.8) clearly 
show, that the following relations hold: 

Q;{F) = Qli^Fl V.Q^iF) = V.Qli^F). (1.11) 
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These two equalities determine a full equivalence from energetic point of 
view between F and *F. This important fact will be substantially used later 
when we'll be writing down the new equations of Extended Electrodynamics 
(EED). 

Prom pure algebraic point of view the £'M-field, i.e. the 2-form F, has 
two invariants 

h = \f^,uF'''' = - E\ h = ^F^.i^Fr" = 2B.E. (1.12) 
These are the coefficients ( up to a sign) of the two 4-forms 
F A *F = -r){F, F)uj, = -^F^^F'^^oOo, 

F AF = -F A**F^ri{F, *F)uj, = ^F^,{*F)''''cOo. 
The following relations hold 

WQ.^Q^'' = + II (47r)2Q^.Q- = + I'^]5;. 

The equations for the eigen values of F and *F 

det\\F^„ - Xr]i^^\\ = 0, det\\{*F)i^^ - X*r]^„\\ 
look as follows 

+ hX^ - hi = 0, (A*)^ - h{X*f - hi 
For the corresponding eigen values we obtain 




Multiplying the equation F^^'^ = A^^ on the left by —F-^ and adding on the 
two sides of the equation obtained ^Fa/sF^^S'i^, we get 
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-FapF'^^ + A=^ 



This shows, that the cigcn vectors of F arc cigcn vectors of Q too, and the 
eigen values A of F and 7 of Q are related by the above condition. The 
corresponding relation between 7 and A* reads 
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As for the eigen vectors of F, *F and Q we mention just the isotropic case, 
i.e. when Ii = I2 = 0. Clearly, if this is the case, then all eigen values are 
equal to zero and it can be shown, that there exists just one and common 
for F, *F, and Q isotropic eigen direction, defined by the isotropic vector 
C, — 0, and all other eigen vectors are space-like. Moreover, there exist 
two 1-forms A and A*, such that the following presentation holds 



;i.i3) 



< 



Obviously, the 1-forms A and A* are defined up to additive factors cohnear 

to C- We show now that these two 1-forms are spacelike, mutually orthog- 
onal, they have equal magnitudes and are orthogonal to C,: A^ = [A*Y 
0, A.A* = 0, A.C = A*.C = 0. In fact, 

= *{A* A A* A C) = *{A* A *F) = -{A*yF^^dx'' = -{AyA^^i^dx", 



< 47tQI = -Fi,F 
1 



= 7i 



- -(*F)4.(*F)^^ 
1 



2 



{A.C.-A^QiA^C-A^C'') 



-iA.Cr^--i*F),^{*Fr^iA\Cf 



We see that in this case of zero inzvariants Ii = I2 — there exists a new 
invariant, namely the square of A and A*. Besides, from the first row of 
equalities is seen that A is an eigen vector of *F and A* is an eigen vector 
of F. Here are two useful relations: 



det\\F,4 = det\\{*F)^4 = -{hf, det\\{F ± = {hf. (1.14) 



Finally we note that all subdeterminants of third order are proportional to 
E.B, so, the highest order non-zero subdeterminants in this case may be 
those of second order. 
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1.1.3 Continuous media 

All considerations made up to here are characterized by the assumption, that 
in regions with p — the EM- field propagates in space without losing energy 
and momentum, i.e. there is no interaction with the medium, which we call in 
such electromagnetic vacuum. CED extends its applicability to media, 

which interact with the field, exchanging energy-momentum at every point 
and moment. The majority of the known really existing media are built 
of a great number of closely connected, i.e. strongly interacting, neutral 
and electrically charged particles. An exact description of the mechanical 
motion of each one of these particles in presence of the external EM-field is 
practically impossible, so we have to assume simplified models of the various 
media. Doing this, it is important not to forget all conditions and scales 
under which the simplifying assumptions of a given model work. 

In order to adapt the already developed mathematical machinery of the 
theory to describe what happens in presence of such media, which are briefiy 
called macroscopic bodies, the approximation physically small volume is in- 
troduced in the following way. If / denotes the average distance among the 
particles, creating a given medium, if AV denotes the physically small vol- 
ume and if L denotes some typical linear scale of the macroscopic object, we 
want the following relations to hold 

P <^AV <^ L\ (1.15) 

Further we assume these conditions satisfied, and all media, satisfying them 
will be shortly called macromedia. 

From practical point of view important are those macromedia, which can 
be electrified and magnetized when placed in external EM-fields. Such media 
are called dielectrics. The additional electrifying is due to the presence of 
bound charges in these media. Subject to the action of the external field these 
charges perform limited in small regions displacements. Such displacements 
cause appearance of additional charges, currents and dipolc moments. After 
an averaging over the volume AV", they arc denoted respectively by pb-bound 
charge density, jb-bound current density, and P -polarization vector. The ad- 
ditional magnetization is due to the circle-like displacements of the charges, 
generating in this way new magnetic moments. The averaging of these new 
magnetic moments over the volume AV defines the magnetization vector M.. 

In analogy with the case free charges in vacuum the following relations 
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among these new quantities are assumed: 

dP 

Pb = -divP, jb = crotM + (1.16) 

After replacing in Maxwell equations (1.1) and (1.2) j and p by (j + jb) 
and (p + Pb) respectively, and having in view (1-16) we obtain the Maxwell's 
equations for continuous media 

IdD 47r 

- =rotH j, dzvB = 0, (1.17) 

c at c 

1 dB 

-— = -rotE, divD = Anp, (1.18) 

c dt ^ ^ ^ 

where 

H = B-A7rM, D = E + A'kP. (1.19) 

When passing from one medium to another, the dielectric properties of 
which strongly differ from those of the first one, it is naturally to expect a 
violation of the continuous properties of H and D. Therefore it is necessary 
to define the behaviour of these quantities on the corresponding boundary 
surfaces. To this end, two new quantities arc introduced: surface density of 
the electric charge a and surface density of the current i. Then the analysis 
of the above equations brings us to the following relations: 

(D„)2 - (£>„)! = 47ra, {E^)2 - {E^)i = 0, 
47r 

{Hn)2 - {Hn)l = —i, {Bn)2 ' {Bn)l = 0, 

c 

where the index " n" denotes the normal to the boundary surface component 
of the corresponding vector at some point. 

Assuming that the quantity of electromagnetic energy, transformed to 
mechanical work or heat during Isec. in the volume V is equal to Jy{j.E)dV, 
and making use of Maxwell's equations for medium, we get 



47r 



E. 



dt J 



dt 



— div 



—E X H 

An 



(1.20) 



Replacing now D — E + AnP and H — B — AttM in this relation we obtain 



d E'^ + B'' 
"di Stt 



— div 



—E X B 
Air 
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ExM . 

These relations describe the local energy-momentum balance. 

The 8 equations (1.17)-(1.18) have to determine 15 functions Ei,Bi,Hi, 
Di,ji. Clearly, more relations among these functions are needed, in order to 
determine them. The usual additional relations assumed are of the kind 

./ . dE^ ■ dB^ \ ■( ■ dE^ ■ dB^ \ 

V ox'^ ox'^ ) \ ax^ ax*^ J 

The most frequently met assumption is P = P{E), Ai = Ai{B) together 
with the requirement P{0) = 0, A1(0) = 0. A series development gives 

P' = K^jE^ + \k)^E^E^ + M' = a)B^ + \ai^B^B'' + ... 

The tensors K^j^., ... are called polarization tensors (of corresponding rank), 
and a*, a*^, ... are called magnetization tensors (of corresponding rank). For 
and iJ* we obtain respectively 

D'^E' + 47r(4E^' + ^ K^E^E'' + ...) = [S] + 47r4)E^' + ... 

H' = B'- 47r(a}B^' + ^ai^B^B'' + ...) = {6} - AT:a))B^ - ... 

If the medium is homogeneous and isotropic and the EM-field is weak, the 
nonlinearities in these developments are neglected, so 

= (1 + AT^K)5iE^ = eiE^ = eS'jE^ 

and 

W = {1- 4Tia)SiB^ = aiB^ = ad^B^ . 

The constants e and fi = are called dielectric and magnetic perme- 
abilities respectively. In case of nonisotropic media the two tensors £*■ and 
/ij — (q^})~^ are used. 

In the relativistic formulation of the electrodynamics of continuous media 
besides the 2-form F, a new 2-form S is introduced, namely 

S = M3dxAdy-M2dxAdz + MidyAdz-PidxAd^-P2dyAd^-P3dzAd^ 



■iJP ..OB 



-\- c.div 
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as well as a new current 



With these notations the equations (1.16) acquire the following compact form 

- (1.21) 



— — T" 



If we introduce now the 2-form G — F — AnS, then equations (1,17)-(1,18) 
look as follows 

5 * F = 0, SG^ Anj. (1.22) 

The two relations — e^jE^ and — n^H^ may be unified in one relation 
of the kind 



G 



(1.23) 

Obviously, -R^"^ = —R;,^^, R^^ = —R^^^. Comparing now the relativistic 
relation (1.23) and the non- relativistic two relations, we obtain 

Rit = 0, R^ = 0, R,t = el, 

= hirXle''^^ xl = (l^T\ k<l,m<n. 



The equations s^- = sj, /jLj — jil lead to R^^ 
verified that 

The (6 X 6) matrix -R"^^ looks as follows: 



^al3 ■ 



It is immediately 



For the invariant R = R^'^^ wc obtain 



xl 


-xl 


xl 








-xl 


xl 


-xl 








xl 


xl 


x\ 

















e\ 


el el 











e\ 


^2 ^2 











4 


2 3 

^3 ^3 



R^2{e\ + el + el + xl + xl + xl)- 

These algebraic properties of the tensor -Rj^"'' are the same as those of the 
Riemann curvature tensor. Since for vacuum we have si — xi — for i?^^^ 
we get 
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or 



which is exactly the induced by rj metric in the bundle of 2-forms over the 
Minkowski space-time. In other words, the presence of an active continuous 
medium, i.e. non-trivial functions sl{x^) and xli^")-! could be described by 
an appropriate curved metric in the space A^(M). 

Now we are going to consider the problem of energy-momentum distri- 
bution of the field in presence of an active medium. Recall that in case of 
vacuum, these quantities are described by the energy-momentum tensor 



1 

47r 



- F ^F^'^S" - F F'' 



1 

Stt 



-F^aF""- - {*F)^,{*Fy 



The natural generalization of this tensor in presence of a new 2-form S, or 
G, looks as follows 



^« = — 



' IJ.cr-' 



;i.24) 



2 ~ ~ M 

Using the identity, which holds for any two 2-forms in the Minkowski space 
1 



(1.25) 



for W!", is obtained 



" Stt 



-F^aG"" - {*F)^,{*Gy 



1 

Stt 



-G^^F'^'^ - {*G),,{*F) 



Obviously, W^i. = W^^, and if Sfj,^, 0, or equivalently, G = F, we get W^— > 
Q)^- Here are the explicit expressions of VF^ by means of the components of 
the 3- vectors E, B, D, H: 



Stt 
1 



EiDj + EjDi + BiHj + BjHi + 5i{B.H - E.D) 



wt = — 

Stt 



1 



4 --(F.D + M) 



{E X E\ + {Bx D), 
It is easily verified the following relation 

F^,{dGy+G^,idFy+{*F)^,{d*Gy+{*G)^,id*Fy 



1 

Stt 



;i.26) 
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If we require ed, j — in (1.22) the following local conservation law to hold 

V.PF; = (1.27) 

and make use of the above introduced definitions, we get the equation 

AnSi^^SSy = F^,{6Sy - {*F)^,{S * Sy. (1.28) 

This relation determines in what way and how much of the field energy- 
momentum is utilized by the medium, characterized by S^u- Since the equa- 
tions are 8, and the unknown functions are 12, the 4 equations (1.28), al- 
though nonlinear, could be used as additional field equations for F^,^ and 
S^^. In this way the number of the equations becomes equal to the number 
of the unknown functions. 

As we mentioned earlier, the symmetry of the energy-momentum tensor 
is necessary to define correctly integral conserved quantities, using the 
10 Killing vector fields of Minkowski space. The tensor FF^, given by (1.24), 
has the following additional properties: 

W(F,G) ^W(G,F), W(F,G) ^W(*F,*G), 

v.w;{F, G) = v.iy;(G, f), v.w^;(f, g) = v.w;{*f, *g). 

Moreover, since the *-operator is apphed only on 2-forms, the tensor is 
conformally invariant. Finally, if the free current j is not equal to zero, then 
from (1.26) and (1.28) and Maxwell's equations it follows that the energy- 
momentum, transferred over to the medium in an unit 4- volume is determined 
by the well known expression F^yj^ . We should not forget also, that if 
G^, = F^, we obtain = Q^. 

All these properties of Wj^ suggest that it is a good candidate for energy- 
momentum tensor of the EM-field in presence of a continuous medium. As 
for the new condition (1.27), its adequacy to the real processes of energy- 
momentum redistribution in the medium under consideration has to be ver- 
ified every time we are going to use it. 



1.2 Solutions to Maxwell's Equations 
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1.2.1 The solutions as models of physical objects and 
processes 

The equations are mathematical relations among several quantities, some of 
which are known and some are unknown. To solve an equation means to 
find those values of the unknown quantities, which make this equation into 
identity. Clearly, the values found, i.e. the solutions, will depend on the val- 
ues of the known quantities, as well as on the very equation. If the number 
of the independent equations is not equal to the number of the unknowns 
we say that the problem is underdetermined or overdetermined depending on 
whether the number of equations is less or larger than the number of the un- 
knowns respectively. One of the important problems is that of uniqueness, i.e. 
under what conditions a given equation has only one solution. This problem 
is based on the fact that an equation may have many, even infinitely many, 
solutions. For example, the equation f'{x) = is satisfied by every constant 
function; the equation ff = 0, where the unknown function / depends on 
two independent variables {x,y), has infinitely many solutions, which are 
parametrized by one differentiable function of one variable: f{x,y) = g{y)- 

In mathematics we frequently write down relations, which have sense in 
various classes of functions, even in various classes of mathematical objects. 
For example, the equation da = admits as solutions differential forms of 
various degree as well as differential forms with compact or noncompact sup- 
port on various manifolds. This is true for a very large class of differential 
operators, therefore, it is necessary the set of objects, where we are going to 
look for a solution, to be pointed out sufficiently in detail. It may happen, 
that in the set of objects, where we look for a solution, the equation consid- 
ered has no solutions. For example, the widely used Laplace equation A/ = 
has no non-constant solutions in the class of finite functions / : TZ^ — > TZ; 
the widely used wave equation □/ = has infinitely many (l+l)-soliton-like 
solutions, and has no (3-|-l)-soliton-like solutions. 

So, every differential equation separates a class of functions (or objects), 
namely those, which have the local property to satisfy this differential equa- 
tion. Not every solution has properties desired by us, so any separation 
of some subclass of solutions with definite additional properties is made by 
means of imposing additional conditions. 

In theoretical physics we make mathematical models of physical objects 
and processes. We assume the idea, that the really existing objects are finite 
and time-stable. This means, that at any moment of their life they occupy 
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finite and sufficiently small regions of the 3- dimensional space and if there 
are no perturbations from outside, they would live sufficiently long time, i.e. 
they are time-stable. From mathematical point of view this requires that, 
with respect to the space variables {x, y, z) the corresponding mathematical 
objects are everywhere smooth and have compact support. The time evolu- 
tion of the object is usually determined by an equation, defining a definite 
relation among the various derivatives of the components of the correspond- 
ing mathematical object with respect to the time and space coordinates, i.e. 
locally. 

The real physical objects have structural properties, as well as properties 
as a whole, i.e. integral properties. If we arc interested only in the inte- 
gral properties and behaviour of the physical object, we talk about material 
points., and the behaviour of such objects is described by ordinary differential 
equations. The local description of an object, having a dynamical structure, 
requires partial differential equations, as well as rules, pointing out the con- 
nection between the local and integral characteristics and properties of the 
object. 

If we are able to identify a physical object during a finite period of time, 
this means, that with respect to our means of identification this object shows 
definite properties of constancy. If these properties of constancy are measur- 
able, then the introduced by corresponding measurement procedures quanti- 
ties will have constant in time values, so we can talk about conservation laws, 
which are specific for the object considered. There exist physical quantities, 
which can be introduced for various physical objects. The importance of a 
physical quantity is in a direct dependence on its universality, i.e. how much 
broad is the class of objects, admitting this quantity as a characteristic. 

At a definite level and under definite conditions it is possible some objects 
to be transformed into other objects. If such a transformation takes place 
the natural question arises: w/ia^ is the behaviour of those physical quantities, 
which characterize the initial objects, as well as the final ones, is there any 
quantitative connection between the initial and final values of these quanti- 
ties?. A positive response to this question would be of great importance for 
theoretical physics since it would allow some definite mathematical relations 
to be written down immediately and some true conclusions, based on these 
relations and concerning the real objects and processes, to be made. 

The study of the real objects and processes from this point of view has 
resulted into formulation of the so called conservation principles for energy, 
momentum and the full angular momentum. (A principle is any rule which 
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is proved to hold for all known situations and which is extended as a hypoth- 
esis for the unknown situations). The physicists have succeeded in defining 
these quantities in such a way, that the time constancy of their values to 
be an adequate expression of the constancy properties of the real objects 
and to point out definite necessary conditions, which all real processes and 
transformations of objects have to respect and obey. We would like to note 
specially, that the universality of these quantities and principles is of primary 
importance. 

Another very important property of these quantities is their additiveness, 
which enables us to build integral conservation laws, making use of the cor- 
responding local conservation laws, i.e. time independent characteristics of 
the system whole. 

Let us now connect these conclusions with the above mentioned circum- 
stance, that the partial differential equations, which we use to describe the 
existence and evolution of the extended natural objects, admit usually many 
solutions. Most frequently, the model differential equation is a mathemati- 
cal expression of one or several characteristic or typical local features of the 
object (or process), but not of all of them. Therefore, some of the solutions 
may posses properties, which do not correspond to the real integral prop- 
erties of the object under consideration. So, additional conditions, such as 
new equations, inequalities, etc., of local or integral character have to be 
formulated. In such a case, in order to separate the non-adequate solutions, 
we shall always try to combine the above introduced notion for finite and 
time stable object and the availability of conserved quantities for any object, 
observing the following rule: 

In order some solution of a given differential equation to be 
considered as a realistic model of a real object it is necessary this 
solution to be finite, time-stable and the corresponding integral 
energy, momentum and full anguleir momentum to be well defined 
finite quantities. 

1.2.2 Spherically symmetric solutions 

As we mentioned in section 1.1.2 the symmetry of a tangent object (i.e. a 
section of some tensor degree of the tangent and cotangent bundles) is a 
smooth map (usually a diffeomorphism) of the base space onto itself, such 
that the object does not change. We recall that in the non-relativistic frame 
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the base space is TZ^, while in the relativistic frame the base space is 7?.^. 
Since the spherically symmetric solutions of Maxwell's equations are very 
important and lead to the notion of elementary charge, we are going to 
consider this point more in detail in nonrelativistic, as well as, in relativistic 
terms. 

The intuitive notion of spherical symmetry consists in, that at the same 
distance from a given point the properties of the object under consideration 
are the same. The more accurate notion in the frame of the mathematical 
model requires: 

1. Pointing out the point of the base manifold, which is the center of 
symmetry. 

2. A definition of distance. 

A procedure for separation of those maps, which keep the symmetry cen- 
ter stable as well as the distance between any two points unchanged. 

In the frame of non- relativistic considerations the base space TZ^, is fur- 
nished by the Euclcdcan distance Al between any two points with standard 
coordinates {xi,yi,Zi) and (a;2, -^2) 



Al = V(^2 - x^y + (1/2 - yir + iz2 - 

or, equivalently, by the metric tensor g, which in these coordinates has the 
well known canonical components g^ = l,gij = if i ^ j. Then the 
local symmetries^ of g are determined by the relation Lxg — 0. This last 
relation defines a system of differential equations for the components of the 
vector field X. Besides the standard translations [d/dx,d/dy,d/dz], these 
equations have in spherical coordinates (r, 9, (p) the following solutions: 

d d d d d 

Xi = sin </?— cot 9 cos </?— , = - cos </?— + cot 9 sin (/?— , = - — . 

o9 dip o9 dip op) 

These, of course, are the generators of the group 5*0(3). So, we want 
our E'M-field, i.e. the couple {E, B) of vector fields, to be symmetric with 
respect to the local isometries Xi, X2, X-^, i.e. the following equations to hold 

Lx,E = [X„ E] = 0, Lx,B = [X,, B] = 

The solutions of these equations are 

d d 
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Since gu = 1, then the corresponding 1-forms are E = q{r;t)dr, B — 
m{r; t)dr. Obviously, dE = 0, dB = 0, i.e. rotE = 0, rotB = 0. From 
Maxwell's equations with zero current it follows now that E and B do not 
depend on time, i.e. all spherically solutions of the current-free Maxwell's 
equations are static. 

Let's now pay some attention to the following intermediate result: Every 
spherically symmetric, i.e. 5*0(3)- invariant , 1-form a on TZ^ is of the kind 
a — f{r)dr, so it is closed: da = 0. Then, if uj is 5'0(3)-invariant 2-form 
on TZ^, and since the action of the isometries commutes with the *-operator, 
defined by the same metric, we conclude that *u! is also 5'0(3)-invariant, 
hence it is closed. Historically, the things have started with a specially chosen 
1-form. namely the Coulomb force F for a unit test charge, so it necessarily 
it is closed, and the special properties of F come from the observation that 
*F is also closed. Prom modern point of view the S'0(3)-invariant 2-form 
approach seems more natural and clearer. 

What has rest to be done is to solve the two equations d*E — 0, d*B — 0. 
For the *-operator in spherical coordinates we get 

d*E ^d* q{r)dr = d[g(r)r^ sin 9d9 A d^p] = ^ sin 9dr Ad9 Ad(p = 0. 
We obtain q{r) — const /r'^ — q^/r'^, so 

^— B = — — 

These results may be given a topological interpretation. In fact, the 
differential forms 

E — ^dr, *E — qo sin 9d9 A dip 

are defined on the space 71^ \ {0}, i.e out of the point r = 0. Maxwell's 
equations require both of them to be closed, but while the 1-form E — q^/r'^dr 
is exact, the two form *E = qosin9d9 A d(p is not exact, since fg2 *E = go- 
Therefore, the 2-form *E represents the nontrivial cohomology class of the 
space TZ^ \ {0}, moreover, it is the only spherically symmetric representative 
of this class. 

If now we favour the topological character of the field configuration ob- 
tained, we should define the electric charge as the (appropriately paramet- 
rized) cohomology class of the space \ {0}, and the field, as a repre- 
sentative of this class (we chose the spherically symmetric representative 
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uj = qosinOdO A d(p of this class, but this is not obhgatory). The available 
eucledean metric g makes it possible to introduce the *-operator and the 
notion of spherical symmetry. Then the 1-form *uj = qor~^dr, which now 
depends on the metric chosen, is exact and coincides with the notion of a 
field of an elementary source in CED. Note that, if we choose another metric 
g on TZ^ \ {0}, then the relation d * uj = may not hold. So, the notion 
of field in CED is strongly connected with the eucledean metric, while the 
notion of charge has topological origin. 

This picture favours the topological aspect of the electric charge and gives 
some benefits for the relativistic formulation of CED. In fact, the closed 
2-form u is pulled back naturally as a closed 2-form on (7^^ \ {0}) x TZ, 
and the pscudocuclcdean metric r] guarantees that *^cj is closed, even exact. 
Of course, the consideration of all solutions of the two equations do; = 0, 
d * a; = as admissible models of real £'M-fields is an additional hypothesis, 
but it is a natural one in the frame of this approach. 

In the frame of the relativistic formalism the spherically symmetrical 
problem requires to solve the equations dF = 0, d * F = together with the 
symmetry relations Lx^F = 0, Ljf. *F = 0, which reduce to the first relation 
only since the *— operator is spherically symmetric. The most general form 
of a spherically symmetric F looks as follows: 

F = f{r, i)dr A d^ + h{r, sin OdO A d(p. 

Maxwell's equations require 

C 

F = -^dr AdC + C2 sin 9d9 A dip, 
where Ci and C2 are constants. The additional requirements 

determine Ci = go/47r, C2 = 0. 

1.2.3 Local and integral interaction energy of spheri- 
cally symmetric fields. The Coulomb's force 

Let us consider the following situation: two electric charges, occupying two 3- 
dimensional regions {Oi} and {O2}, considered as "open balls". The bound- 
aries of these balls are the two non-overlapping spheres 5*1 and ^2, and the 
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distance between the centers of these two spheres is R. The non-trivial 
topology of the space, where the two representatives Ui and uj2 of the corre- 
sponding cohomology classes are defined is TZ^ \ {Oi, O2}. Our purpose is to 
define interaction between the two fields cui and 002 in such a way, that the 
calculated integral interaction energy to be equal to the well known classical 
expression W = qiq2/r. Clearly, the local interaction energy expression has 
to be a 3-form and symmetric with respect to the two interacting fields. We 
define this 3-form as follows: 

w = -^cui A *u}2. (1-29) 
47r 

The two fields uji and uj2 are the only spherically symmetric representatives 
of the two cohomology classes. Since the eucledean metric is also spherically 
symmetric, the corresponding 1-forms *u!i and *u!2 are also closed 1-forms 
and even exact. Therefore, there are two functions /i and /2, such that 
d/i = *u!i and d/2 = *u!2- The 3-form w turns out to be exact, in fact 

1 1 

4:7rw — uji A *UJ2 = -cJi A *UJ2 + -0J2 A *UJ\ — 

= ^[a;i A d/2 + a;2 A d/i] = ^[d(/2a;i) + d(/ia;2)] = \^^\h^\ + h'^-^- 

We integrate now this expression over the region D out of the two open balls. 
Obviously, the boundary Sd of D \^ So = S"^ U SiU 82- Then, making use 
of the Stokes theorem, we obtain 

W ^ [ w = ^ f d[/2u;i + fiu}2] = 
Jd 871 Jd 



1 

'Stt 



We introduce two spherical coordinate systems (r, 9, (p) and (r, 9, (f) originat- 
ing at the centers of Si and 5*2. We can write 

uui = qi sin 9d9 A dip, fi — — — , 

r 

- - q2 
1^2 = q2 sin 9d9 A dip, /2 = — — 

r 
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Since /i and /2 get zero values on the infinite sphere, we have to calculate 
the integrals over the two spheres and S2. Having in view the relations 



/ a;i = 0, / UJ2 = 0, fi\si = const, f2\s2 = const 

J S2 ^ Si 



we obtain 
1 



Stt 



/ f2(^l+ / = ^ / 

JSi JS2 J OTT Us 



sin 9d9 Adip ^ f sin 9d9 A dip 



iSi r 



i 



S2 r 



Q1Q2 
Stt 



R\ Jsi f R2 JS2 r 



dSo 



where dSi — R\ sin 9d9 A dip and dS2 — R^ sin 9d9 A dip are the two surface 
elements. Since the function 1/r is harmonic out of the point r — 0, then 
using the mean- value theorem for a harmonic function, we get 



AnRj Is 



dSi 

Si f 



Finally, 



1 

R' 



W = 



— / 

kRo J So 



dS2 
A-kR^ Js2 r 



1 

R' 



11(12 

R ' 



in correspondence with the purpose we set. This result shows that the 
differential dW , where W is considered as a function of the two centers 
Ti = 0,r2 = 0, has nothing to do with *u!i, or *u!2, which are not defined at 
these points at all. 

In order to introduce the Coulomb's force in a correct way we first note 
that its real sense is to define what part of the full energy and momentum of 
the two fields is transformed into mechanical energy and momentum of the 
two particles as a consequence of the local interaction of the two fields. This 
is an integral effect, since the very concept of electric charge has an integral 
character. Formally, we attack this problem in the following way. First we 
consider the trivial bundle {TZ^ x qiS'^, tt, TZ^, S'^). In the natural coordinates 
{x, y, z; 9, ip) we consider the 3-form 



^^12 = qisin9d9 A dip A 7r*(*a;2) = qisin9d9 A dip A 



q2{xdx + ydy + zdz) 
^(x2 + y^ + ^2)3 
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defined on the whole bundle space. Now the Coulomb's force can be defined 
as the fiber integral of flu as follows: 



(xdx + ydy + zdz) 
ili2 = qm- 



3 



The so obtained Coulomb's force could be interpreted as a characteristic of 
the field UJ2 when qi = 1, but it shouldn't be identified in any way with *u!2, 
although they look the same. 

Let's consider briefly the problem of general spherically symmetric solu- 
tions of Maxwell's equations. Since in this case the time derivatives vanish 
we have the system 

rotE = 0, rotB = 0, divE = 0, divB = 0, 

or, all the same 

dE = 0, dB = 0, 5E = 0, SB ^ 0. 
Now, the Laplace operator A is defined by 

A = d5 + 5d, 



so, we get 

AE = 0, AS = 0. 

(We should not forget that these are just necessary conditions, i.e. the last 
two equations may have solutions, which do not satisfy the static Maxwell 
equations.) We see, that the components of E and B are harmonic func- 
tions. According to the theory of these functions, they can not have local 
exteremums inside the regions of harmonicity, and if the region is the whole 
space 71^ and the function is in addition bounded, it is a constant. These two 
properties of the solutions of Laplace equation do not recommend them as 
possible models of free real objects, which are finite and have to be described 
by finite functions of the three space variables {x,y, z), i.e. such, that neces- 
sarily to achieve their maximum values, since they arc zero out of some finite 
subregions of TZ'^. Only in some topologically non-trivial regions they could 
be of some interest in this respect. 

Finally we note, that from relativistic point of view the concept "static 
field" is not invariant, since the Lorentz transformations transform any "static 
field" into "non-static" one. 
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1.2.4 Wave equations and wave solutions 

Prom physical point of view when we talk about waves we mean propaga- 
tion of some disturbance, or perturbation, in a given medium. It is also 
assumed that the perturbation does not alter the characteristic properties 
of the medium, and the time-evolution of the perturbation depends on the 
medium properties as well as on the specificities of the very perturbation. 
The waves are divided to 2 classes: elementary (linear) and intrinsically co- 
ordinated (nonlinear). The elementary waves are observed in homogeneous 
media and are generated by perturbing the equilibrium state of the medium 
through small quantities of external energy and momentum. The important 
properties of linear waves come from the condition, that during the propa- 
gation of the initial disturbance throughout the medium the structure of the 
medium does not change irreversibly, and the various such propagating per- 
turbations do not interact with each other substantially. Prom mathematical 
point of view this means that the evolutionary equations, which are partial 
differential equations, describing such phenomena, are linear, so any linear 
combination with constant coefficients of solutions gives again a solution. In 
other words, the set of solutions of such equations is a real (finite or infinite 
dimensional) vector space. 

The intrinsically coordinated, or nonlinear, waves disturb more deeply the 
medium structure, but the corresponding changes of the medium structure 
stay reversible. When subject to several such perturbations, the medium re- 
sponses to the various disturbances is different in general, so the medium 
reorganization requires more complicated intrinsic coordination. All this 
demonstrates itself in various ways, depending on the medium properties 
and the initial perturbation. What we observe from outside is, that some 
important properties of the initial perturbations are changed in result of the 
interaction. In some cases we observe a time-stable coordination among the 
responding reactions of the medium and if the corresponding formation is 
finite, we may consider it as a new object. If this object keeps its energy and 
momentum we frequently call the corresponding medium vacuum. Clearly, 
such objects can exist only in appropriate media. In such cases, studying the 
objects, we get some information about the medium itself. From mathemat- 
ical point of view these waves are described by nonlinear equations, so that 
a linear combination of solutions is not, as a rule, a new solution. The huge 
variety of various such cases could hardly be looked at from a single point of 
view, except some most general features. 
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It is important to note, that in the both cases, hnear and nonhnear, the 
perturbations are bearable for the medium in the sense, that they do not 
destroy it. We are not going to consider here unbearable perturbations. 

One common for every kind of waves characteristic is the polarization. 
The polarization determines the relation between the direction of propaga- 
tion (at some point of the medium) and the direction of deviation from the 
equilibrium state of the medium point considered. If these two directions 
are parallel we say that the polarization is longitudinal, and if these direc- 
tions are orthogonal we say the polarization is transverse. In general the 
polarization depends on the space-time point, i.e. it is a local characteris- 
tic. When the wave passes through some region of the medium, the points 
inside this medium commit some displacements along some (usually closed) 
trajectories. If these trajectories are straight lines we say that the polariza- 
tion is linear, if they are circles we say the polarization is circular, etc. It is 
important to note that the polarization is an intrinsic property of the wave, 
therefore it ia an important for the theory characteristic. In particular, the 
mathematical character of the object (scalar, tensor, spinor, differential form, 
etc.), describing the wave, depends substantially on it. If the wave is linear, 
and the corresponding equation admits solutions with various polarizations, 
then summing up solutions with appropriate polarizations we can obtain a 
solution with a beforehand defined polarization. 

Other common characteristics of the waves are the propagation velocity, 
determining the energy transfer from point to point of the medium, and the 
phase surface, built of all points, being in the same state with respect to the 
equilibrium state at a given moment. We are going to make use of these 
characteristics further. 

Let us consider now the Maxwell's equations in regions far from sources: 
p — 0. We have 

1 OE 1 OB 

-_ = rotB, divB = 0, -— = -rotE, divE = 0. 

COT COT 

From the first and the second equations we obtain 

15 1 0^ B 1 0^ B 

rot(rotB) - -—rotE = grad(divB) - AB + -^^—r- = - AB + r = 0. 
cot c^ at^ c^ ot^ 

Prom the third and the fourth equations we obtain 

1 d^E 
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These are the well known wave equations, and we are going to consider some 
of their properties. First we note, that these equations are just necessary 
conditions for every solution of the vacuum Maxwell's equations. Therefore, 
they may have solutions, which do not satisfy Maxwell's equations. Second, 
every component of E and B satisfies the same equation and does not depend 
on the rest ones. Third, these are second order differential equations of 
hyperbolic type. 

We are interested in the following: Do the vacuum Maxwell's equations 
admit finite and time-stable solutions, so that such solutions to serve as mod- 
els of real objects?. The positive answer to this question would be a serious 
virtue from the point of view of their adequacy as model equations for an 
important class of real objects, while the negative answer would make us 
searching for new equations, having solutions with the desired properties. At 
the beginning of the last century (about 1818), i.e. more than 40 years be- 
fore the appearance of Maxwell's equations this problem has been essentially 
solved by Poisson, and because of its importance we shall consider it in some 
more detail. 

Let's denote by u any component of the vector fields E and B. Then u 
satisfies the wave equation. We are interested in the behaviour oi u at t > 0, 
ii at t — the function u satisfies the initial conditions 



Further we assume that the functions f{x,y,z) and F{x,y,z) are finite, i.e. 
they are different from zero in some finite region D C TZ^, which corresponds 
to the above introduced concept of a real object. Besides, we assume also 
that / is continuously differentiable up to third order, and F is continuously 
differentiable up to the second order. Under these conditions Poisson proved 
that a unique solution u{x,y, z;t) of the wave equation is defined, and it is 
expressed by the initial conditions / and F by the following formula: 



where P is a point on the sphere S centered at the point {x, y, z) and a radius 

r = at, da,, is the surface element on Sr=ct- 

The above formula (1.30) shows the following. In order to get the solution 
at the point {x, y, z), being at an arbitrary position with respect to the region 



u\t=o = f{x,y,z), 



du 

dt t=o 



F{x,y,z). 
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D, where the initial condition, defined by the two functions / and F, is 
concentrated, it is necessary and sufficient to integrate these initial conditions 
over a sphere S, centered at {x, y, z) and having a radius of r = ct. Clearly, 
the solution will be different from zero only if the sphere Sr=ct crosses the 
region D at the moment i > 0. Consequently, if ri = cti is the shortest 
distance from (a;, z) to and r2 = ct2 is the longest distance from (x, y, z) 
to -D, then the solution will be different from zero only inside the interval 
(^1,^2). 

Prom another point of view this means the following. The initially con- 
centrated perturbation in the region D begins to expand radially , it comes to 
the point {x,y,z) at the moment t > 0, makes it "vibrate" ( i.e. our devices 
show the availability of a field) during the time interval At = t2 — ti, after 
this the point goes back to its initial condition and our devices find no more 
the field. Through every point out of D there will pass a wave, and its fore- 
front reaches the point {x, y, z) at the moment ti while its backfront leaves 
the same point at the moment t2. Roughly speaking, the initial condition 
"blows up" radially and goes to infinity with the velocity of light. 

This mathematical result shows that every finite nonstatic solution of 
Maxwell's equations in vacuum is time-unstable, so these equations have no 
smooth enough time-dependent solutions, which could be used as models of 
real objects. As for the static solutions, as it was mentioned earlier, they 
also can not describe real objects. 

These explicit results state clearly, that if we want to describe 3-dimensio- 
nal time-dependent soliton-like electromagnetic formations (or configurations) 
it is necessary to leave off Maxwell's equations and to look for new equations 
for E and B, or for F^^^. 

In relativistic notations the vacuum Maxwell's equations 

dF = 0, (5F = 

naturally, require F to satisfy the equation 

(d(5 + 5d)F = AF = 0, 

which in standard coordinates {x^ , x'^ , x^ , x"^) = {x,y,z,^ = ct) gives the 
usual wave equations for the components F^j^: 

^ dx'^dxf^ dt"^ dx"^ dy^ dz^ 

Of course, not every solution of AF = is a solution to Maxwell's equations. 
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1.2.5 Plane electromagnetic waves 



The exact solutions of Maxwell's equations in the whole 3-space, known as 
plane electromagnetic waves, are interesting not as some models of really ex- 
isting objects, but as an convenient way to introduce some important charac- 
teristics of a class of EM-fields. The standard (i.e. the most widely spread) 
way to define such a solution is the following: there is a rectangular coordinate 
system {x,y,z), in which this solution depends on one space variable only. 
The time-dependence of the solution is determined by the equations. Right 
now we note, that such a solution, if it exists, will be infinite] In fact, if z is 
the only coordinate, on which the solution depends, even if the dependence 
on z is finite, i.e. localized and without singularities, with respect to the rest 
two coordinates {x, y) this solution is constant; even at x — > oo, y — > oo the 
values of the components of E and B, or F^^,, do not change. This simply 
means that the initial condition occupies the whole 3-space, or its infinite 
subregion, with finite values for the components of E and B. Since the in- 
tegral energy W of every solution does not depend on time, we calculate it, 
making use of the initial condition, and obtain 



W 



^ ^ f {E^ + B'^)dxdydz = oo. 



(1.31) 



Let us now see how the plane wave looks like in the corresponding coordinate 
system, where E and B depend on z and t only. Since the derivatives with 
respect to x and y will be zero, from the wave equations we get 



E 



Ei{ct + ez), E2{ct + ez), E^{ct ez) 



B 



±1. 



Bi{ct -\- ez),B2{ct -\- ez),B^{ct -\- ez) , 

Now the equations divE = and divB = require E^ = const and S3 = 
const. Let us put these constants equal to zero since we do not interest in 
constant solutions. The Maxwell's equations reduce to 

IdBi dE2 ldE2 dBi 



c dt 
Id El 



dz ' c dt 
dB2 1 dB2 



c dt dz ^ c dt 

These equations have the following solution: 



dz ' 
dEi 
dz 



E 



Ei{ct + ez), E2{ct + ez),0 = u{ct + ez),p{ct + ez),0 
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B 



B-i{ct -\- ez),B2{ct -\- ez),Q — ep{ct + ez) , —eu{ct + ez),0 



For the Poynting's vector we obtain AnS — 0, 0, —£c{u'^+p'^) . This solution, 
obviously, has the properties 

E.B = 0, -B^ = 0, 

i.e. the field has zero invariants. Now wc show in relativistic terms how 
the requirement for zero invariants Ii = I2 = determines the solution plane 
wave. According to subsection 1.1.2 at zero invariants the energy-momentum 
tensor Q^u, defined by (1.7), has only zero eigen values and unique isotropic 
eigen direction, defined by the couple of opposite isotropic vectors ±V — eV. 
Making use of the representation (1.13) for F and *F we obtain for Q^i, 

Ql = -A^A'^V^V' = -{A*\{A*YV^Vr 

On the solutions of Maxwell's equations we shall have 

- v,q; = -A^vv^v^ - v^wM^y")- 

This relation shows, that the integral lines of the vector field V are isitropic 
geodesies, i.e. straight lines. Let's now choose the coordinates {x,y,z,^) in 
such way that the integral lines of V to lie entirely in the plane {z,^)- Since 
V4 always, we can suppose V4 — 1. Then in these coordinates we shall 
have V — (0, 0, £, 1) and 

F12 = F34 = 0, Fi3 = eFi4, F23 = eF24, 

A = (Fi4, F24, 0, 0), A* = (-F23, Fi3, 0, 0) = {-eA2, sAi, 0, 0). 

Clearly, in these notations A and A* are the relativistic equivalents of E and 

B respectively. 

Denoting F14 = u, F24 = p, for dF = and 6F = we get 

dF — e{px — Uy)dx A dy A dz + {px — Uy)dx Ady A d^+ 

e{u(^ — euz)dx A dz A d^ + e{p^ — epz)dy A dz A dC, = 0, 

5F = {u^ - euz)dx + (p^ - epz)dy + e{ux + Py)dz + {u^ + Py)di = 0. 

From these equations it follows + Py = and Uy — p^ = and from these 
last relations we get the equations u^x + Uyy = 0, Pxx + Pyy = 0, i.e. u 
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and p are harmonic functions with respect to the variables (a:, y) . Since we 
are searching for finite solutions in the whole space, from the well known 
properties of the harmonic functions it follows that u and p do not depend 
on {x, y, ). Thus, 

u^u{^ + ez), p = p{^ + ez). 

It is seen that the dependence of the field components on the unique space- 
variable z stays arbitrary, i.e. can not he determined by the Maxwell's equa- 
tions. Obviously, if {k, Lpo) and (/, tpo) arc two couples of real numbers, 
then u{k^ + ekz + (po) and p{l^ + elz + ipo) define again a solution. On 
the other hand the numbers k, I define two constant isotropic vectors k = 
(0, 0, £/c3, /C4), 1 = (0, 0, £^3, ^4), k^ — ki — k, I3 — U — I, so we can write 
down u{k$, + ekz + ipo) = u{k^x^^ + v?o); vi^i + ^^-^ + "^o) = pilfj.x'^ + V'o)- 
The vectors (k. A;), (1, /), or just their space-like parts k and 1 are called wave 
vectors of the two independent solutions 

F{u) = eudx /\dz + udx A d^, F{p) — epdy Adz + pdy A dC,, 

and the quantities ^po,ipo are called phase constants. The quantities {k^x^^ + 
(po), {l^x^^+ipo) are called phases. Every of the two solutions is called linearly 
polarized plane wave. 

It seems important to note that the general plane wave is a sum, or a 
linear combination, of two entirely independent ( except the common direction 
of propagation) linearly polarized plane waves. Because of the linearity of 
Maxwell's equations a sum of linearly polarized plane waves with different 
directions of propagation is also a solution, but it is no more a plane wave. To 
propagate as a whole along a definite direction is one of the specific properties 
of the plane waves but this is not enough to use it as a model solution for 
real objects because of their infinity: infinite volume, infinite energy, infinite 
momentum and angular momentum. 

A special interest for the theory is the choice of the two functions u and 
p as elementary periodic functions, namely 

u = Uocos{k^x^ + v?o), P = PgCos{1^x^ + ^0), 

since most of the real EM-fields show definite properties of periodicity. In 
such case the quantities ki/c and /4/c are usually called frequency u, (or 
circular frequency cu — 27ri/), and the quantity A = |k|~^ is called wave 
length. So we can write down (e.g. for u) 

^ = ^ = - = <^^' 

k 1/ 
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where T = 1/v is called period. Such waves are usually called harmonic. 
We'd like to note specially, that this time-periodicity is a consequence of 
the specially chosen initial condition, namely that the function u, considered 
as a function of one variable, is a periodic function of the space-variable z. 
This periodicity is admissible by the Maxwell's equations, but it is not a 
necessary consequence of these equations. So, the introduced "wave" char- 
acteristics come from a special class of initial conditions and nothing more. 
Finally we note that the general plane wave is determined by 4 real parame- 
ters /c, /, ipQ, tpo and 2 arbitrary functions of one (and the same) independent 
variable. 

Two linearly polarized harmonic plane waves of the kind 



Uocos{i't + £- + (po), 0, 
A 



0, -eUocos{i't + e- + (po), 
A 



z 



0, Posin(ut + £- + ipo), 
A 



Bo = 



z 



£Posin(i't -|- £-- -|- ipo), 0, 
A 



are called consistent. Summing them up we obtain again a harmonic plane 
wave 



E^ 



z z 
Uocoslvt + e- + (fio), Posinlut + e- + (fio), 
A A 



B 



z z 
ePQsin{i't + e- + ipo), -eUQCos{ut + e- + ipo), 
A A 



We see that the consistent linearly cycling E and B of the two harmonic and 
linearly polarized waves create an illusion for circulating E and B of their 
superposition, i.e. the couple of orthogonal vectors {E, B) takes part in two 
motions: advancing along z and circulating in the plane (x, y) left-wise or 
right-wise with the frequency of v. In this way we get an impression about 
an elliptically (or circularly at Uq = Pq) polarized plane wave. 

Note that the harmonic plane wave occupies the whole 3-space, its energy- 
density is the constant quantity w ~ (t/^ + P^), its full energy is infinite since 
the 3-dimensional integral of a constant over the whole space TZ^ is infinity. 

Finally we note the transverse character of any plane wave, which is seen 
from the transverse direction of E and B with respect to the direction of 
propagation. 

1.2.6 Potentials 



According to the general notion of potential this is a scalar or some tensor 
field, from which by differentiating (once or more) the physical field, i.e. the 
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field that is physically measured, is obtained. The most frequent case is a 
"single" differentiating, but there are examples of "double" differentiating. 
Such is the case of gravitational field in the frame of General Relativity, where 
the curvature field, which is identified with the physical field, is obtained 
from the metric tensor through double differentiating. As we shall see now, 
some of the solutions of Maxwell's equations can also be obtained by double 
differentiating of some of the solutions of the wave equation. This is the 
method of Hertz potentials. In fact, if Ai and A2 are two solutions (vector 
fields) of the vector wave equation DA — 0, then the expressions 

d d 
E = rot{rotAi) - —rotA2, B = rot{rotA2) + ^'^o^^i 

define a solution to Maxwell's equations. In relativistic notations we have 

{UG)^^ = - [{d5 + 5d)G\^^ = 0. 

Then 

F^dSG^ -5dG 

is a solution to Maxwell's equations 6F = 0, dF = 0. In fact, since dod = 0, 
SoS^O and Dod^doD, Do^^^on, then 



SF = S{dSG) = (□ - dS)SG = SaG = 0, dF = d{dSG) = 0. 

These solutions are used for description of the Hertz vibrator's radiation. 
The spherical wave 

ai{ct — r) a2{ct — r) as{ct — r) 
) ) ) 

rjT' ry 

which is a solution of the wave equation, but is not a solution to the Maxwell's 
equations, is used in the above shown way to build a solution of Maxwell's 
equations. Clearly, the standard choice of the components Oj as elementary 
harmonic functions brings a singularity at the point x = y = z = 

The general solution of Maxwell's equations with non-zero and not de- 
pending on the field electric current, is obtained by means of introducing 
the 4-potential, i.e. an 1-form A = A^dx^ and defining F = dA, which is 
always possible, since in CED we have in all cases dF = 0. The 1-form A is 
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determined up to a term of the kind df, since A and A + df give the same 
F. Then 

5F = 5dA = (□ - d5)A - n{A + df) - d5{A + df). 

We see that choosing / as a solution to the equation — 5d/ = □/ = SA we can 
redefine A as A' = A + df. Clearly A and A' define the same F, and 6 A' = 0. 
So, the 1-form A' satisfies the equation SdA' — OA' — Anj, the solutions of 
which are well known when the current j does not depend on A'. Of course, 
when the change in the mechanical energy of the charge- carriers is taken into 
account, and these carriers do not "appear" and "disappear", then according 
to subsection (1.1.2), equation (1.9), the 4-current = pu^ depends on the 
field F, and the equations become nonhnear. Therefore, it is an illusion to 
think, that the 4-potential approach solves the problem completely: in all 
really significant cases the 4-current depends substantially on the field F 
in accordance with equation (1.9), which, in turn, follows from the energy- 
momentum conservation law and could hardly be put into some doubt. This 
fact shows some inadequacy of the thesis for the universal character of the 
4-potential approach, its (always possible) introduction does not lead to a 
complete solution of the entire problem. In all cases the energy-momentum 
conservation requires nonlinear inter-relations between the current and the 
field. 

1.3 Amplitude and Phase 

1.3.1 Amplitude and phase of a plane wave 

The importance of the concepts of amplitude and phase in the electromagnetic 
theory is out of any doubt, but sufficiently general and universal definitions 
of these concepts in CED are still missing. Our purpose in this section is 
to consider some new ways to introduce these concepts into theory through 
a pure algebraic and coordinate free approach in both, nonrelativistic and 
relativistic formalisms. We consider first the case of a plane wave solution. 
In the corresponding coordinate system the plane wave solution is 

F — eUoCos{ki^x'^ + (fo)dx Adz + Uocos{knx'^ + (fio)dx A d^+ 

ePQsin{lfj_x^ + ipo)dy f\dz + Posin{lnX^ + ipQ)dy A d^. 
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or, in terms of E and B 



E 



B 



Uocos{k^x^' + (^o), Posin{l^x^ + ^o), 



ePQsin{lij_x^ + V^o); ~£ UoCos{knX^ + ipo), 



As we mentioned, the quantity k/^x'^ + fo — vt + ez / X + (po is called phase 
of the plane wave. As for the amplitude, according to the usual sense of this 
concept, it is the magnitude, or the maximum value of a given quantity. In 
our case we have a couple {E, B) of vector fields, so it seems natural to define 
the amplitude by the relation 



i.e. a square root of the energy-density. For the vector product E x B we 
obtain 

ExB= [(0, 0, -e{U^ + Po') 

This is a constant vector. 

Now, the triple {E, B,ExB) defines a basis of the tangent (or cotangent) 
space at every point, where the field is different from zero (we assume E ^ 
0,B 0). Moreover this is an orthogonal basis. We denote this basis by TZi, 
so we can write TZi = {E, —eB, —eE x B). From the properties of the plane 
wave solutions we obtain \E\ = \B\. But, the physical dimension of the third 
vector E X B is different from that of the first two. So, we introduce the 
factor a 
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Making use of a, we introduce the basis 

7^ = \aE, -eaB, -eo?E x B . 

\ ience, at every point we've got two bases: and the coordinate basis — 
dy, d}^ , as well as the corresponding co-bases TV and 7?.q = [dx, dy, dz). 
We are interesting in the invariants of the corresponding transformation ma- 
trix M from TZl to TZ* . It is defined by the relation TZ^.M = TZ*. So, we 
obtain 

aEi -saBi -£a^{E x B)i 

M. — aE2 —eaB2 —ea'^{E x B)2 
aEs -saBs -ea'^{E x B)^ 
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We shall try to express the amplitude and the phase of the plane wave as 
functions of the invariants of this matrix. So, in all cases , where this is 
possible, the invariant character of the so defined phase and amplitude will 
be out of doubt. As it is well known, in general, every square n x n-matrix 
jC has n invariants Ii, I2, In, where Ik is the sum of all principle minors of 
order k. The invariant Ii{C) = Cn + ... + C^n is the sum of all elements on 
the principle diagonal, and the invariant /„ = det{C) is the determinant of 
the matrix. In our case n = 3, so for the invariant I2 we get 



det 


mil 


mi2 


+ det 


mil 


"713 


+ det 


"722 


"723 




m2i 


"7-22 




"731 


"733 




"732 


"733 



We compute h{TZ). 

I^ijl) = det 



aEi —eaBi —ea'^{E x B)i 
aE2 —eaB2 —£a^{E x 3)2 
aEs -eaBs —ea^{E x 5)3 



a\E X Bf. 



Now the amplitude A of the plane electromagnetic wave may be defined as 
follows: 

A = ^a-^h{n) = ^a^{E x Bf = a\E x B\. 



1.3.2 Amplitude and phase of a general field 



The invariant character of the above given definition of the plane wave am- 
plitude suggests its natural extending to an arbitrary field. So, if the cou- 
ple {E,B) represents the field, we introduce the matrix M.{TZ) of the basis 
aE, —aB, —a^{E x B) and define the amplitude M of the field by 



7^ = 



A{E, B) = ^a-'^h{n) = a\E x B\. 



(1.32) 



We go further now to define the phase in the general case. We'll make 
use of the matrix of the basis 



n 



aE,-aB,-a'^E x B 



defined by the general field {E, B). The invariants 

h{n) = aEi - aB2 - a^{E x B)^, 
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I^(n) = -a^{E X B% + a^[E{E.B)-B{E.E)\^ + a" 

h{n) = a\E X B)\ 
obviously, are physically dimensionless. When the inequality 



E{B.B) - B{E.B) 



h{n) - 1 



< 1, 



holds, then the function arccos is defined on the expression on the left. In 
these cases, by definition, the phase cp of the field {E, B) shall be defined by 



(/? = arccos 



h{n) 



(1.33) 



For the plane wave solution 
E 



we get 



u{ct + £z),p{ct + £z), 

/l(7^) = /2(7^) 



B 



2u 



p{ct + ez), —u{ct + ez), 
2Ei 



+ 1 



1^1 



+ 1, 



and for a circularly polarized plane monochromatic wave we get ip = k^x^ + 
const. 

Let's now see when the basis TZ is normed, i.e. when 

\aE\ = 1, \aB\ = 1, a^\E x B\ = 1. 

From the first two equations it follows \E\ = \B\, and from the third equation 
it follows E.B = 0. Moreover, the relations \E\ = \B\, E.B = follow from 
the third equation only: a'^\E x B\ = 1. So, the normed character of TZ 
leads to its orthonormal character, consequently, detM.{7l) — 1. Vice versa, 
the requirement detA4(7?.) = 1 leads to the orthonormal character of TZ. 
We obtain, that the requirement det (TV) = 1 is equivalent to the null 
character of the field: = B^ - E^ = 0, E.B = 0. 

The relations obtained suggest to define and consider the following 4- 
linear map: 

xi yi {v X w)i 



R{x, y, V, w) — det 



X2 2/2 {v X w)2 
X3 ys {v X w)3 
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The following relations hold: 

R{x,y,v,w) = {x X y).{v x w) = y x [v x w) 



.X 



[v X w) X X 



R{x, y, V, w) = -R{y, x, v, w), 

R{x, y, V, w) = -R{x, y, w, v), 

R{x, y, V, w) + R{x, V, w, y) + R{x, w, y, v) = 0, 

R{x,y,v,w) = R{v,w,x,y), 

R{x,y,x,y) = {x x yf. 

We note that this 4-linear map has all algebraic properties of the Rieman- 
nian curvature tensor, therefore in the frame of this section, we shall call it 
algebraic curvature. For the corresponding 2-dimensional curvature K{x,y), 
determined by the two vectors {x, y) we obtain 



K{x.y) = 



R{x,y,x,y) ^ _ 

x^y"^ — [x.yY x'^y'^{l — cos'^{x,y)) x'^y'^sin'^{x,y) 



{x X yf 



x^y'^sin'^{x,y) ^ 



Let (61,62,63) be a basis. We compute the corresponding Ricci tensor Rik 
and the scalar curvature R. 

Rijki = R{ei, Cj, Ck, ei) = {si X ej).{ek x ci), 

Rik = ^R\kl = (^i X ei).(6fc X 61) + {Ci X 62). {Ck X 62) + {Ci X 63).(6fc X 63), 

I 

R = E^i = 2[(ei X 62)=^ + (61 X 63)' + (62 X 63)' 

i 

For our basis TZi we obtain the following non-zero components: 



R 



12,12 



R 



13,13 



R 



2^2 



sin\E,B), 



2B' 



12,12- 



^2 ^2' -^23,23 -Rl2,12-^2 _^ ^2' 



and for the scalar curvature we get 



^iE,B) = 24^^^^^s^n\E,B) 
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After this short retreat let's go back to the quantities phase and amplitude. 
The above mathematical consideration suggests to try to relate these two 
concepts with the notion of curvature in purely formal sense, namely as a 
2-form with values in the bundle I/r(7e3) of linear maps in the tangent bundle. 
Most generally, a 2-form it! with values in the bundle I/t(-r,3) looks as follows 



R 



^Ryx' A dx^ ® ^ 



dx^. 



We have to determine the coefficients Rfji, i.e. we have to define a (3 x 3)- 
matrix TZ of 2-forms. We define this matrix in the following way: 



aEidy A dz aBidy A dz a'^{E x B)idy A dz 
aE2dz A dx aB2dz A dx a^{E x B)2dz A dx 
aE^dxAdy aB^dxAdy a^{E x B)^dx Ady 



The columns of this matrix are the 2-forms *B, *[E x B), multiplied 
by the factor a at some degree in order to obtain physically dimensionless 
quantities. 

We aim to define the amplitude and the phase of the field {E, B), making 
use of this matrix. The amphtude TZ of the field we define by 



«4 — —RijkiR^^'^'' — — 
6a 6a 



1 + 2' 



{E X Bf 



(1.34) 



(E2 + B2)2_ 

In order to define the phase we first consider the 2-form tr o R. We get 

tr oTZ = aEidy Adz + aB2dz Adx + a^{E x B)y,dx A dy. 
The square of this 2-form is 



{tr o nf = a^ 



{E^f + {B2f + a\ExB)^\ 



Now the phase </? of the field we define by 



</? = arccos 





{tr oRf-l 











;i.35) 



whenever the right-hand expression is well defined. 
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Let's compute this quantity for a plane wave, moving along the 2;-coor- 
dinate from — oo to +00: 



E 



u{z — ct), p{z — ct), 



B 



—p{z — ct), u{z — ct), 



We get 



troR^ 



u 



dyAdz+ 



0, 0, {u'+p') 



1 



a — 



u 



dzAdx-\-dx/\dy, {troRf = 1 + 



2u' 



+p^ 



A 



+ p^, if — arccos 



\u\ 



In the case of a plane harmonic wave = k^^x^ + </7o- Note that since in this 
coordinate system the components of the plane wave do not depend on the 
coordinates x and |/, the corresponding 2-form troRis closed: d{tr oR) = 0. 

If we work in relativistic terms, it is necessary to introduce some changes 
in the matrix of 2-forms. First, we add one more column and one more row. 
Second, instead of *E, *B, *[E x B) it is more convenient to use their dual 
with respect to the pseudoeucledean *-operator *4{*E), *4{*B), *4{*E x B). 
So, the matrix TZ takes the form 

aEidxAd^ oBidxAd^ a'^{E x B)idx A d$, a^i^E x B)idx A d^ 

aE2dyAd^ aB2dy A d^ a^{E x B)2dy A d^ a'^{E x B)2dy A d^ 

aEsdzAd^ aB^dzAd^ a^{E x B)^dz Ad^ a'^{E x B)sdz A d^ 



Respectively, we obtain 

tr o R = aEidx A d^ + aB2dy Ad^ + a^{E x B)^dz A d^, 



1 
2 

2 „ 2 



(tr o Ry = -a 



{Elf + {B2f + a\E X B) 



In these terms the definitions for the amplitude A and the phase </? will look 
as follows 
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A = -j^R^uapR^"""^, V = arccos 



\{troRy 



:i.36) 



For the solution plane wave the 2-form tr o R is 
troR 



u u 
—===dx AdC+ . . , ., dy AdC + dzAdC 

and in the general case it is not closed. Since costp — uj + p^, then the 
equation d(tr o i?) = requires (in this coordinate system) the following 
conditions to be fulfilled: (^a; = V'y, '{'z = 0. 



With these remarks on the amplitude and phase of the £'M-field in CED 
we come to a close of our short review of the basic principles, concepts and 
relations in Maxwell's theory. As it is clear from what we presented so far 
our purpose is not to describe positive prescriptions for getting results in 
this theory. We have tried to pay attention to those moments of the theory, 
which show directly or indirectly, its frame of applicability. Doing this, we 
get a clearer notion of how and what to alter in view of what kind of objects 
we want to describe. In the next section we summarize those points of of 
Maxwell's theory in order to have clearer and more definite motivation for 
the appropriate to our aims changes in the theory. 



1.4 Why and What to Change in 
Classical Electrodynamics 

If the question why do we want to change something in CED is raised, we 
respond shortly in the following way: because we want to enrich CED with 
new areas of applicability, extending in a natural way the class of admissible 
solutions, aiming to describe (3+1)- dimensional soliton-like objects in the 
pure field case, as well as in the presence of active external fields (media). 

At the end of the last and the beginning of this century it has become 
clear that some experimentally established facts can not be understood and 
explained in the frame of Faraday-Maxwell's electrodynamics. For example. 
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1. Why the initiation of photochemical reactions depends on the color 
and not on the intensity of light? 

2. Why the velocities of the photoelectrons does not depend on the in- 
tensity of light? 

3. Why the shortwave radiation is emitted from bodies at high tempera- 
ture? 

4. Why the shortwave radiation is chemically more active than the long- 
wave one? 

More generally: why the influence of light on matter depends qualitatively 
on its color and not on its intensity! 

These and other experimentally established facts motivate a serious anal- 
ysis of Maxwell's electrodynamics. In result, Planck and, later, Einstein set 
fort the discrete point of view on the structure of the electromagnetic field. 
The later experiments of Compton proved the truth of this viewpoint and 
the notion of photon as an elementary electromagnetic formation (natural 
object) was created. Soon the photons were provided with integral char- 
acteristics like frequency, energy, mom,entum, spin. The Planck constant h 
turned into an omnipresent parameter, serving to separate the real photons 
from other theoretically admissible electromagnetic formations. The Planck's 
formula E — hu is the essential criterion for reality of the objects considered. 
This formula clearly says that only those elementary electromagnetic for- 
mations can really exist, the existence of which is strongly bound up with 
the availability of an intrinsic periodic process with a period of T = l/i/, 
and the corresponding to this process integral action E.T is exactly equal 
to the Planck constant: h — E.T. May be this limitation seems to be very 
strong, but we have no reasons not to trust it so far. Anyway, it is clear 
that the Planck formula separates a class of natural objects being character- 
ized by the elementary action of h. Moreover, this intrinsic periodic process 
and the fact, that every photon moves as a whole uniformly by the same 
velocity c, no matter what its frequency is, support the notion that they are 
finite soliton-like objects. Otherwise it is hardly understandable where the 
quantity frequency can come from as a characteristic of a free and uniformly 
moving point-like, i.e. structureless, object. 

These and other circumstances set a clear chalenge before the those days 
theoretical physics: description of 3-dimensional, finite soliton-like objects, 
having all integral properties of the free photons. The established through 
time quantum-probabilistic approach does not solve this problem since it is 
built on other principles and purcues other objectives. Modern quantum 
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electrodynamics, although its serious achievements in describing the atomic 
phenomena, also works with the assumption "structureless photon" and is 
interested mainly in its integral propeties. 

In view of our intention to build a soliton-like model of the free photon a 
serious analysis of the initial and basic principles of classical electrodynamics 
had to be made. The purpose of this ananlysis is to find what to preserve 
and what to change, i.e. to find those points of the Faraday-Maxwell theory, 
the appropriate change of which will be of use and will not bring us to any 
undesirable complications. Being fully aware of the significance of Faraday- 
Maxwell theory in physics and in all our knowledge of the natural world, we 
present our conclusions fairly and in a transperant way as far as we are able 
to do this. As we said earlier, in doing this we follow the rule that the respect 
and esteem paid to the creators can not be honest and genuine if they are not 
in correspondence with the respect and esteem paid to the Truth. 

1. The elementary spherically symmetric and topologically nontrivial solu- 
tion (see 1.2.2) for the electric field E {B = 0), which could be identified 
with the only spherically symmetric representative of the corresponding co- 
homology class, defines in fact the electric charge as a topological invariant. 
The important point here, we'd like to mention, is the static character of 
the solution-field obtained, so where a test particle, placed somewhere in the 
field, should take momentum from in order to change its own momentum 
according to the momentum conservation law, i.e. the equations of motion, 
is not quite clear. Since -6 = 0, the Poynting vector S — E x B, traditionally 
considered as describing the local momentum-transport of the field, is equal 
to zero. Clearly, the static character of the field is an illusive feature and 
does not give an adequate picture of a real situation. But it is an exact, so 
a trustworthy, solution! Because of the radial direction of the field strength, 
i.e. of the particles' momentum change, it seems hardly possible to avoid 
the notion, that some real objects move radially and carry momentum to and 
from the source. Since there is no momentum accumulation at the source 
object the same momentum has to be carried "to" and "from" for a unit 
time. As for the intrinsic mechanism of momentum exchange between the 
field and the particle CED tells nothing, it gives just the final integral effect. 
For the general static case, according to (1.2.3), because of the well known 
properties of the solutions to Laplace equation, they are not suitable for 
models of real objects. So, our general conclusion is that no static solution 
of Maxwell's equations presents a sufficiently adequate picture of real objects 
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and processes 



2. The computation of the full interaction energy for two spherically sym- 
metric fields in (1.2.3) clearly shows, that while the interaction of the two 
fields is a local fact and takes place at every point where the two fields are 
different from zero, the interaction energy is an integral characteristic and 
in no case should be identified with some potential. From this point of view 
the Coulomb force, describing this interaction as a change of the mechanical 
momentum of the charge- carriers, can not be a local characteristic, since it 
takes into account the contribution to this process of the interaction at all 
points. The "close" form of *uj and dW, where W is considered as a func- 
tion of the points, where the two charges stay, is not a sufficient ground the 
m^e^ra/ characteristic "interaction force" to be identified with *u!. As for the 
potential U, introduced by the relation dU — it is a local characteristic 
too, so dU 7^ dW always. 

3. The considerations in (1.2.4) are of basic significance for us: every localized 
finite initial condition determines unique solution of the wave equation, the 
future time-behaviour of which could shortly be characterized as a "radial 
blow-up". So, the same is true for the pure field Maxwell equations. An 
important feature of the 3-dimensional case is the availability of "forefront" 
and "backfront", which simply means, that every point of the space will 
"feel" the field a finite period of time, after which it will "forget" what 
happened. This result leaves no chance and hopes for making use of Maxwell 
equations to obtain soliton-like solutions, they have no such solutions in the 
whole space. Although undesirable, this conclusion is unavoidable. This is 
the mathematical reality and we have to accept it with the corresponding 
respect. 

4. As we noted in (1.2.5) the popular and well liked solution of the pure field 
Maxwell equations plane wave is unphysical, unrealistic, since it is infinite: it 
possesses infinite integral energy and occupies an infinite 3-volume. Besides, 
all periodic-wave sohitions arc defined by appropriate initial conditions, so 
the wave characteristics like period, frequency and wave vector come from 
these initial conditions, i.e. they are admissible by the equations but are 
not necessary for all solutions. So, other kind of solutions, having no such 
characteristics, are also admissible. But, electromagnetic radiation without 
such characteristics, has never been observed and is hardly possible. If this is 
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true, what to do with such solutions in view of the desired adequacy between 
the theory and the reahty. Note that the representation of the plane wave 
u{z — ct) through simple monochromatic plane waves (the so called wave 
packets) does not solve the problem since these packets are not time-stable 
objects. 

5. The 4-potcntial approach is not of interest from our point of view, since 
for the pure field case it reduces the problem to a subclass of solutions of the 
wave equation for the components of this 4-potential, namely those, sat- 
isfying the additional condition V fxA'^ — 0, and so sohton-hke solutions are 
excluded. In the presence of a current the problems with the nonlinear depen- 
dence of the current on the field, mentioned in (1.2.6), appear and are hardly 
avoidable in general. From our viewpoint the 4-potential approach is impor- 
tant rather to legalize the gauge fields in theoretical fields, although at an 
elementary (linear) abelean level U{1). Unfortunately, even for nonabelean 
gauge fields, leading to the nonlinear Yang-Mills equations, (3-|-l)-soliton- 
like solutions are not found, moreover, there are some studies, showing that 
in the pure field case there are no such solutions. Let's not forget, that the 
so called instantons, i.e. solutions to the equations *F — F, are possible 
only at positively defined space-time metric and zero energy-momentum ten- 
sor, so they could be hardly considered as models of real objects. As for the 
monopoles, they want additional field, interacting with the gauge field in a 
special way. 

6. In the case of continuous media CED adds to the free current an addi- 
tional current j^, called bound. The two new vectors P and Ai arc introduced 
(see 1.1.3) in the same way as the vectors E and B are connected to the free 
current. So, the number of the unknown functions becomes greater than 
the number of the equations, which makes the things difficult to overcome 
as the historic development shows. The introduction of coefficients-material 
constants by means of series developments seems to be a very useful practi- 
cal skill, but it can not serve as a promising theoretical idea. The inherited 
inertia in thinking that energy-momentum exchange with the field can oc- 
cur only if charge carrying particles are present still bears upon our minds. 
The hypothesis that the Faraday induction law is universally (i.e. for all me- 
dia, including vacuum) valid has turned almost into a dogma, which, by the 
way, forbids energy- momentum exchange through *F. Even if some brave 
investigators admit an energy-momentum exchange between the field and the 
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medium to occur through *F, they begin to talk about magnetic charges (in 
analogy with the electric case) having similar to the electric charges proper- 
ties. 

From our point of view the real and important moment is the very energy- 
momentum exchange, and how this exchange is realized is an additional prob- 
lem, depending on the special case considered. So, the most important step 
in our approach is to find an appropriate model equation for this exchange, 
because this is the universal characteristic property of every interaction in 
nature. All further specializations depend on the case under consideration. 

Hence, in view of what was said so far and in view of the purposes we set, 
we come to the following conclusion. The algebraic construction a couple of 
vector fields {E,B) on TZ^, or a differential 2-form on the Minkowski space- 
time, is in general adequate to the field as far as it reflects well enough its 
algebraic and general polarization properties. Maxwell's equations do not 
reflect adequately enough the local properties of the field, therefore not all 
solutions can represent satisfactory models of real objects. So, our choice 
is to preserve, though in an altered form, the basic algebraic picture of the 
field in relativistic terms, but we'll replace Maxwell's equations with new 
nonlinear equations, the physical sense of which is to define how the local 
energy-momentum exchange between the field and some other continuous 
physical object (medium or field) is carried out. The reason to turn to the 
local energy- momentum conservation laws reflects our point of view that they 
are more hopeful and more universal. 

Let us outline in few words our notion of the objects we want to describe. 
As we mentioned in (1.2.1) they must be extended, but finite and time-stable. 
Besides, their existence must be strongly connected with an internal and 
intrinsic dynamics, in particular - periodic process. The characteristics of the 
internal dynamics must be in strong relation with the integral characteristics 
of the object. It is necessary to find invariant quantities, separating the 
really existing objects from all theoretically admissible. If some interaction 
processes take place, a transformation of these objects to each other or to 
new ones, obeying definite conservation laws, should be possible. Various 
internal structures at different levels, as well as creation of stable structures 
out of these objects, should be also possible. A stability with respect to 
external disturbances must be available, so that such a disturbance to result 
finaly in the behaviour of the object as a whole: uniform motion when there 
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are no external disturbances and accelerated motion in presence of a bearable 
external disturbance. From mathematical point of view this means that at 
any moment the functions, describing our object(s), shall be different from 
zero inside a finite connected 3-dimensional region with trivial or non-trivial 
topology, while the time-evolution should be determined by the dynamic field 
equations. 

Passing to the formulation of the basic principles and their mathematical 
adequacies of what will be called further Extended Electrodynamics (EED), 
we are fully aware of, that at this initial stage of our investigation the fol- 
lowing two things should be obeyed. First, from purely pragmatic point of 
view, it seems better to preserve as admissible all solutions to Maxwell's 
equations as exact solutions to the new equations and to use them when- 
ever it is possible. Second, the local energy-momentum conservation laws 
will probably define "weak" equations, so additional conditions (initial data 
or some new equations), reflecting some new speciflc features of the objects 
under consideration, have to be imposed on the solutions. In such cases we 
shall make use of the considerations and conclusions in (1.2.1). We prefer to 
work in rclativistic terms, since we consider this language more adequate to 
the physical situations we want to describe. 

Finally, we assume the general covariance principle, i.e. the understand- 
ing that physical sense may have those concepts and statements, which do not 
depend on the local coordinates used. Accordingly, we'll aim at a coordinate- 
free formulation of the basic statements and equations in all cases when this 
is possible, paying no attention to the simple Minkowski space background 
used. 

1.5 Extended Electrodynamics 

1.5.1 Physical conception for the EM-field in EED 

As it was mentioned, the mathematical models in CED of the real elec- 
tromagnetic fields in vacuum arc "infinite", or if they arc finite, they are 
strongly time-unstable. These models are not consistent with a number of 
experimental facts. A deeper analysis resulted in the new conception for a 
discrete character of the field. This physical understanding of the field is the 
true foundation of EED and it shows clearly the principle differences between 
CED and EED. For the sake of clarity we shall formulate our point of view 
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more explicitly. 

The electromagnetic field in vacuum is of discrete character and consists 
of single, not-interacting (or very weakly interacting) finite objects, called 
photons. All photons move uniformly as a whole by the same velocity c, carry 
finite energy E, momentum p and intrinsic angular momentum. These fea- 
tures imply structure and internal periodic process of period T , which may 
he different for the various photons. The quantity E.T , called "elementary 
action", has the same value for all photons and is numerically equal to the 
Planck constant h. The invariance of c and h means nondistinguishability of 
the photons, considered as invariant objects. The integral value of the intrin- 
sic angular momentum is equal to ±h. For the topology of the 3-dimensional 
region, occupied by the photon at any moment, there are no experimental 
data, so it is desirable the model-solutions to admit arbitrary initial data. 

We'd like to stress once more: EED considers photons as real objects, 
and not as convenient theoretical concepts, and it aims to build adequate 
mathematical models of these entities. So, the first important problem is to 
point out the algebraic character of the modeling mathematical object for a 
single photon. The corresponding generalization for a number of photons is 
easily done (subsec.2.4). 

1.5.2 Choice of the modeling mathematical object 

According to the non-relativistic formulation of CED the electromagnetic 
field has two aspects: "electric" and "magnetic". These two aspects of the 
field are described by two 1-forms on TZ^ and a parametric dependence on 
time is admitted: the electric field E and the magnetic field B. The consid- 
erations made in (1.1.1) brought us to the conclusion that these two fields 
can be considered as two vector com,ponents of a new object, 1-form Q, tak- 
ing values in a real 2- dimensional vector space, naturally identified with TZ^ . 
This mathematical object unifies and, at the same time, distinguishes the two 
sides of the field: there is a basis in TZ'^, in which the electric and magnetic 
components are delimited, but in an arbitrary basis the two components 
mix (superimpose), so the difference between them is deleted. The physi- 
cally important quantity, energy density, is given by the sum E^ + . This 
point of view seems appropriate and relevant in view of the pointed out in 
(1.1.1) invariance of Maxwell's equations with respect to some linear trans- 
formations, mixing E and B. All unimodular such transformations keep the 
energy-density unchanged. 
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In the relativistic formulation of CED the difference between the electric 
and magnetic components of the field is already quite conditional, and from 
invariant-theoretical point of view there is no any difference. However, the 
2-component character of the field is kept in a new sense and manifests itself 
at a different level. In fact, as we mentioned above, some linear combina- 
tions of the electric and magnetic fields generate a new solution to Maxwell's 
equations. In particular, such is the transformation, defined by the matrix 



1 
-1 



This matrix, defining a complex structure in TZ^, transforms a field of the 
kind {E, 0) into a new field of the kind (0, E) and a field of the kind (0, B) 
into a field of the kind {—B, 0), i.e. the electric component into magnetic and 
vice versa. This observation draws our attention to looking for a complex 
structure J, = —id in the bundle of 2-forms on the Minkowski space, such 
that if F presents the first component of the field, then J(-F) to present the 
second component of the same field. Such complex structure truly exists and, 
according to (1.2.1), it coincides with the restriction of the Hodge *-operator, 
defined by the pseudometric to the space of 2-forms: **2 = —id. So, 
the non-relativistic vector components B) are replaced by the relativistic 
vector components {F,*F). The following considerations support also such 
a choice. 

The relativistic Maxwell's equations in vacuum dF = 0, d * F = are, 
obviously, invariant with respect to the interchange F *F. Moreover,if F 
is a solution, then an arbitrary linear combination aF + b * F is again a so- 
lution. More generally, if {F, *F) defines a solution, then the transformation 
(F, *F) — > {aF + b* F, mF + n* F) defines a new solution for an arbitrary 
matrix 

" a m 
b n 

Now, using the old notation fl for the new object fl = F ^ ei + *F ® 62, 
Maxwell's equations are written down as dJl = 0, or equivalently SH, — 0. 
Clearly, an arbitrary linear transformation of the basis (ei, 62) keeps Q as a 
solution. 

Recall the energy- momentum tensor in CED, defined by (1.7) 



1 



1 

Stt 



-F,,F^''-{*F),,{*Fy 
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It is quite clearly seen, that F and *F participate in the same way in Q^, 
and the full energy-momentum densities of the field are obtained through 
summing up the energy-momentum densities, carried by F and *F. Since 
the two expressions F^^F"'^ and {*F)^a-{*Fy^ are not always equal, the 
distribution of energy-momentum between F and *F may change in time, 
i.e. energy- momentum may be transferred from F to *F, and vice versa. So 
we may interpret this phenomenon as a special kind of interaction between F 
and *F, responsible for some internal redistribution of the field energy. Now, 
in vacuum it seems naturally to expect, that the energy-momentum, carried 
from F to *F in a given 4- volume, is the same as that, carried from *F to F in 
the same volume. However, in presence of an active external field (mcdlTim), 
exchanging energy-momentum with Q, it is hardly reasonable to trust the 
same expectation just because of the specific structure the external field 
(medium) may have. So, the external field (further any such external field 
will be called medium for short) may exchange energy-momentum preferably 
by F or *F, as well as it may support the internal redistribution of the 
field energy-momentum between F and *F, favouring F or *F. From the 
explicit form of the energy-momentum tensor it is seen that the field may 
participate in this exchange by means of any of the two terms F^F'^"' and 
{*F)ij,a{*Fy^. Moreover, from the expression (1.8) for the divergence of the 
energy-momentum tensor 



F^{SFr + (tF)^(S * FY 



it is also clearly seen that the quantities of energy-momentum, which any 
of the two components F and *F may exchange in a unit 4- volume are 
invariantly separated and given respectively by 

F^^idFY, {*F)^,{d*Fy. 

But, in CED the exchange through *F is forbidden, the expression 
{*F)^i,{6 * FY is always equal to zero. Of course, we do not reject the exis- 
tence of such media, but we do not share the opinion that all media behave 
in this same way just because this can not be verified. On the other hand, 
in case of vacuum, we can not delimit F from *F, these are two solutions 
of the same equation and it is all the same which one will be denoted by F 
(or *F), i.e. CED does not give an intrinsic criterion for a respective choice. 
Only in regions with non-zero free charges and currents, when dF = and 
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5F = J 7^ 0, the choice can be made , but this presupposes (postulates) that 
the field is able to interact, i.e. to exchange energy- momentum, only with 
charged particles. This postulate we can not assume ad hoc. 

Having in view these considerations we assume the following postulate in 
EED in order to specify the algebraic character of the modehng mathematical 
object: 

In EED the electromagnetic field is described by a 2-form fl, defined on 
the Minkowski space-time and valued in a real 2- dimensional vector space W 
and such, that there is a basis (ei, 62) ofW in which takes the form 

Vt = F ®ei + *F ®e2. (1.37) 

Since W is isomorphic to Ti?, further we shall write only Ti? and all 
relations obtained can be carried over to W by means of the corresponding 
isomorphism. In particular, every W will be considered as being provided 
with a complex structure J, so, the group of automorphisms of J is defined. 
Our purpose now is to prove that the set of 2-forms of the kind (1-37) is 
stable under the invariance group of J. 

First we note, that the equation aF + 6 * F = requires a = 6 = 0. In 
fact, if a 7^ then F = -J*-F. From aF+h*F = we get a*F-bF = and 
substituting F, we obtain (a^ + 6^) *F = 0, which is possible only if a = 6 = 
since *F 7^ 0. In other words, F and *F are linearly independent. Let now 
{ki, ^2) be another basis of 7^^ and consider the 2-form ^ = G(E)ki + *G®k2- 
We express {ki, ^2) through (ci, 62) and take in view what we want 

G ^ ki + *G ^ k2 ^ G ^ (aei + 662) + *G ^ (mei + 7162) = 

= {aG+m*G)(^ei + {bG+n*G)(S>e2 = {aG+m*G)(S>ei + *{aG+m*G)(^e2. 

Consequently, bG + n*G — a*G — mG, i.e. (6 -|- m)G -\- {n — a) * G — 
which requires m — —b, n — a, i.e. the transformation matrix is 

a —b 
b a 

Besides, if Qi and Q2 are of the kind (1.37), it is easily shown that the hnear 
combination XQi+fiQ2 is of the same kind (1.37). These results show that the 
2-forms of the kind (1.37) form a stable with respect to the automorphisms 
of {W^, J) subspace of the space A^{M, U^). 
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Wc note that the following 2-forms: F®ei + *F(g)e2, F^ki + *Fi^k2, are 
different, i.e. it is important which basis will be used. In order to separate 
the class of bases we are going to use, first we recall the product of 2 vector 
valued differential forms. If $ and ^ are respectively p and q forms on 
the same manifold N, taking values in the vector spaces Wi and W2 with 
corresponding bases (ei,...,e^) and (/ci, /c„), and ip : Wi x W3 is 

a bilinear map into the vector space W^, then a (p + g)-form ip ^) on 
with values in W3 is defined by 

if = ^ A (g) (p{ei, kj). 

In particular, if Wi — W2 and W3 — TZ, and the bilinear map is scalar (inner) 
product we get 

($,*) = ^ A 
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Let now X and F be 2 arbitrary vector fields on the Minkowski space M, 
be of the kind (1.37), Q/^i/ be the energy tensor in CED and g be the 
canonical inner product in Ti?. Then the class of bases we shall use will be 
separated by the following equation 

Q^^X^'Y'' = ]^*g(i{X)n,*i{Y)VL). (1.38) 

We develop the right hand side of this equation and obtain 

]^*g(i{X)n,*i{Y)Q) = 

^ * g{i{X)F (g) ei + i{X) * F (g) 62, H{Y)F (g) ei + H{Y) * F (g) 62) = 
[i{X)F A *i(y)F)^(ei, ei) + [i{X)F A *i(y) * F)g{e^, 62)+ 
+ (i(X) * F A *i(r)F)^(e2, ei) + (i(X) * F A *i{Y) * F)c/(e2, 62 



1 

— * 
2 



— X'^y^ 
2 



+ (f^,(*F)^ + (*F)^,F;)5(ei, 62) 



F^,F^g{e,,e,) + (*F)^,(*F)^^(e2, 62)+ 
1 



2 



F^,F; + (*F)^.(*F)^ 
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In order this relation to hold it is necessary to have 

g{ei,ei) = 1, 9(62,62) = 1, 9(61,62) = 0, 

i.e., we are going to use orthonormal bases. So, the stability group of the 
subspace of forms of the kind (1.37) is reduced to S0(2) or U(l). So, in 
this approach, the group 80(2) appears in a pure algebraic way, while in 
the gauge interpretation of CED this group is associated with the equation 
dF = 0, i.e. with the traditional and not shared bu us understanding, 
that the £'M-field can not exchange energy-momentum with some medium 
through *F. 

1.5.3 Differential equations for the field 

We proceed to the main purpose, namely, to write down differential equations 
for our object f2, which was chosen to model the i?M-field. We shall work 
in the orthonormal basis (61,62), where the field has the form (1.37). The 
two vectors of this basis define two mutually orthogonal subspaces {ci} and 
{62}, such that the space TZ^ is a direct sum of these two subspaces: TZ^ — 
{ci} ® {62}. So, wc have the two projection operators tti : Ti? — > {ei}, 7r2 : 
'R? — > {62}. These two projection operators extend to projections in the 
72.^-valued differential forms on M: 

TTiQ = 1:1(9} (8) /ci + (8) k2) ® TTiki + (8) 7riA;2 = 

= 9^ ® 7ri(a6i + 662) + 9'^ <^ 7ri(m6i + 7162) = (a9^ + m9'^) ® 61. 
Similarly, 

-K29 = (b9i + n92) ® 62. 
In particular, if 9 is of the form (1.37), then 

7ri(F (2) 61 + *F (g) 62) = F (g) 61, n2(F (g) 61 + *F (g) 62) = *F (g) 62. 

Let now our FM-field 9 propagates in a region, where some other contin- 
uous physical object (external field, medium) also propagates and exchanges 
energy-momentum with 9. We are going to define explicitly the local law 
this exchange obeys. 

First we note, that the external field is described by some mathematical 
object(s). From this mathematical object, following definite rules, refiecting 
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the specific situation under consideration, a new mathematical object Ai is 
constructed and this new mathematical object participates directly in the 
exchange defining expression. The ii^M-field participates in this exchange 
defining expression directly through fl, and since takes values in TZ^, then 
Ai must also take values in TZ"^. 

We make now two preliminary remarks. First, all operators, acting on 
the usual differential forms, are naturally extended to act on vector valued 
differential forms according to the rule D ^ D x id. In particular. 



Second, in view of the importance of the expression (1.8) for the divergence 
of the CED energy-momentum tensor, we give its explicit deduction. Recall 
the following algebraic relations on the Minkowski space: 

Sp — (— l)^*~^d* — *d*, 5*p — *dp forp — 2k+l, S*p — —*dp forp = 2k. 



(1.39) 



If q; is a 1-form and F is a 2-form, the following relations hold: 



* (a A *F) 




In particular, 




Having in view these relations, we obtain consecutively: 



- F »F"^5'' - F F""^ = 



^F"^[(dF),^^ - V«F^^ - VpF^] - {y,F^,)F^'^ - F^V^F 
]^F^P{dF)^p^ - F^^V^F-^ = -{*F)^,V,{*Fy - F^^V ,F^ 
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^{*F)^,{S*Fy + F,,{SFY. 

Let now our field Q interact with some other field. This interaction, i.e. 

energy- momentum exchange, is performed along 3 "channels". The first 2 
channels are defined by the 2 (equal in rights) components F and *F of 
fl. This exchange is real in the sense, that some part of the £^M-energy- 
momentum may be transformed into some other kind of energy-momentum 
and assimilated by the external field or dissipated. Since the two components 
F and *F are equal in rights it is naturally to expect that the corresponding 
2 terms, defining the exchange in a unit 4-volume, will depend on F and *F 
similarly. The above expression for Vj^Q)! gives the two 1-forms 

F^^iSFydx'^, {*F)^,{5*Fydx^ 

as natural candidates for this purpose. As for the third channel, it takes into 
account a possible infiuence of the external field on the intra-field exchange 
between F and *F, which occurs without assimilation of energy-momentum 
by the external field. The natural candidate, describing this exchange is, 
obviously, the expression 

F^,{5 * Fydx>' + {*F)^,{8Fydx''. 

It is important to note, that these three channels are independent in the sense, 
that any of them may occur without taking care if the other two work or don't 
work. A natural model for such a situation is a 3-dimensional vector space 
K, where the three dimensions correspond to the three exchange channels. 
The non-linear exchange law requires some valued non-linear map. Since 
our fields take values in Ti? this 3-dimensional space must be constructed 
from 'R? in a natural way. Having in view the bilinear character of V^/QJ^ 
it seems naturally to look for some bilinear construction with the properties 
desired. These remarks suggest to choose the symmetrized tensor product 
Symin"^ (g) 7^2) = 7^2 V TZ'^. So, from the point of view of the £;M-field, the 
energy-momentum exchange term should be written in the following way: 

y{6n,*n). (i.4i) 

In fact, in the corresponding basis (ei, 62) we obtain 

V(5Q, *n) = y(5F ® ei + 5 * F ® e2,*F ® ei + * * F ® 62) ^ 
= (5FA*F)(g)ei Vei-F(5*FA**F)(8)e2Ve2 + (-5FAF-F5*FA*F)(8)eiVe2. 
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This expression determines how much of the EM-field energy-momentum has 
been carried irreversibly over to the external field and how much has been 
redistributed between F and *F by virtue of the external field infiuence in a 
unit 4-volume. 

Now, this same quantity of energy-momentum has to be expressed by new 

terms, in which the external field "agents" should participate. Let's denote by 
$ the first agent, interacting with ttiQ, and by ^ the second agent, interacting 
with 7T2ft. Since the corresponding two exchanges are independent, we may 
write the exchange term in the following way: 



According to the local energy-momentum conservation law these two quan- 
tities have to be equal, so we obtain 



This is the basic relation in EED. It contains the basic differential equa- 
tions for the EM-iie\d components and requires additional equations, spec- 
ifying the properties of the external field, i.e. the algebraic and differential 
properties of $ and ^f. The physical sense of this equation is quite clear: 
local balance of energy-momentum. Further we shall study this relation in 
various cases, and in the first place - the vacuum. 



V *7riO) + V(*,*7r20). 



(1.42) 



(1.43) 
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Chapter 2 

Extended Electrodynamics in Vacuum 



2.1 Vacuum Equations in EED 
2.1.1 Vacuum in EED 

In CED the term "vacuum" is used in the sense, that in the region, where we 
consider an EM-field, there is no free or bound electric charges. This notion 
of vacuum in CED comes from the conception, that the field may exchange 
energy-momentum only with electric charge carriers and only through the 
component F. EED extends the possibilities for energy-momentum exchange 
assuming that the full quantity of exchanged energy-momentum in a unit 4- 
volume is given by the general expression (1.42). That's why in EED we talk 
about a field Q in vacuum every time when this expression (1.42) is equal 
to zero. FormaUy this means that every external field (medium) , which does 
not exchange energy-momentum with Q, can be treated as vacuum as far as 
the £'M-field is concerned. Talking about exchange, we mostly have in mind 
that some energy-momentum is transferred from the field to the medium, 
however, we do not formally exclude the reverse process. 

Assume now that in some region we have an £'M- field Q. In the corre- 
sponding basis (ei, 62) we can write 

n = F ei + *F IS) 62, $ = ei + (g) 62, ^' = Ci a'' 62- 

Remark. Further the 1-forms a*, i = 1,...4, as well as the corresponding 
through the pseudometric r] vector fields, will be called shortly currents. 
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Of course, this terminology should not be associated with some charged 
particles, or with some before given specific structure. In the general case 
these currents are just the tools of the external field to gain some energy- 
momentum from the field Q. 

The expression (1.42) looks as follows: 

V($, ^TTiil) + V(^, *n2^) = A *F (g) d V ei + a"^ A * * F (g) 62 V 62+ 

+(a^ A * * F + A *F) (8) ei V 62. 
Clearly, the necessary and sufficient condition for vacuum is 

A *F = 0, A * * F = 0, A * * F + a'^ A *F = 0, (2.1) 

or, in components 

= 0, {*F)^,al = 0, (*F)^,a^ + F^,a^ = 0. (2.2) 

We see that there are various possibilities, i.e. relations among F and a^, to 
realize a vacuum situation. The strongest condition is, of course, $ = \I' = 0, 
i.e. all currents are equal to zero: a* = 0, i = 1,...,4. If, at least one of 
the two currents and is different from zero, then the above relations 
(1.45) require ci6t||F^,^|| = 0, i.e. F A F = 0, or orthogonality between E 
and B. If a*^ = «^ = 0, but the other two currents are different from zero 
then F A F 7^ in general. In accordance with our interpretation of the 
equations (1.42)this means that the medium affects the exchange between F 
and *F without gaining and keeping any energy-momentum. In such a case 
the orthogonality between E and B is not needed. It seems important to 
note, that (2.1) requires every couple o;*, to be in the kernel of F and *F, 
i.e. 

F{a\ a^) = i*F){a\ a^) = 0. (2.3) 

So, in case of vacuum, the currents are strongly dependent on the field Q. 
The above relations define equations, connecting the 16 components of the 
4 currents with the components F^j,. Of course, in curvelinear coordinates 
{y'^} these equations will depend strongly on the metric coefficients r]^y{y'^) 
in these coordinates. Finally we note, that and may be considered as 
eigen vectors respectively for F and *F at zero eigen values, which, according 
to the formulas in (1.1.2), is always possible if I2 = 2E.B = 0. 
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2.1.2 Equivalent forms of the equations 

According to the last subsection an external field is called vacuum with re- 
spect to the EM-fie\d fl if the right hand side of (1-43) is equal to zero. Then 
the left hand side of (1-43) will also be equal to zero, so we get the equations 

V {Sn, *Q) = 0. (2.4) 

From this coordinate free compactly written expression we obtain the follow- 
ing equations for the components of in the basis (ei, 62): 



{6F A *F) (g) ei V ei + (5 * F) A * * F) (g) 62 V 62 + (2.5) 
+ {SF A * * F 5 * F A *F) O 61 V 62 = 0. 

So, the field equations, expressed through the operator S, look as follows 

(5FA*F = 0, 5*FA**F = 0, S * F A *F - SF A F ^ 0. (2.6) 

These equations, expressed through the operator d, have the form 

*FA*d*F = 0, FA*dF = 0, F A *d * F + *F A *dF = 0. (2.7) 

Using the components F^j^, we obtain for the equations (2.6) 

F,,{6Fy = 0, i*F),,{6 * FY = 0, F^,{6 * Ff + {^F),,{5Fy = 0. (2.8) 

In the same way, for the equations (2.7) we get 

(*F)^'^(d*F)^., = 0, F'^''(dF)^., = 0, (*F)^'^(dF)^., + F'''^(d*F)^., = 0. 

(2.9) 

Now we give the 3-dimensional form of the equations in the same order: 
rotB - —j - EdivE = 0, E. irotB - —j = 0, (2.10) 



(dB \ I dB \ 

rotE + — j - BdtvB = 0, B. irotE + —j = 0, (2.11) 



^ dB\ ( dE\ 

rotE + j X 5 + WotB - — j X E + BdivE + EdivB = 0, 
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B. [rotB -^)-E. {rotE + ^) = 0. (2.12) 

From the second equations of (2.10) and (2.11) the well known Poynting 
relation follows 

div {ExB) + — = 0, 

and from the second equation of (2.12), if E.B = g{x,y,z), we obtain the 
known from Maxwell theory relation 



B.rotB = E.rotE. 



The explicit form of equations (2.10) end (2.11) should not make us conclude, 
that the second (scalar) equations follow from the first (vector) equations. 
Here is an example: let divE = 0, divB = and the time derivatives of E 
and B are zero. Then the system of equations reduces to 

E X rotE = 0, S X rotB = 0, B.rotE = 0, E.rotB = 0. 

As it is seen, the vector equations do not require any connection between 
E and B in this case, therefore, the scalar equations, which impose such 
a connection, can not follow from the vector ones. The third equations of 
(2.7) and (2.8) determine ( in equivalent way) the energy-momentum quan- 
tities, transferred from F to *F, and reversely, in a unit 4-volume, with the 
expressions, respectively 

F,,(5 * FY = -{*F),,{5Fr, F^'^id * F),,, = -(*F)'''^(dF)^.,. 

From these relations it is seen, that the 1-forms 6F and d*F play the role of 
"eoctemal" currents respectively for *F and F. In the same spirit we could 
say, that the energy-momentum quantities F^i,{SFy and {*F)^,^{S * F)^ , 
which F and *F exchange with themselves, are equal to zero. And this 
corresponds fully to our former statements, concerning the physical sense of 
the equations for the components of Vt. 



2.1.3 Conservation laws 

From the first two equations of (2.8) and from the earlier given expression for 
the divergence ViyQJ^ of the Maxwell's energy-momentum tensor Q'i in CED 



64 



it is immediately seen that on the solutions of our equations (2.8) this diver- 
gence is also zero. In view of this we assume the tensor Q'j^, defined by (1.7) 
to be the energy-momentum tensor in EED. We shah be interested in finding 
explicit time-stable solutions of finite type, i.e. F^i, to be finite functions of 
the three spatial coordinates, therefore, if it turns out that such solutions 
really exist, then integral conserved quantities can be easily constructed and 
computed, making use of the 10 Killing vectors on the Minkowski space-time. 
We recall that in CED such finite and time-stable solutions in the whole space 
are not allowed by the Maxwell's equations. 

2.2 General Properties of the Equations and 
Their Solutions 

2.2.1 General properties of the equations 

We first note, that in correspondence with the requirement for general covari- 
ance, equations (2.4), given above and presented in different but equivalent 
forms, are written down in coordinate free manner. This requirement is uni- 
versal, i.e. it concerns all basic equations of a theory and means simply, 
that the existence of real objects and the occurrence of real processes can 
not depend on the local coordinates used in the theory., i.e. on the conve- 
nient for us way to describe the local character of the evolution and structure 
of the natural objects and processes. Of course, in the various coordinate 
systems the equations and their solutions will look differently. Namely the 
covariant character of the equations allows to choose the most appropriate 
coordinates, reflecting most fully the features of every particular case. A 
typical example for this is the usage of spherical coordinates in describing 
spherically symmetric flelds. Let's not forget also, that the coordinate-free 
form of the equations permits an easy transfer of the same physical situation 
onto manifolds with more complicated structure and nonconstant metric ten- 
sor. Shortly speaking, the coordinate free form of the equations in theoretical 
physics reflects the most essential properties of reality, called shortly objective 
character of the real phenomena. 

Since the left hand sides of the equations are linear combinations of the 
first derivatives of the unknown functions with coefficients, depending linearly 
on these unknown functions, (2.8) present a special type system of quasilin- 
ear first order partial differential equations. The number of the unknown 
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functions F^^^ = —F^^ is 6, and in general, the number of the equations is 
3.4 = 12, but the number of the independent equations depends strongly on 
if the two invariants Ii = ^F^^^F^^ and I2 = \{*F)ni,F^^" are equal to zero or 
not equal to zero. If /2 7^ then det{F^^) ^ and the first two equations of 
(2.8) are equivalent to = and (5 * F = 0, which automatically eliminates 
the third equation of (2.8), i.e. in this case our equations reduce to Maxwell's 
equations. 

It is clearly seen from the (2.9) form of the equations, that the metric ten- 
sor essentially participates (through the ^-operator applied to 2-forms only) 
in the equations. If we use the (5-operator, then the metric participates also 
through the *-operator, applied to 3-forms, but this does not lead to more 
complicated coordinate form of the equations. It is worth to note that in 
nonlinear coordinates the metric tensor will participate with its derivatives, 
therefore, the very solutions will depend strongly on the metric tensor cho- 
sen. This may cause existence or non-existence of solutions of a given class, 
e.g. soliton-like ones. In our framework such additional complications do 
not appear because of the opportunity to work in global coordinates with 
constant metric tensor. 

We note 2 important invariance (symmetry) properties of our equations. 

Property 1. The transformation F *F does not change the system. 

The proof is obvious, in fact, the first two equations interchange, and the 
third one is kept the same. In terms of fl this means that if is a solution, 
then *Q, is also a solution, which means, in turn, that equations (2.4) are 
equivalent to equations 

V(Q,*dQ) = 0. (2.13) 

Property 2. Under conformal change of the metric the equations do not 
change. 

The proof of this property is also obvious and is reduced to the notice, that 
as it is seen from (2.9), the *-operator participates only with its reduction to 
2-forms, and as we noted in subsec.(1.1.2), *2 is conformally invariant. 

Summing up the first two equations of (2.10) and (2.11) we obtain how 
the classical Poynting vector changes in time in our more general approach: 

d 

— (ExB)^ EdivE + BdivB -Ex rotE -Bx rotB. 
In CED the first and the second terms on the right are missing. 
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Here is an example of static solutions of (2.4), which are not solutions to 
Maxwell's equations. 

E — {asinaz, acosaz, 0), B — {bcosaz, —bsinaz, 0), 

where a, b and a are constants. We obtain 

rotE = {aasinaz, aacosaz, 0), rotB = {bacosaz, —basinaz, 0) 

Obviously, E x rotE = 0, B x rotB = 0, E.rotB = 0, B.rotE = 0. For 
the Poynting vector we get E x B = (0, 0, —ah), and for the energy density 
w = |(a^ + 6^). Considered in a finite volume, these solutions could model 
some standing waves, but we do not engage ourselves with such interpre- 
tations, since we do not accept seriously that static E'M-fields may really 
exist. 

2.2.2 General properties of the solutions 

It is quite clear that the solutions of our equations are naturally divided into 
two classes: linear and nonlinear. The first class consists of all solutions 
to Maxwell's vacuum equations, where the name linear comes from. These 
solutions are well known and won't be discussed here. The second class, 
called nonlinear, includes all the rest solutions. This second class is natu- 
rally divided into two subclasses. The first subclass consists of all nonhnear 
solutions, satisfying the conditions 

5F^0, 5*F^0, (2.14) 

and the second subclass consists of those nonlinear solutions, satisfying one 
of the two couples of conditions: 

SF=0, S*F^O; SF^O, 6 * F = 0. 

Further we assume the conditions (2.14) fulfilled, i.e. the solutions of the 
second subclass will be considered as particular cases of the first subclass. 
Our purpose is to show explicitly, that among the nonlinear solutions there 
are soliton-like ones, i.e. the components F^i, of which at any moment are 
finite functions of the three spatial variables with connected support. We 
arc going to study their properties and to introduce corresponding charac- 
teristics. First we shall establish some of their basic features, proving three 
propositions. 
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Proposition 1. All nonlinear solutions have zero invariants: 



h = \f^..f^'' = 0, h 



0. 



Proof. Recall the field equations in the form (2.8): 

F^,{5FY = 0, {*F)^,{5 * FY = 0, F^,{5 * Ff + {*F)^,{5Fy = 0. 

It is clearly seen that the first two groups of these equations may be con- 
sidered as two linear homogeneous systems with respect to and 5 * F^^ 
respectively. In view of the nonequalities (2.14) these homogeneous systems 
have non-zero solutions, which is possible only if det^F^y) = det{*F)^y) = 0, 
i.e. if I2 — 2E.B — 0. Further, summing up these three systems of equations, 
we obtain 

(F + *F)^,(5F + 5*F)" = 0. 
If now {8F + 5* FY ^ 0, then 



= det{F + *F) 



hF + *F),y{*F-Fr' 



- [-2F,yFn' = (/i 



If dF'' = —{6* Fy 7^ 0, we sum up the first two systems and obtain 
(*F - F)^y{5 * FY = 0. Consequently, 



= det{*F - F)^, = ^(*F - F)^,{-F - ^FY""^ ' 



This completes the proof. 

Recall that in this case the energy-momentum tensor Q/^i, has just one 
isotropic eigen direction and all other eigen directions are space-like. Since 
all eigen directions of F^j, and are eigen directions of Q^i, too, it is clear 
that Ffj^i, and {*F)^y can not have time-like eigen directions. But the first 
two systems of (2.8) require 5F and 5 * F to be eigen vectors of F and *F 
respectively, so we obtain 



{SF).{SF) < 0, {S* F).{S * F) < 0. 



(2.15) 



Proposition 2. All nonlinear solutions satisfy the conditions 

{SF)i,{6*FY = 0, \SF\ = \S*F\ 



(2.16) 
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Proof. We form the inner product i{S * F){SF A *F) — and get 

{5 * FY{5F)i,{*F) -5FA{5* FY{*F)^^dx'' = 0. 

Because of the obvious nuhfication of the second term the first term will be 
equal to zero (at non-zero *F) only if {SF)n{S * FY = 0. 

Further we form the inner product i{5*F) (5F AF — 5 * F A *F) — and 
obtain 

{5 * FY{SF)^F -5FA{S* F^F^^dx"- 
-(5 * F)2(*F) + 5*FA(5* FY(*F)i,^dx'' = 0. 

Clearly, the first and the last terms are equal to zero. So, the inner product 
hy SF gives 

{SF^iS * FYF^^dx"" - [{5FY{5 * FfF^,] 5F+{5* Ff{8FY{*F)^,dx'' = 0. 

The second term of this equality is zero. Besides, {5 * FYF^i^dx'^ — 
-{5FY{*F)^,dxr So, 

\8Ff - (5 * Ff 

Now, if {5FY{*F)^i,dx'' ^ 0, then the relation \6F\ = \6 * F\ follows imme- 
diately If {6FY{*F)^^dx'' = = -(5 * FYF^ydx'' according to the third 
equation of (2.8), we shall show that {5FY — [5 * FY — 0. In fact, forming 
the inner product i{SF){5F A *F) = , we get 

{5FY *F-5FA {5FY{*F)i,^dx'' = {5FY * F = 0. 

In a similar way, forming the inner product i{5 * F)5 * F A F — we have 

(5 * FYF - S{*F) A{S* FYF^.udx'' = {5* FYF = 0. 

This completes the proof. 

We just note that in this last case the isotropic vectors 5F and 5*F are eigen 
vectors of Q^^ too, and since Q^^ has just one isotropic eigen direction, we 
conclude that 5F and 5 * F are colinear. 

In order to formulate the third proposition, we recall from subsec. (1.1.2) 
that at zero invariants Ji = J2 = the following representation holds: 

where C, is the only (up to a scalar multiple) isotropic eigen vector of Q^. 
Also, the relations A.C, =0, A* .C, — Q are in force. Having this in view we 
shall prove the following 



{5FY{*F)^,dx'' = 0. 
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Proposition 3. All nonlinear solutions satisfy the relations 



C^(5F)^ = 0, C^(5*F)^ = 0. 



(2.17) 



Proof. We form the inner product i{Q{SF A *F) = : 



[C^{5F)^] A* AC- {SF A CK'{A*), + (SF A A*)C^C, = 0- 



Since the second and the third terms are equal to zero and *F ^ 0, then 
C^(5-F)/i = 0. Similarly, from the equation {6 * F) A F — we get 
(^(5 * F)ij, — 0. The proposition is proved. 

2.2.3 Algebraic properties of the nonlinear 
solutions 

Since all nonlinear solutions have zero invariants Ii — I2 — we can make 
a number of algebraic considerations, which clarify considerably the struc- 
ture and make easier the study of the properties of these solutions. As we 
mentioned earlier, all eigen values if F, *F and Q^^, in this case are zero, 
and the eigen vectors can not be time-like. There is only one isotropic direc- 
tion, defined by the isotropic vectors ±( and the representations F — AA(, 
*F = A* A C hold, moreover, we have A.A* = 0, = {A*y < 0,A.C = 
A*.( = 0. Recall that the two 1-forms A and A* are defined up to isotropic 
additive factors, colinear to (. The above representation of F and *F through 
C shows that these factors do not contribute to F and *F, therefore, we as- 
sume further that, these additive factors are equal to zero. 

We express now Q^,^ through A, A* and (. First we normahze the vector (. 
This is possible, because it is an isotropic vector, so its time-like component 
(■4 is always different from zero. We divide by (4 and get the vector 
V = (V^, V^, V^, 1), defining, of course, the same isotropic direction. Now 
we make use of the identity (1.25), where we put F^^, instead of G^,^. Having 
in view that h = ^F^^F^"" = 0, we obtain F^F''" = {*F)^^{*Fy^. So, the 
energy-momentum tensor looks as follows 



O'" = —F F' 





V 



1 



(2.18) 
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This choice of C = V determines the following energy density 

AtiQI = \A\' = \A*\\ 

We consider now the influence of the conservation law V^Q^^ = on V. 

V,Ql = -A VV.V^ - V^V. {A^Y"") = 0. 

This relation holds for every /i = 1, 2, 3, 4. We consider it for // = 4 and get 
V'V^l) = V'^a^(l) = 0. Therefore, V4V^ (A^V^) = {A^W") = 0. Since 
7^ 0, we obtain that V satisfies the equation 

V'V^V'^ = 0, 

which means, that V is a geodesic vector field, i.e. the integral trajectories 
of V are isotropic geodesies, or isotropic straight lines. Hence, every nonlin- 
ear solution F defines unique isotropic geodesic direction in the Minkowski 
space-time. This important consequence allows a special class of coordinate 
systems, called further F-adapted, to be introduced. These coordinate sys- 
tems are defined by the requirement, that the trajectories of the unique V, 
defined by F, to be parallel to the (z, ^)-coordinate plane. In such a coor- 
dinate system we have = (0,0, e, 1), e = ±1. Further on, we shall work 
in such arbitrary chosen but fixed F-adapted coordinate system, defined by 
the corresponding F under consideration. 

We write down now the relations F — A AV, *F — A* AV component- 
wise, take into account the values of in the F-adapted coordinate system 
and obtain the following explicit relations: 

Fi2 = F34 = 0, Fi3 = £Fi4, F23 = BF24, 
(*F)i2 = (*F)34 = 0, (*F)i3 = £(*F)i4 = -F24, (*F)23 = e{*F)24 = Fi4, 

A = (Fi4, F24, 0, 0) , A* = (-F23, Fi3, 0, 0) = {-eA2, sA^, 0, 0) . (2.19) 

Clearly, the 1-forms A and —A* can be interpreted as electric and magnetic 
fields respectively. Only 4 of the components Q'j^ are different from zero, 
namely: Qf = —Ql = eQl = —eQl — Introducing the notations 

Fi4 = u, F24 = p, we can write 

F — eudx Adz-\- udx Ad^-\- epdy Adz -\- pdy A d^ 

*F = —pdx Adz — epdx Ad^-\- udy Adz -\- eudy A d^. 
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In the important for us spatially finite case, i.e. when the functions u and p 
are finite with respect to the spatial variables {x, y, z), for the integral energy 
W and momentum p we obtain 

W — j Q\dxdydz — j {u^ + p^)dxdydz < oo, 

p^(o,0,e^y ^c^\p\^-W^ = Q. (2.20) 

Now wc show how the nonlinear solution F defines at every point a pseu- 
doorthonormal basis in the corresponding tangent and cotangent spaces. The 
nonzero 1-forms A and A* are normed to A = and A* = A*/|y4*|. Two 

new unit 1-forms R and S are introduced through the equations: 

R2 = -1, A^R^ = 0, (A*)^R^ = 0, V^R^ = £, S = V + eR. 

The only solution of the first 4 equations is R^^ — (0, 0, —1, 0). Then for S we 
obtain — (0, 0, 0, 1). Clearly, R^ = —1 and = 1. This pseudoorthonor- 
mal (co-tangent) basis is carried over to a (tangent) pseudoorthonormal basis 
by means of the pseodometric rj. 

We proceed further to introduce the concepts of amplitude and phase 
in a coordinate-free manner. Of course, we shall use the considerations in 
subsec. (1.3.1). First, of course, we look at the invariants, we have: I1 — I2 — 
0. But in our case we have got another invariant, namely, the module of 
the 1-forms A and A*: \A\ = \A*\. Let's begin with the amplitude, which 
shall be denoted by 0. As it's seen from the above obtained expressions, the 
magnitude of |A| coincides with the square root of the energy density in any 
F-adapted coordinate system. As we noted in (1.3.1) this is the sense of the 
quantity amplitude. So, we define it by the module of \A\ — \A*\. We give 
now two more coordinate-free ways to define the amplitude. 

Recall first, that at every point, where the field is different from zero, we 
have three bases: the pseudoopthonormal coordinate basis {dx,dy,dz,d^), 
the pseudoorthonormal basis x° = (A, sA*, R, S) and the pseudoorthogonal 
basis X ~ R, S). The matrix x^u of X with respect to the coordinate 

basis is 

u -p 
_ p u 
Xu^^- Q -10' 

1 
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We define now the amplitude of the field by 



\det{Xfj,u)\- 

We consider now the matrix TZ of 2-forms 



(2.21) 



7^ 



udx A d$, —pdx A 



pdy A d^ 





udy A d^ 








—dy A dz 








dz A d^ 



or, equivalently: 

TZ — udx Ad^<S> {dx <S> dx) — pdx Ad^ <S) {dx <S) dy) + pdx Ad^ <^ {dy <S> dx)+ 

+udy Ad^(S) {dy dy) - dy Adz® {dz ® dz) + dz Ad^® {d^ ® d^}- 
Now we can write 



We proceed further to define the phase of the nonlinear solution F . We 
shall need the matrix of the basis x° with respect to the coordinate basis. 
We obtain 



X 



p 








u 











The trace of this matrix is 



Hxl.) 



Obviously, the inequality ||t?"(x^i/)| < 1 is fulfilled. Now, by definition, the 
quantity </? = |^'"(x°i/) will be called phase function of the solution, and the 
quantity 

1 



6 — arccos{ip) — arccos f -tr{x 



(2.22) 



will be called phase of the solution. 

Making use of the amplitude and the phase function (/? we can write 



u 



.LP, p 



l-(^2 



(2.23) 
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We note that the couple of 1-forms A = udx + pdy^ A* — —pdx + udy 
defines a completely integrable Pfafi^ system, i.e. the following equations 
hold: 

dAAAAA* = 0, dA*AAAA* = 0. 

In fact, AaA* = {u^ + p'^)dx Ady, and in every term of dA and dA* at least 
one of the basis vectors dx and dy will participate, so the above exteriour 
products will vanish. 

Remark. These considerations stay in force also for those linear solutions, 
which have zero invariants Ii = I2 = 0. But Maxwell's equations require 
u and p to be running waves, so the corresponding phase functions will be 
also running waves. As we'll see further, the phase functions for nonlinear 
solutions are arbitrary bounded functions. 

We proceed further to define the new and important concept of scale 
factor L for a given nonlinear solution. It is defined by 

\A\ \A*\ 



\SF\ \S*F\ 



(2.24) 



Clearly, L can not be defined for the linear solutions, and in this sense it is 
new and we shall see that it is really important. 

From the expressions F = A A \ and *F = A* A V it follows that the 
physical dimension of A and A* is the same as that of F. We conclude that 
the physical dimension of L coincides with the dimension of the coordinates, 
i.e. [L] — length. From the definition it is seen that L is an invariant 
quantity, and depends on the point, in general. The invariance of L allows to 
define a time-like 1-form (or vector field) /(L)S, where / is some real function 
of L. So, every nonlinear solution determines a time-like vector field on M. 

If the scale factor L, defined by the nonlinear solution F, is a finite and 
constant quantity, we can introduce a characteristic finite time-interval T{F) 
by the relation 

cT(F) = L(F), 
as well, as corresponding characteristic frequency by 

u(F) = 1/T(F). 

In these "wave" terms the scale factor L acquires the sense of "wave length", 
but this interpretation is arbitrary and we shall not make use of it. 

It is clear, that the subclass of nonlinear solutions, which define constant 
scale factors, factors over the admissible values of the invariant f{L). This 
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makes possible to compare with the experiment. For example, at constant 
scale factor L if we choose f{L) = L/c, then the scalar product of (L/c)S with 
the integral energy-momentum vector, which in the F-adapted coordinate 
system is (0, 0, eW, W), gives the invariant quantity W.T, having the physical 
dimension of action, and its numerical value could be easily measured. 



2.3 Nonlinear Solutions. Description of 
photon-like objects 

2.3.1 Explicit solutions in canonical coordinates 

As it was shown in the preceding section with every nonlinear solution F of 
our nonlinear equations a class of F-adapted coordinate systems is associated, 
such that F and *F acquire the form respectively 

F — eudx f\dz -\- udx A + epdy /\dz -\- pdy A d^ 

*F = —pdx Adz — epdx Ad^ + udy Adz + eudy A d^. 
After some elementary calculations we obtain 

5F = (u^ - £Uz)dx + (p^ - epz)dy + e{ux + Py)dz + {u^ + Py)d^, 

S * F ^ -e{p^ - £pz)dx + £{u^ - epz)dy - {px - Uy)dz - {p^ - Uy)d^, 
F^^iSFYdx'' = {*F)^,{S * FYdx'' = 
= e \p{p^ - epz) + u{u^ - euz)] dz + [p{p^ - ep^) + u{u^ - euz)] d^, 
{SF)^ = {S * Ff = -K - euzf - ip^ - epzf 
A simple direct calculation shows, that the equation 

F^,{5*Fy + {*F)^,{5Fy = Q 

is identically fulfilled for any such F with arbitrary u and p. We obtain that 
our equations reduce to only 1 equation, namely 

p{Pi - ePz) + u{u^ - euz) = ^ \{u^ + p')^ - e{u^ + p^)^] = o. (2.25) 
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The obvious solution of this equation is 



u'+p' = ^\x,y,i + ez). (2.26) 

The solution obtained shows that the equations impose some limitations 
only on the amplitude function (j) and the phase function tp is arbitrary exept 
that it is bounded: \(p\ < 1. The amplitude is a running wave along the 
specially chosen coordinate which is common for all F-adapted coordinate 
systems. Considered as a function of the spatial coordinates, the amplitude 
is arbitrary^ so it can be chosen spatially finite. The time-evolution does not 
affect the initial form of 0, so it will stay the same in time. This shows, that 
among the nonlinear solutions of our equations there are (3+1) soliton-like 
solutions. The spatial structure is determined by the initial condition, while 
the phase function (/? can be used to define internal dynamics of the solution. 
Recalling the substitutions (2.23) 

U = (j).(p, p = 0^1 - (p^, 
and the equality \A\ — (f), we get 



\5F\ ^\5.F\^ M^i^, L = (2.27) 
For the induced pseudoorthonormal bases (1-forms and vector fields) we find 



A = (pdx + — 'p'^dy, eA* = 1 — Lp'^dx + (pdy, R = —dz, S = d^, 
d I d I d d d d 

Hence, the nonlinear solutions in canonical coordinates are parametrized 
by one function of 3 parameters and one bounded function of 4 parameters. 
Therefore, the separation of various subclasses of nonlinear solutions is made 
by imposing additional conditions on these two functions. Further in this 
subsection we are going to separate a subclass of solutions , the integral 
properties of which reflect well enough the well known from the experiment 
integral properties and characteristics of the free photons. These solutions 
will be called photon-like and will be separated through imposing additional 
requirements on (p and L in a coordinate-free manner. 
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We note first, that we have three invariant quantities at hand: 0, ip and 
L. The ampHtude function (f) is to be determined by the initial conditions, 
which have to be finite. So, we may impose additional conditions on L and (p. 
These conditions have to express some intra-consistency among the various 
characteristics of the solution. The idea, what kind of intra-consistency to 
use, comes from the observation that the amplitude function is a first 
integral of the vector field V, i.e. 

-e-Q^ + ^1 (0) = -^^'^(^' ^ + + ^'^(^' y^^ + = 0- 

We want to extend this available consistency between V and (p, i.e. we shall 
require the two functions ip and L to be first integrals of some of the available 
vector fields. Explicitly, we require the following: 

1°. The phase function ip is a first integral of the three vector fields A, A* 
and R: A{(p) = 0, A*{ip) = 0, R{(p) = 0. 

2^. The scale factor L is a non-zero finite first integral of the vector field 
S: S(L) = 0. 

The requirement R(v5) = just means that in these coordinates ip does 
not depend on the coordinate z. The two other equations of 1" define the 
following system of differential equations for (p: 

d(p r 2^(p r 2^(p d(p 

Noticing that the matrix 

—(p —VI — 

has non-zero determinant, we conclude that the only solution of the above 
system is the zero-solution: 

^ = ^ = 
dx dy 

We obtain that in the coordinates used the phase function (p depends only 
on ^. Therefore, in view of (2.27), for the scale factor L we get 

m\ 
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Now, the requirement 2°, which in these coordinates reads 



0, 



just means that the scale factor L is a pure constant: L — const. In this way 
we obtain the differential equation 

| = 4yrr^. (2.28) 

The obvious solution to this equation reads 

ip{i) = cos [kj^ + const^ , (2.29) 

where k — ±1. It is worth to note that the characteristic frequency 

U^y (2.30) 

has nothing to do with the frequency in CED. In fact, the quantity L can 
not be defined in Maxwell's theory. 

Finally we note, that the so obtained phase function ip(^) leads to the 
following. The 2-form tr{TZ'^), where 7?.° is the matrix of 2- forms, formed 
similarly to the matrix TZ, but using the basis (A,£A*,R, S) instead of the 
basis (^4, eA*, R, S), is closed. In fact, 

^r(7^°) = ipdx Ad^ + ipdy A - dy A dz + dz A 

and since (p — (p{C): we get dtr(7^°) ~ 0. Note also that the above explicit 
form of tr{TZ^) allows to define the phase function by 



|^r(7^' 



o^|2 



Remark. If one of the two functions u and p, for example p, is equal to zero: 
p = 0, then formally we again have a solution, which may be called linearly 
polarized by obvious reasons. Clearly, the phase function of such solutions 
will be constant: (f — const, so, the corresponding scale factor becomes 
infinitely large: L oo, therefore, condition 2° is not satisfied. The reason 
for this is, that ai p = the function u becomes a running wave and we get 
\6F\ = \6 * F\ = 0, so the scale factor can not be defined by (2.24). 
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2.3.2 Intrinsic angular momentum (spin, helicity) 

The problem for describing the intrinsic angular momentum (lAM), or in 
short helicity, spin of the photon is of fundamental importance in modern 
physics, therefore, wc shall pay a special attention to it. In particular, wc 
are going to consider two approaches for its mathematical description. But 
first, some preceding comments. 

First of all, there is no any doubt that every free photon carries such an 
intrinsic angular momentum. Since the angular momentum is a conserved 
quantity, the existence of the photon's intrinsic angular momentum can be 
easily established and, in fact, its presence has been experimentally proved 
by an immediate observation of its mechanical action and its value has been 
numerically measured. Assuming this is so, we have to understand its origin, 
nature and its entire meaning for the existence and outer relations of those 
natural entities, called shortly photons somewhere in the first quarter of this 
century. 

So, we begin with the assumption: efer?/ free photon carries an intrin- 
sic angular momentum with integral value equal to the Planck's constant h. 
According to our understanding, the photon's lAM comes from an intrinsic 
periodic process. This point of view undoubtedly leads to the notion, that 
photons are not point-like structureless objects, they have a structure, i.e. 
they are extended objects. In fact, according to one of the basic principles 
of physics all free objects move as a whole uniformly. So, if the photon is 
a point-like object any characteristic of a periodic process, e.g. frequency, 
should come from an outside force field, i.e. it can not be free: a free point- 
like (structureless) object can not have the characteristic frequency. 

This simple, but true, conclusion sets the theoretical physics of the first 
quarter of this century faced with a serious dilemma: to keep the notion of 
structurelessness and to associate in a formal way the characteristic frequency 
to the microobjects, or to leave off the notion of structuelesness, to assume 
the notion of extendedness and availability of intrinsically occurring periodic 
process and to build corresponding integral characteristics, determined by 
this periodic process. A look back in time shows that the majority of those 
days physicists had adopted the first approach, which has brought up to life 
quantum mechanics as a computing method , and the dualistic-probabilistic 
interpretation as a philosophical conception. If wc set aside the wide spread 
and intrinsically controversial idea that all microobjects are at the same time 
(point-like) particles and (infinite) waves, and look impartially, in a fair- 
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minded way, at the quantum mechanical wave function for a free particle, we 
see that the only positive consequence of its introduction is the legalization of 
frequency, as an inherent characteristic of the microobject. In fact, the prob- 
abilistic interpretation of the quantum mechanical wave function for a free 
object, obtained as a solution of the free Schroedinger equation, is impossible 
since its square is not an integrable quantity (the integral is infinite). The 
frequency is really needed not because of the dualistic-probabilistic nature 
of microobjects, it is needed because the Planck's relation E = hv turns 
out to be universally true in microphysics, so there is no way to avoid the 
introduction of frequency. The question is, if the introduction of frequency 
necessarily requires some (linear) "wave equation" and the simple complex 
exponentials of the kind const. e,xp[z(k.r — z/t)], i.e. running waves, as "free 
solutions". Our answer to this question is "no". The classical wave is some- 
thing much richer and much more engaging concept, so it hardly worths to 
use it just because of the attribute of frequency. In our opinion, it suffices to 
have a periodic process at hand. 

These considerations made us turn to the soliton-like objects, they realize 
the two features of the microobjects {localized spatial extendedness and time- 
periodicity), simultaneously, and, therefore, seem to be more adequate theo- 
retical models for those microobjects, obeying the Planck's relation E — hu. 
Of course, if we are interested only in the behaviour of the microobject as a 
whole, we can use the point-like notion, but any attempt to give a meaning 
of its integral characteristics without looking for their origin in the consis- 
tent intrinsic dynamics and structure, in our opinion, is not a perspective 
theoretical idea. And the "stumbling point" of such an approach is just the 
availability of an intrinsic mechanical angular momentum, which can not be 
understood as an attribute of a free structureless object. 

Having in view the above considerations, we are going to consider two 
ways to introduce and define the intrinsic angular momentum as a local 
quantity and to obtain, by integration, its integral value. So, these two ap- 
proaches will be of use only for the spatially finite nonlinear solutions of 
our equations. The both approaches introduce in different ways 3-tensors 
(2-covariant and 1-contravariant). Although these two 3-tensors are built of 
quantities, connected in a definite way with the field F, their nature is quite 
different. The first approach is based on an appropriate tensor generaliza- 
tion of the classical Poynting vector. The second approach makes use of the 
concept of torsion, connected with the field F, considered as 1-covariant and 
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1-contravariant tensor. The first approach is pure algebraic, while the sec- 
ond one uses derivatives of F^y. The spatially finite nature of the solutions 
F allows to build corresponding integral conserved quantities, naturally in- 
terpreted as angular momentum. The scale factor L appears as a multiple, 
so these quantities go to infinity for all linear (i.e. for Maxwell's) solutions. 

In the first approach we make use of the scale factor L, the isotropic 
vector field V and the two 1-forms A and A* . By these four quantities we 
build the following 3-tensor H: 

H ^ ®{A^A*). (2.31) 

The connection with the classical vector of Poynting comes through the ex- 
teriour product of A and A*, the 3-dimensional sense of which is just the 
Pointing's vector. In components we have 

H^,^K^Y^{AyAl-A,Al). 
In our system of coordinates we get 

H ^ K- l-e-^ + ^1 (g) {e(j)^dx A dy), 
c \ oz a^J 

so, the only non-zero components are 

c c 
It is easily seen, that the divergence V f_iH^^ — V^i/f2 is equal to 0. In fact, 

" ' 



because (fP' is a running wave along the coordinate z. Since the tangent 
bundle is trivial we may construct the antisymmetric 2-tensor 



Hj,^ = / HAy^dxdydz, 



the constant components of which are conserved quantities. 

L „ W 



H12 = — H21 



f L W 

/ H4 udxdydz = ke—W = keWT = ke — . 

JR? ' c V 
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The non-zero eigen values of Hi^^r are pure imaginary and are equal to ±iWT. 
This tensor has unique non-zero invariant P{F), 



The quantity P{F) will be called Planck's invariant for the finite nonlinear 
solution F. All finite nonlinear solutions Fi, F2, satisfying the condition 



where h is the Planck's constant, will be called further photon-like. The 
tensor field H will be called intrinsic angular momentum tensor and the 
tensor H will be called spin tensor or helicity tensor. The Planck's invariant 
P{F) = WT, having the physical dimension of action, will be called integral 
angular momentum, or just spin or helicity. 

The reasons to use this terminology are quite clear: the time evolution of 
the two mutually orthogonal vector fields A and A* is a rotational-advancing 
motion around and along the 2;-coordinatc (admissible are the right and the 
left rotations: k = ±1) with the advancing velocity of c and the frequency 
of circulation u — c/L. We see the basic role of the two features of the solu- 
tions: their soliton-like character, giving finite value of all integral quantities, 
and their nonlinear character, allowing to define the scale factor L correctly. 
From this point of view the intrinsic angular momentum h of a free photon 
is far from being incomprihensible quantity, connected with the even more 
incomprihensible duality "wave-particle", and it looks as a quite normal in- 
tegral characteristic of a solution, presenting a model of our knowledge of 
the free photon. 

We proceed to the second approach by recalling the definition of torsion 
of two (1,1) tensors. If G and K are 2 such tensors 



S{G, K) (X, Y) = [GX, KY\ + [KX, GY\ + GK[X, Y\ + KG[X, Y]- 
-G[X, KY] - G[KX, Y] - K[X, GY] - K[GX, Y], 




(2.32) 



P(Fi) = P{F2) = ... = h, 




their torsion is defined as a 3-tensor Sf^^, — —S^ by the equation 
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where [, ] is the Lie-bracket of vector fields, 



GX = G':.X''^, GK = G^.Kttdx'' ^ ^ 



and X, Y are 2 arbitrary vector fields. liG — K, in general S{G, G) ^ Q and 

S{G, G){X, r) = 2 {[GX, GY] + GG[X, Y] - G[X, GY] - G[GX, Y]} . 

This last expression defines at every point x E M the torsion S{G, G) — Sq of 
G with respect to the 2-dimensional plane, defined by the two vectors X{x) 
and Y{x). Now we are going to compute the torsion Sp of the nonlinear 
solution F with respect to the intrinsically defined by the two unit vectors 
A and sA* 2-plane. In components we have 



('S'f)L - 2 



dF" 

^ dx<^ dx" " dxi" " dx'^ 

In our coordinate system 

d I d I d d 

ox ^ ay ^ ox oy 

so, 

{SFr^^A^sA*" = {SpYA^'eA*' - A'eA*'). 
For {Sf)i2 we get 



Remark. In our case {Sf)i2 — {S*f)i2, so further we shall work with Sf 
only. 

It is easily seen that the following relation holds: A^eA*^ — A'^eA*^ — 1. 
Now, for the case 

u — (f){x, y, ^ + ez) cos | k-^ + const j , p — (j){x, y, ^ + ez) sin f k-^ + const j 



we obtain 



{Sf)12 = -e{SF)t2 = -Se^^', 
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0,0,-25-^0^2- 



Since 0^ is a running wave along the 2;-coordinate, the vector field Sf{A, eA*) 
has zero divergence: Vi, [Sf{A, eA*)Y — 0. Now we define the helicity vector 
for the solution F by 

^F^^SF{A,eA*). 
zc 

Since L — const, then "Ep has also zero divergence, so the integral quantity 

{'SF)4dxdydz 

does not depend on time and is equal to kWT. The photon-like solutions are 
separated in the same way by the condition WT = h. Here are three more 
integral expressions for the quantity WT. We form the 4-form 

1 K 
— — S A *Ei? = —(p^uJo 
L c 

and integrate it over the 4-volume TZ^ x L, the result is kWT. Besides, we 
verify easily the relations 

-/ \A^A*\uJ, = — ( \5F ^5*F\uJ, = WT. 

C Jb?xL C Jb?xL 

Since we separate the photon-like solutions by the relation WT — h, the 
last expressions suggest the following interpretation of the Planck's constant 
h. Since |A A A*| is proportional to the area of the square, defined by the 
two mutually orthogonal vectors A and eA*, the above integral sums up 
all these areas over the whole 4-volume, occupied by the solution F during 
the intrinsically determined time period T, in which the couple {A^eA*) 
completes a full rotation. The same can be said for the couple {5F, S * F) 
with some different factor in front of the integral. This shows quite clearly 
the "helical" origin of the full energy W = hp oi the single photon. 



2.3.3 Solutions in spherical cordinates 

The so far obtained soliton-like solutions describe objects, "coming from 
infinity" and "going to infinity". Of interest are also soliton like solutions 
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"radiated" from, or "absorbed", by some central "source" and propagating 
radially from or to the center of this source. We are going to show, that our 
equations admit such solutions too. We assume this central source to be a 
small ball with radius To, and put the origin of the coordinate system at 
the center of the source-ball. The standard spherical coordinates (r, ^,</7, ^) 
will be used and all considerations will be carried out in the region out of 
the ball R^. In these coordinates we have 

The *-operator acts in these coordinates as follows: 

*rfr = PsinOdO A d(p A *{dr Ad6 A dcp) = {r"^ sin9)~^ d^ 

*dd = —sinOdr A dip A d^ *{dr A dO A d^) = sinOdcp 

*d(p = {sinOy^dr A dOd^ *{dr A dip A d^) = -{sin9)'^d9 

*d^ = r^sinOdr A dOdip *{d9 A dip A d() = {r^sin9)^^dr 

*{dr A d9) = -sin9dip A d^ *{d9 A dip) = -{r'^sin9)~'^dr A d^ 
*{dr A dip) = {sin9)~^d9 A d^ *{d9 A d^) = -sin9dr A dip 
*{dr A d^) = r'^sin9d9 A dip *{dip A d^ = {sin9)~'^dr A d9. 

We look for solutions of the following kind: 

F = eudr Ad9 + ud9 Ad^ + epdr A dip + pdip A d^, (2.33) 

where u and p are spatially finite functions. We get 

p p 

*F = -dr Ad9 + e -d9 Ad^ - usin9dr A dip - esin9dip A d^. 

sin9 sin9 

The following relations hold: 

F AF = 2e{up - up)dr A d9 A dip A d^ = 0, 

F A*F = \-u^sin9 + u^sin9 - + ) dr A d9 A dip A d^ = 0, 
\ sm9 sin9 ) 

i.e. the two invariants are equal to zero: {*F)^,^F^^'^ — 0, F^^.F^'^ — 0. 
After some elementary computation we obtain 

5FAF^5*FA*F^s[u {spr +P^) - P (sur + u^)] dr Ad9 A dif+ 

+ [u {eur + u^) — u {epr + p^)] d9 Adip A d^, 
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F A *dF = e 



u {eur + u^) sinO + 
u {eUr + u^) sinO + 



P i^Pr + P() 



sinO 
p(spr+p^) 



u {eUr + u^) sin6 + 



sind 

p {epr + p^) 



sinO 



dr Ad9 A dip— 
dO Ad(pA d^, 

dr Ad6 A d(p- 



u {eUr + u^) sinO + 



p {epr + Pi) 
sinO 



d9 Ad(j)A d^. 



So, the two functions u and p have to satisfy the equation 

i , \ ■ Q , Pi^Pr+Pi) n 
u [eUr + u^) smU H ^-p, — ^ = U, 



which is equivalent to the equation 



u sind + 



y SinO J ^ y 

The general solution of this equation is 



sinO 



e u sind + 



(2.34) 



p" 



sinO 



0. 



(2.35) 



u sind + 



P' 



sinO 



(2.36) 



For the non-zero components of the energy-momentum tensor we obtain 



Qi — Qi — Q4 — Q, 



Anr^sinO 



u sinO -\ — : 



p^ 



sin9 



(2.37) 



It is seen that the energy density is not exactly a running wave but when we 
integrate to get the integral energy, the integrand is exactly a running wave: 



W 



An Jr3-ro 



u sind + 



p" 



sinO 



\dr Ade A dcf). 



Since the functions u and p are spatially finite, the integral energy W is finite, 
and from the explicit form of the energy-momentum tensor it follows the well 
known relation between the integral energy and momentum: — c^p^ = 0. 
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2.4 Interference, Nonlinearity and 
Superposition 



2.4.1 Introductory remarks 

Having at hand the photon-hke solutions a natural next step is to try to 
describe the situation when two photons occupy the same (or partially the 
same) 3-region in some period of time. It is clear, that if these two photons 
meet somewhere, i.e. their cilinder-like world-tubes intersect non-trivially, 
the interesting case is when they move along the same spatial straight line 
and in the same direction. Since they move by the same velocities they will 
continue to overlap each other until some outer agent causes a change. What 
kind of an object is obtained in this way, is it a photon or not, what kind of 
interaction takes place, what is its integral energy , its momentum and its 
angular momentum? Many challenging and still not answered questions may 
be set in this direction before the theoretical physics. And this section is 
devoted to consideration of some of these problems in the frame of Extended 
electrodynamics. 

Almost all experiments set to find some immediate mutual interaction 
of two (or more) electromagnetic fields in vacuum, causing some observ- 
able effects (e.g. frequency or amplitude changes), as far as we know, have 
faild, exept when the two fields satisfy the so called coherence conditions. 
In the frame of CED and working with plane waves this simply means, that 
their phase difference must be a constant quantity. The usual way of con- 
sideration is limited to coszne-Iikc running waves with the same frequency. 
The physical explanation is based on the linearity of Maxwell's equations, 
which require any linear combination of solutions to be again a solution, so 
the "building points" of the medium, subject to the field pressure of the 
two independent fields, go out of their equilibrium state obeying simoulta- 
neously the two forces applied in the overlaping 3-region. After getting out 
of this overlaping 3-region the fields stay what they have been before the 
interaction. In order to describe the interaction, i.e. the observed redis- 
tribution of the energy-momentum density inside the overlaping 3-region, 
CED uses the corresponding mathematical expressions in Maxwell's theory 
and gets comparetivcly good results. Most frequently the Poynting vector 
S* ~ [{El + E2) X {Bi + B2)] is used and the cross-terms {EiX B2) + {E2X Bi) 
are held responsible for the interaction, in fact, the very interference is de- 
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fined by the condition that these cross-terms, usually called "interference 
terms", are different from zero. But, the interference takes really place only 
when the coherence conditions are met, while CED permits interference al- 
most always. In other words, from Maxell's theory we can not obtain these 
coherence conditions as necessary conditions for some interference to take 
place. 

In EED we have no superposition principle, so we have to approach this 
physical situation in a new way. First, let's specify the situation more in 
detail and in terms of the notion for E'M-field in EED. Roughly speaking, 
this notion is based on the idea for discreteness, i.e. the real electromagnetic 
fields consist of many noninteracting, or very weakly interacting, photons, 
moving in various directions. Because of the great velocity of their straight 
line motion it is hardly possible to observe and say what happens when two 
photons meet somewhere. The experiment shows that in the most cases they 
pass through each other and forget about the meeting. As we mentioned 
above, the interseting case is when they move along the same direction and 
the regions, they occupy, overlap nontrivially. 

The nonlinear solutions we obtained in the preceding section can not 
describe such set of photons, moving in various directions. Even if we choose 
the amplitude function to consist of many " 3-bubles" (since is a running 
wave, all these bublcs, arc parts of the same running wave), they all have 
to move in the same direction, which is a special, but not the general, case 
of the situation we consider here. So, in order to incorporate for description 
such situations, some perfection of EED is needed. As before, this perfection 
shall consist of two steps: first, elaboration of the algebraic character of the 
mathematical field, second, elaboration of the equations. The second step, 
besides its dynamical task, must define also the necessary conditions for 
interference of photon-like solutions, which should coincide with the above 
mentioned, experimentally established and repeatedly confirmed coherence 
conditions. We shall see that this is easily achievable in EED. 

2.4.2 Elaborating the mathematical object 

Recall that our mathematical object that represents the field is a 2-form 
with values in TZ"^. We want to elaborate it in order to reflect more fully 
the physical situation. The new moment is that inside the 3-region under 
consideration we have many photons. Each of these photons, considered 
as independent objects, is described by a pfoton-like solution as given in the 
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preceding section, i.e. each of them has its own spatial structure, its own scale 
factor (or frequency) and its own direction of motion as a whole. Of course, 
the velocity of motion is the same for all of them. To this physical situation 
we have to juxtapose one mathematical object, which have to generalize in a 
natural way our old object Q. The idea for this generalization is very simple 
and consists in the following. With every single photon, we associate its own 
7^^-space, so if the number of the presenting photons is A^, we'll have N such 
spaces, or a new A^- dimensional vector space. Denoting this vector space by 
jV, our object becomes a 2-form Q, with values in the vector space TZ^ ® H: 
Q, e {M,TZ^ We recall now how this vector space M is exphcitly 

built. 

If /C is an arbitrary set, finite or infinite, we consider those mappings 
of this set into a given field, e.g. TZ, which arc different from zero only for 
finite number of elements of K,. These mappings will be the elements of the 
space M. A basis of this space is built in the following way. We consider the 
elements / e A/", having the property: if a G /C then /(a) = 1 and / has zero 
values for all othe elements of /C. So, with every element a G /C we associate 
the corresponding element /„ G A/", therefore, an arbitrary element / G A/" is 
represented as follows: 

N 
i=l 

where A', i — 1,2, ...,N, are the values, aquired by F, when i runs from 1 to 

(of course, some of the A's may be equal to zero). The linear structure in 
A/" is naturally introduced, making use of the linear structure in TZ in the well 
known way. The linear independence of /„. is easily shown. In fact, assuming 
the opposite, i.e. that there exist such A*, among which at least one is not 
zero and the following relation holds 

N 

1=1 

then for any j = 1,2, ...,N we'll have 

f:A7a,(%)=A^ = 0, 

1=1 

which contradicts the assumption. Hence, fa^ define really a basis of A/". Now 
we form the injective mapping : N ^ Af, defined by 

^iv(a) = fa: a G N, 
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so the set N turns into a basis of N". If such a construction is made, then 
A/" is called a free vector space over the set N. Further onthe corresponding 
basis of our set of photons will be denoted by Ea- So, our mmathematical 
object will look as follows (summing up over the repeating index a) 



O = O'^ £;„ = [F" + {*Fy ® e^] ® Ea, 



(2.38) 



where (e", 63) is the associated with the field basis. If we work in an 
arbitrary basis of 7?.^, the full writing reads {i—1,2) 



Q = Q''(^Ea = Vl'l^kl(» Ea. 



(2.39) 



Following the several times used already method we define the product of two 
2-forms of the kind (2.39). For erxample, if $ = ^'^®k\®Ea, * = '^]®li®Ei,, 
we'll have 

(V, V)($, *) = (V, V)($f A;> Ea, *5 ® li E,) = 

N 
a=l 

+($^ A + A *^) ®kly ®Eay Ea+ 

N 

+ A *i + $1 A ^^) ^k^y ll + ($^ A + $^ A ^>l) ®kly ll+ 



+($? A + A + A + $^ A ^?) ^k'aV lt\ ^EaV Eb. 

Let now be of the kind ^2 = (F" ® + *F"' ® e^) ® F^. Then, forming 
*fi and for (V, V)((5f2, *fi) we obtain 



V, 



N 

, V)(5f}, *n) = J2 [SF" A *F'' ely el + 5 * F" A * * F" ^ elv el+ 

a=l 



+(5F" A**F" + 5*F"A *F") ® V e^J ® F„ V F„+ 

AT 

+ [((^^"A*F''+5F''A*F")(8)e^Ve^+(5*F"A**F''+5*F''A**F")®e^Ve 



a<6=l 



+((5F"A**F'' + 5*F"A*F'' + (5F''A**F" + (5*F''A*F")®e>e^ ^EaVEb. 
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2.4.3 Elaborating the field equations 

If we want to consider a set of independent solutions, then in the above 
expression we take the trace tr over the indeces of EaV Ef,. The compact 
writing of this condition reads 

tr(V, V)(5Q,*Q) = 0, (2.40) 

which is equivalent to the equations 

5F"A*F" = 0, 5*F"A**F" = 0, A * * + 5 * A *F'' = 0. (2.41) 

Clearly, in this case the full energy-momentum tensor Q"^ will be a sum of 
all energy tensors {Q"')'^ of the single solutions. 

The general equations are written down as follows: 

(V, V)(5(],*(]) = 0. (2.42) 

The equivalent (component- wise) form of (2.42) reads 

A *F" = 0, 5 * A * * F" = 0, SF" A * * F" + S * F" A *F" = 0, 

SF" A *F^ + ^F*" A *F" = 0, 5 * F" A * * F*" + 5 * F'' A * * F" = 0, 

(5F" A * * F^ + 5 * F^ A * * F" + 5 * F" A *F^ + 5 * F^ A *F" = 0. 

Let now F", a = 1,2,. ..,N define a solution of the above system of 
equations (2.42). We are going to show that the linear combination with 
constant coefficients Aa 

N 

F = 53AaF" 

a=l 

satisfies the equations: 5F A *F = 0, 5 * F A * * F = 0, 
SF A**F + S*F A*F ^0. In fact 

N N 

5FA*F = ^(A„)2(5F"A*F'^)+ KKiSF" A*F'' + dF'' A*F^), 

a=l a<b=l 
N N 

S*FA**F = ^(A„)2((5*F"A**F")+ J2 A„A6((5*F"A**F''+(5*F''A**F"), 

a=l a<b=l 

N 

SF A**F + S*F A*F ^ Y.iKTiSF" A * * F")+ 

a=l 
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N 



a<b=l 

N 
a<b=l 

N 



N 



a=l 



a=l 



N 



+ Yl A„Ab((5F" A * * F^ + 5F^ A * * F"" + 5 * F'' A *F^ + 5 * F^ A 

a<6=l 

Obviously, the component-wise writing down of the equations (2.42) shows 
that every addend is equal to zero. This result can be interpreted as some 
particular "superposition principle" , i.e. if we have finite number of solutions 
of the system 

5FA*F = 0, 5*FA**F = 0, 6F A * * F + S * F A *F = 0, (2.43) 
which solutions satisfy additionally the equations 

5F" A *F'' + ^F'' A ^F'' = 0, 5 * F" A * * F'' + 5 * F^ A * * F" = 0, (2.44) 

5F" A * * F*" + 5 * F*" A * * F" + 5 * F" A *F'' + 5 * F'' A *F" = 0, (2.45) 

then the 2-form F — J2a=i ^aF"' is again a solution of (2.43). Then, clearly, 
if F and G are 2 solutions of (2.43) and satisfy (2.44) and (2.45), the new 
solution (F + G) of (2.43) is naturally endowed with the following energy- 
momentum tensor 



1 

In the general case we'll have 



-{F + G),,{F + G): 



Q 



An 



\a=l / \a=l / y 



In this way we can compute the corresponding "interference terms". In 
particular, the "interference" energy density is obtained proportional 
to -2F4^G^'^. 
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2.4.4 Coherence and interference 



We consider now two photon- like solutions Fi and F2 of equations (2.43), 
propagating along the same direction. We choose this direction for the ^-axis 
of our coordinate system. We are going to find what additional conditions 
on these solutions come from the additional equations (2.44) and (2.45). 
We assume also, that the 3-regions, where the two amplitudes 0i and 02 
are different from zero have non-empty intersection, because otherwise, the 
interference term is equal to zero. 

Fi — SiUidx Adz + uidx Ad^ + Sipidy Adz + pidy A d^ 
F2 — e2U2dx Adz + U2dx A -|- e2P2dy Adz -\- P2dy A d^, 



where 



ui = (picos [ —^t, + oij , Pi = (pism \ —-^ + ^1 



U2 = 02COS {-^i + ^2^ , P2 = (t)2sin ^^^^ + . 
After an elementary computation we obtain 

5Fi A *F2 + SF2 A *Fi = - -^^^ {uip2 - U2Pi)dx A dy A dz+ 

( K\V\ K2t^2\ , \ , , ,^ 

+ ^^~c~) ^^^^ ~ U2Pi)dx Ady A dt,+ 

+ [Pl (Ulx + Ply) + P2 {U2x + P2y)\ (£1^2 " '^)dx A dz A d^^ 

+ [ui {uix + Ply) + U2 {u2x + P2y)] (1 " £i£2)dy A dz A d^. 

Since 

, , . \fK2l^2 \ ^ , , 1 / n 

U1P2 - U2P1 = 0102'Sm I— —j ^ + 62 - Oi 7^ 0, 

the coefficient before dx Ady A dz will be equal to zero only if kiUi — K2i'2- 
But Vi and z/2 are positive quantities, so it is necessary ki = K2, ui = 1/2- 
Now, the coefficient in front of dx Ady A d^ will become zero if Si = £2- From 
this last relation it follows that the other two coefficients, obviously, are also 
zero. A corresponding computation shows that the so obtained conditions 

ui = f2, Si = £2, 1^1 = 1^2 (2.46) 
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are sufficient for Fi and F2 to satisfy the rest two equations of (2.44) and 
(2.45). Hence, if the 2-form 

Q = (Fi ® ei + *Fi ® 62) ® -El + (F2 ^ki + *F2 (g) ^2) ® E2 

satisfies equations (2.42), then the 2-form, F = Fi + F2 is a solution of our 
initial equations 

5FA*F^0, 5 * F A** F ^0, 5F A * * F + 5 * F A *F ^ 0, 

but not of photon-like type, the coherence conditions (2.42) are satisfied and 
the interference of the two fields Fi and F2 is possible. As for the "interfer- 
ence" energy density we obtain the well known from CED expression 

Wi2 = 01 + 02 + 20102COS(62 " ^l) 

from which the classical interference picture is readily obtained. 
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Chapter 3 

Extended Electrodynamics in Media 



3.1 Basic equations 
3.1.1 Preliminary remarks 

Recall from subsec. (1.5.3) that when we talk about a medium in EED, we 
mean any continuous, i.e. spatially distributed, physical object, exchanging 
energy- momentum with the available in the same region £'M-field Q. For- 
mally, the medium is described by some mathematical object and, when this 
object is chosen, wc talk about external or outer field. When interaction 
between Q and the outer field takes place, of basic importance for the the- 
ory is how Q and the external field participate in the expression, defining 
the exchanged energy-momentum in an unit 4-volume. According to the 
hypotheses we made in subsec. (1.5.3) the E'M- field Q participates directly 
in this expression, while the cxteriour field participates in this expression 
through specially constructed two 7^^-valued 1-forms, and these 1-forms may 
depend on the derivatives of the external field too. The vector-components 
of these 7?.^- valued 1-forms were called currents, since the classical current, 
considered as a 1-form, may be considered as a particular case. Taking into 
account the mathematical model of the £'M-field, as well as the experience 
with Maxwell's theory, we postulated the expression (1.42) 

V($, *7ril)) + V(^', *7r2f)) 
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as sufficiently general and adequate to describe large enough class of exchange 
processes, i.e. interaction of Q with outer fields. 

Such an approach, when the algebraic character and the differential equa- 
tions for the exteriour field are not known, needs a new, general enough and 
sufficiently adaptable viewpoint, expressing a definite comprehension for the 
character of the physical processes considered, as well as general enough and 
adequate enough mathematical facts. Such an adequacy must give reasons 
for definite hypotheses, and, finally, to result in writing down definite equa- 
tions for the currents, no matter what the particular nature of the currents 
is. We note, that we do not require and do not forbid the four currents a* 
to have zero divergence. In our opinion, the new facts to be used as fully as 
possible, are that their number is more than one , that every current realizes 
a separate energy exchange channel, and that there should be some correla- 
tions among them. It is naturally to expect such correlations among the four 
currents to exist since the exchange processes occur locally and on the other 
side of this exchange stays just one physical object - the £^Af -field Q. These 
correlations must be organized in such a way, that to incorporate the special 
case of only one current different from zero, as it is in CED. Of course, we 
have to remember that in CED the current is a vector field, while our currents 
are 1-forms. These are different objects although the available isomorphism 
through the pseudometric rj, and this difference will be explicitly taken into 
account in our approach. 

3.1.2 Maxwell's theory as a particular case 

Let's recall some important features of CED. First, the EM-field exchanges 
energy-momentum only through F and does not exchange energy-momentum 
through *F. Formally, this is accounted by the equation S * F = 0, which 
contains (in differential form) the Faraday's induction law. In other words, 
CED assumes, that this experimentally established fact for some media, holds 
for all media. As we know, this assumption legalizes the C/(l)-gauge inter- 
pretation of CED, where the equation S * F = acquires the geometrical 
interpretation of Bianchi's identity for the abelean group U{1). 

Second, an energy-momentum exchange may be realized only with free or 
bound electric charges. The formal description of this exchange was explained 
and commented in subsec. (1.1.3). The second equation of CED, 

SF = 4TT{j + jh), 
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identifies the pure field quantity 5F witli tlie outer quantity current, which 
characterizes the distribution and mechanical behaviour of the charge carri- 
ers. To what extent such an identification is admissible is a personal view, 
and we are not going to comment it. According to us much more natural is 
to write down a relationship having the sense of local energy-momentum bal- 
ance. In other words, the same (exchanged) quantity of energy-momentum 
to be written down in two ways: in the first way, by means of the components 
F^y and their derivatives only, and in the second way, through expressions 
where the outer field components necessarily take part. Then, according to 
the local conservation law, the two expressions are equalized. Namely such 
an approach we have realized in EED. 

The third feature we recall is the lack of some general enough and com- 
mon approach for determination of the full current (j/ree + jbound) ■ As we 
mentioned earlier, the series developments of the polarization vector P and 
magnetization vector A4 (or the 2-form S) with respect to E and B and 
their derivatives may be a felicitous working skill, but it is not a perspective 
theoretical idea. For example, this approach is not applicable for strongly 
nonhomogeneous media, while a local energy-momentum balance equation 
can be written down always and to be particularized and made more precise 
in the course of work. 

We see now how, from formal point of view, CED is incorporated in EED. 
The mathematical expression for the exchanged energy-momentum in CED 
is 

47rFf,^{SS+jydxi'. 

In our approach, using the defined by the field fl basis (ei, 62), this expression 
is represented consecutively as follows: 

A-n-Fi^^iSS + jydx^" (g) Ci V ei = *[A'k{5S + J) A *F] O ei V ei = 

= *V [47r((5S' + J) (g)ei,*F(8)ei] = * V [tti*, tti * fi], 

where 7ri$ = 4:n{6S + J) ® Ci. In this way, using the notations of EED, we 
get 

y(5n, *n) = V(7ri$, 7ri*n) + V(7ri^', 7:2*^), ni^ = TTi*. 

We see that there are no terms, describing an exchange through *F, and the 
redistribution energy-momentum equation reduces to one of the real exchange 
equations. So, CED reflects the following additional requirements to the 
equations of EED: $ = (g) ei, ^' = (g) ei, = = 47r{SS + J), = 
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a'^ = 0. The very Maxwell's equations S * F = 0, 5F ^ A7r(5S + J) may 
written also as: 7r2(50 = 0, ni5il — A7r{5S + J) ei. 

3.1.3 Component form of the equations 

The coordinate free written relationship (1.43) 

V {Sfl, *n) =\/{^,*nin) + \/{-^,*n2^}) 

is equivalent to the following relations: 

5F A*F ^a^ A*F, 5*F A**F ^a^ A**F, (3 
5FA**F + 5*FA*F^a^A**F + a^A*F, 

or, in components 

F^,(5F)- = {*F)^,{5 * FY = (*F)^,(a4)^ (3 

Moving everything on the left, we get 

{SF - a^) A *F = 0, (5 * F - a^) A * * F = 0, 

(5F - a^) A * * F + (5 * F - a^) A *F = 0, 
or in components 

F^,{5F - a^y = 0, {*F)^,{5 * F - a^)'^ = 0, 

F,,{5 * F - a^y + (*F)^,(5F - a^)- ^ q. 
Summing up the two equations 

we obtain 

F,,{5Fy + {*F),,{5 * Fy = v.g;: = F^,{a'y + (*f)^,(«^)^ 

This relation shows that the sum 

F^Ac^'y + {*F)^Aa^y 
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is a divergence of a 2-tensor, which we denote by — -P^- In this way we obtain 
the local conservation law 

v.(g; + p;) = 0. (3.3) 

Thus, we get the possibility to introduce the full energy-momentum tensor 

where is interpreted as interaction energy-momentum tensor. Clearly, 
can not be determined uniquely in this way. 

So, according to (3.2), for the 22 functions F^,^, (a*)^ we have 12 equa- 
tions, and these 12 equations are differential with respect to F^^ and algebraic 
with respect to (a*)//. Our purpose now is to try to write down differential 
equations for the components of the 4 currents. The leading idea in pur- 
suing this goal will be to establish a correspondence between the physical 
concept of nan- dissipation and the mathematical concept of integrability of 
Pfaff systems. The suggestion to look for such a correspondence comes from 
the following considerations. 

Recall from the theory of the ordinary differential equations (or vector 
fields), that every solution of a system of ordinary differential equations 
(ODE) defines a local (with respect to the parameter on the trajectory) group 
of transformations, frequently called local flow. This means, in particular, 
that the motion along the trajectory is admissible in the two directions: we 
have a reversible phenomenon, which has the physical interpretation of lack of 
losses (energy-momentum losses are meant). Assuming this system of ODE 
describes fully the process of motion of a small piece of matter (particle), we 
assume at the same time, that all energy-momentum exchanges between the 
particle and the outer field are taken into account, i.e. we have assumed that 
there is no dissipation. In other words, the physical assumption for the lack 
of dissipation is mathematically expressed by the existence of solution - local 
flow, having definite group properties. The existence of such a local flow is 
guaranteed by the corresponding theorem for existence and uniqueness of a 
solution at given initial conditions. This correspondence between the mathe- 
matical fact integrability and the physical fact lack of dissipation in the simple 
case "motion of a particle" , we want to generalize in an appropriate way, 
having in view possible applications in more complicated physical systems, 
in particular, the physical situation we are going to describe: interaction of 
the field VL with some outer field, represented in the exchange process by the 
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four 1-forms a*. This will allow to write down equations for a* in a direct 
way. Of course, in the real world there is always dissipation, and following 
this idea we are going to take into account its neglecting as conditions (i.e. 
equations) on the currents. As it is well known, the mathematicians have 
made serious steps towards studding and formulation of criteria for integra- 
bility of partial differential equations, so it looks unreasonable to close eyes 
before the available and represented in an appropriate form mathematical 
results. 

On the other hand it is interesting, and may be suggesting, to note the 
following . In physics we have two universal things: dissipation and gravita- 
tion. We are going now to establish a correspondence between the physical 
notion of dissipation and the mathematical concept of non-integrability. As 
we know, the mathematical non-integrability is measured by the concept of 
curvature. General theory of Relativity describes gravitation by means of 
Riemannian curvature. The circle will be closed if we connect the universal 
property of any real physical process to dissipate energy-momentum with the 
only known so far universal interaction in nature, the gravitation. 



3.2 The Frobenius Integrability and 
Dissipation 

3.2.1 Integrable distributions and curvature 

The problem for integration of a system of partial differential equations of 
the kind 

^ = f^(x',y'), t,k = l,...,p; a,b=l,...,q, (3.4) 

where f^{x^,y'') are given functions, obeying some definite smoothness con- 
ditions, has brought about to the formulation of a number of concepts, which 
in turn have become generators of ideas and directions, as well as have shown 
an wide applicability in many branches of modern mathematics. A particular 
case of the above system (nonlinear in general) of equations is when there is 
only one independent variable, i.e. when all ,t* are reduced to x^, which is 
usually denoted by t and the system acquires the form 

^ = r{y\t), a,b = l,...,q. (3.5) 
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We are going to give now the system of concepts used in considering the 
integrabihty problems for the equations (3.4) and (3.5), making use of the 
geometric language of manifolds theory. 

Let X be a vector field on the g-dimensional manifold M and the map 
c : I ^ M, where / is an open interval in TZ, defines a smooth curve in M. 
Then if X" arc the components of X with respect to the local coordinates 
{y^, ■■■,y'^) and the equality c'{t) — V{c{t)) holds for every i e /, or in local 
coordinates, 

f - ^°(^')- 

c{t) is called integral curve of the vector field X. As it is seen, the differ- 
ence with (3.5) is in the additional dependence of the right side of (3.5) on 
the independent variable t. Mathematics approaches these situations in an 
unified way as follows. The product TZx M is considered and the important 
theorem for uniqueness and existence of a solution is proved: For every point 
p e M and point t there exist a vicinity U oi p, a positive number e and 
a smooth map $ : (r — e, r + e) x f/ M, $ : (t, y) ^tiv), such that for 
every point y & U the following conditions are met: i^riy) = V, t ^ Vtiy) is 
an integral curve of X, passing through the point y G M; besides, if two such 
integral curves of X have at least one common point, they coincide. More- 
over, if (t', I/), {t-\-t\ y) and {t, '^{y)) are points of a vicinity U' of {0} x 7^ in 
TZ X M, we have ipt+t'{y) = Vt{'^t'{y))- This last relation gives the local group 
action: for every t E I we have the local diffeomorphism (pf : U ^ Vt{U). So, 
through every point of M there passes only one trajectory of X and in this 
way the manifold M is foliated to non-crossing trajectories - 1-dimensional 
manifolds, and these 1-dimensional manifolds define all trajectories of the 
defined by the vector field X system of ODE. This fibering of M to non- 
intersecting submanifolds, the union of which gives the whole manifold M, 
together with uniting t and y{t) in one manifold, is the leading idea in treat- 
ing the system of partial differential equations (3.4), where the number of 
the independent variables is more than 1, but finite. For example, if we con- 
sider two vector fields on M, then through every point of M two trajectories 
will pass and the question: when a 2-dimcnsional surface, passing through a 
given point can be built, and such that the representatives of the two vector 
fields at every point of this 2-surface to be tangent to the surface, naturally 
arises. Mathematics sets this question for p-dimensional surfaces, builds the 
necessary concepts and proves the corresponding theorems. These theorems 
formulate the criteria for integrabihty of (3.4), and are known in the litera- 
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ture as Frobenius theorems. For simplicity, further we consider regions of the 
space TZ^ x 7^^, but this is not essentially important since the Frobenius theo- 
rems are local statements, so the results will hold for any (p-h g) -dimensional 
manifold. 

Let [/ be a region in W x 7^^ and (x^ x^, = = xP^") are 

the canonical coordinates. We set the question: for which points (a;o,yo) of 
U the system of equations (3.4) has a solution = Lp'^{x''), defined for points 
X, sufficiently close to Xq and satisfying the initial condition ip{xq) = y^i The 
answer to this question is: for this it is necessary and sufficient the functions 
/f on the right hand side of (3.4) to satisfy the following conditions: 

-^(^^ y) + y)-fA^^ y) - -^(^^ y) + y)-f'i^^ y)- (3-6) 

This relation is obtained as a consequence of two basic steps: first, equalizing 
the mixed partial derivatives of with respect to and x^ , second, replacing 
the obtained first derivatives of with respect to on the right hand side of 
(3.4) again from the system (3.4). This second step means, that everywhere 
in (3.6) y are considered as functions of i.e. there is no explicit dependence 
on y. If the functions /f satisfy the equations (3.6), the system (3.4) is called 
completely integrable. In order to give a coordinate free formulation of (3.6) 
and to introduce the curvature, as a measure for non-integrability of (3.4), 
we shall first sketch the necessary terminology. 

Let M be an arbitrary n = p + q dimensional manifold. At every point 
X & M the tangent space Tx{M) is defined. The union of all these spaces with 
respect to the points of M defines the tangent bundle. On the other hand, 
the union of the co-tangent spaces T*{M) defines the co-tangent bundle. 
At every point now of M we separate a p dimensional subspace A^(M) 
of Te(M) in a smooth way, i.e. the map x A^; ix smooth. If this is 
done we say that a p-dimensional distribution A on M is defined. From the 
elementary linear algebra we know that every p-dimensional subspace A^, of 
Tj,(M) defines unique {n — p) = q dimensional subspace A* of the dual to 
Tx{M) space T*{M), such that all elements of A* annihilate (i.e. send to 
zero) all elements of A.^-. In this way we get a g-dimcnsional co- distribution 
A* on M. We consider those vector fields, the representatives of which 
at every point are elements of the distribution A, and those 1-forms, the 
representatives of which at every point are elements of the co-distribution 
A*. Clearly, every system of p independent vector fields, belonging to A, 
defines A equally well, and in this case we call such a system a differential 
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p-system V on M. The corresponding system V* of q independent 1-forms 
is called g- dimensional Pfaff system. Clearly, ii a & V* and X E V, then 
a{X) = 0. 

Similarly to the integral curves of vector fields, the concept of integral 
manifold of a p-dimensional differential system is introduced. Namely, a 
p-dimensional submanifold of M is called integral manifold for the p- 
dimensional differential system V, or for the j9-dimensional distribution A, to 
which V belongs, if the tangent spaces of at every point coincide with the 
subspaces of the distribution A at this point. In this case is called also 
integral manifold for the gr-dimensional Pfaff system V*. If through every 
point of M there passes an integral manifold for V, then V and V* are called 
completely integrable. 

Now we shall formulate the Frobenius theorems for integrability. 

A differential system V is completely integrable if and only if the Lie 
bracket of any two vector fields, belonging to V, also belongs to V . 

So, if generate the completely integrable differential system 

7^, then 

{X,,X,\ = C%X^, (3.7) 

where the coefficients depend on the point. 

This criterion is not quite convenient to use because its usage presupposes 
the knowledge of the functions C\y It turns out that the corresponding 
criterion for Pfaff systems does not require any additional information. In 
fact, let the 1-forms (a^, a'^) define the g-dimensional Pfaff system V* . 
The following criterion holds (the dual Frobenius theorem): 

The Pfaff system V* is completely integrable if and only if 

^a" = Kl^ct A a", h< c. (3.8) 

It is easily shown, that the above equations are equivalent to the following 
equations: 

{da"-) A A ... A a" A . . . A = 0, a = 1, q, (3.9) 
which, obviously, do not depend on any coefficients. 
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When a given Pfaff system V*, or the corresponding differential system 
V, are not integrable, then the relations (3.7)-(3.9) are not fulfilled. From 
formal point of view this means that there is at least one couple of vector fields 
(X, Y), belonging to V, such that the Lie bracket [X, Y] does not belong to V. 
Therefore, if at the corresponding point x E M we choose a basis of Tx{M) 
such, that the first p basis vectors to form a basis of A^, then [X, F] will 
have nonzero components with respect to those basis vectors, which belong 
to some complimentary to A^. subspace : T^^M) = A^; ©F^. If the spaces 
Fj.,^ e M are beforehand given we may consider the projection of the Lie 
bracket [X, Y] onto these comphmentary subspaces. Then this projection is 
defined uniquely by the choice of the distribution A, so it is a natural measure 
for the non-integrability of A. If M has the structure of bundle space, which 
means that a base-space B, diniB = p, is given, and a smooth map tt : M —>■ 
B of maximal rank, i.e. rank{dn) = p < n, is given, then the subspaces T4 = 
Ker{dn)a; C T^{M) are naturally separated. It is easily shown that these 
subspaces. called usually vertical, form an integrable distribution. If wc orient 
our interest towards distributions A(M), which are complimentary to vertical 
distributions and are usually called horizontal, then the non-integrability of 
A(M) will be determined entirely by the vertical projection v : T{M) 
V{M), defined by the definition of A(M) and considered on the various Lie 
brackets of horizontal vector fields. It is clear now, that the curvature JC of 
the horizontal distribution A is defined by 

}C{X, Y) = v{[X, Y]), X,Y eA. (3.10) 

It is also clear, that the curvature /C is a 2-form on M with values in the 
tangent bundle T{M), and IC is reduced to identity on Lie brackets of vertical 
vectors. In fact, if / is a smooth function and X, Y are two horizontal vector 
fields, then 

}C{X, fY) = v{\X, fY]) = v{f[X, Y] + XU)Y) = 

= fv{[X, Y]) + X{f)v{Y) = flC{X, Y) 

because v{Y) = 0, F — horizontal. 

The corresponding g-dimensional co-distribution A* (or Pfaff system V*) 
is defined locally by the 1-forms (9^, ...,9''), such that 9°'{X) = for all 
horizontal vector fields X. In this case the 1-forms 9"' are called vertical, and 
clearly, they depend on the choice of the horizontal distribution A. In these 
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terms the non-integrability of A means 



^ Klf A r, h< c. 

This non-equahty means that at least one of the 2-forms dO"- is not vertical, 
i.e. it has a nonzero horizontal projection H*d6'^, which means that it does 
not annihilate all horizontal vectors. In these terms it is naturally to define 
the curvature by 

//Mr = - K^^e'' AO", b< c. 

Further we shall see how this picture is defined by the equations (3.4). 

Let's begin with the remark that the consideration of bundle spaces only, 
is not a limitation and does not bounds the results, because we want through 
every point of M to pass locally only one integral manifold of a given hor- 
izontal distribution and this integral manifold will be diffeomorphic to an 
open set in the base manifold B. And this is just what is guaranteed by the 
bundle structure of M: every point b & B has an open vicinity U, and 7i~^{U) 
is diffeomorphic to the direct product U x N, where is a g-dimensional 
manifold, called standard fiber. This bundle structure allows canonical (or 
bundle adapted) local coordinates (x',y") to be introduced, reflecting the 
local-product nature of M: = o tt, where z'^ are local coordinates on 
U C B, and y"' are local coordinates on N. In these coordinates a local basis 
of the vertical vector fields is given by 

9 

— , a = p + l,...,p + g = n. 

Now, making use the equations of the system (3.4), i.e. the functions /f , we 
have to define the horizontal spaces at every point of it~^{U) C M, i.e. local 
linearly independent vector fields Xi,i = 1, ...,p. The definition is 

^. = |T + /f^. (3.11) 

The corresponding Pfaff system shall consist of 1-forms 9'^,a — p+1, ...,p + q 
and is defined by 

O"" = dy"" - f^dx\ (3.12) 

In fact. 
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^0 + fi5^-f^5i-0^0. 

Remark. The coordinate 1-forms dy"' are not vertical with respect to the so 
defined horizontal distribution. 

In this way in the bundle-adapted coordinate systems the system (3.4) 
defines unique horizontal distribution. On the other hand, if a horizontal 
distribution is given and the corresponding vertical Pfaff system admits at 
least i basis, then this basis may be chosen of the kind (3.12) always. More- 
over, in this kind it is unique. In fact, let {6'^, ...,9"^) be any local basis of 
A*. Then in the adapted coordinates we'll have 

where Al, Bf are functions on M. We shall show that the matrix Al has 
non-zero determinant, i.e. it is non-degenerate. Assuming the opposite, we 
could find scalars A^, not all of which are equal to zero, and such that the 
equality XaAl = holds. We multiply now the above equality by A" and sum 
up with respect to a. We get 

Note now that on the right hand side of this last relation we have a hori- 
zontal 1-form, while on the left hand side we have a vertical 1-form. This 
is impossible by construction, so our assumption is not true, i.e. the inverse 
matrix {A^)~^ exists, so multiplying on the left O'"' by {A^)~^ and putting 
(^a)-i^/6 ^ Qa obtain 

e'' = dy''+{Al)-'^B^dx\ 

We denote now {A^)^^B'^ = — /f and get what we need. The uniqueness part 
of the assertion is proved as follows. Assume there is another basis {6'^)" of 
the same kind. So, there must be a non-degenerate matrix C^, such that 
(^«)" = C^eK We get 

dy'' - {ft)"dx' = Ctdy' - Ctf,dx\ 

and from this relation it follows that the matrix is the unit one. 

Let's see now the explicit relation between the integrability condition 
(3.6) of the system (3.4) and the curvature if*d^" of the defined by this 
system horizontal distribution. We obtain 

a fa a fa 

de" = -^dx' A dx' - ^dy'' A dx\ 
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In order to define the horizontal projection of d^" it is necessary to separate 
the horizontal part of dy"". Since O"" are vertical, from their explicit form 
is seen that the horizontal part of dy°' is just f^dx^. That's why for the 
curvature H*d9°^ we get 



dx^ dx^ ' dy'' '' •' dy^ * 



dx'' A da^ , i < j. 



(3.13) 



It is clearly seen that the integrability condition (3.6) coincides with the 
requirement for zero curvature. The replacing of dy^, making use of the 
system (3.4), in order to obtain (3.6), acquires now the status of "horizontal 
projection" . 

We verify now that the curvature, defined by v{[Xi,Xj]) gives the same 
result. 

— + r— — + f^— 




dx^ dy^ dx^ dy°- 



dy"- dx^ 

' J i Ci r, C\ h J i 



J ■ 



dx^ dx^ dy 



Qya Qyh 



J Qyb 

Qyh Qya ^ 

d 



b •'1 



Qyb 3 J Qya 



The interchange of the indices i and j does not impact the equivalence of 
this result to the above obtained for H*d9°'. 



3.2.2 Physical interpretation 

As it was noticed in the preceding subsection we want to connect the phys- 
ical concept of dissipation with the mathematical concept of Frobenius non- 
integrability of Pfajf systems. The availability of a well defined mathematical 
quantity as the curvature, which has been extensively used from the begin- 
ning of this century in mathematics (differential geometry and differential 
topology) and theoretical physics (General Relativity (GR) and Gauge the- 
ories) makes the things more attractive in view of its wider use in physics. 
In General Relativity curvatures of Riemanncan connections are used, and 
because of the stress on the metric tensor as a potential of the curvature, 
there has not been paid enough attention to the original meaning of the cur- 
vature, namely as a measure for non-integrability. Moreover, the definitive 
and physically not motivated assumption of the Riemannean curvature as a 
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mathematical adequate of the gravitational field in GR does not contribute 
to a full comprehension of why the computations in the theory meet the 
experiment in the Solar system (and even out of it) so well. This circum- 
stance, being so charming in the early days of the theory, may generate some 
hesitations, because 80 years seem to be long enough time for the clarifi- 
cation of this fundamental for GR problem. Together with the well known 
"energy-momentum problems", this may lower the authority of GR. 

In gauge theories the curvature, considered as generated by connections 
on principle bundles for some groups and their representations, leading to 
linear connections in vector bundles, is also a leading concept. For exam- 
ple, the energy-momentum tensor in these theories is a quadratic expression 
of the curvature, and the significance of energy and momentum in modern 
classical and quantum gauge theories and in the whole physics at all is out 
of any doubt. Except electrodynamics, where we have much enough expe- 
rience, in other gauge theories there is also no enough motivation for using 
namely the curvature as a mathematical adequate of a physical field. The 
consideration of connections as basic mathematical objects in physical the- 
ories we do not consider as sufficiently legalizing move for the introduction 
of curvatures, although the connections define derivative laws for the sec- 
tions of vector bundles. In our opinion, a more basic analysis of the question 
why the curvature works well in physical theories from the point of view of 
the Frobenius integrability theory would contribute to a more complete and 
detailed understanding of this important moment in the field theories. 

In mathematics the curvature defines those conditions, at which, given 
differential relations determine integral objects, or what are the obstacles for 
building these integral objects. In physics, from Newton's time on, a basic 
quantity is the force, i.e. the quantity of transferred energy-momentum from 
the outer field to the object under observation (usually test particle (s)) in an 
unit 4-volume. From the point of view of the outer field this means loss of 
energy-momentum, and that's why it does not seem reasonable to expect a 
Frobenius integrability for the equations, describing the outer field, if these 
equations do not take into account these losses. It is quite illogical and 
unreasonable to expect that the expressed through differential relationships 
properties of a physical object should define it entirely (or integrally) if there 
is some "flowing out" of substance towards other physical objects, if this 
"flowing out" takes some energy- momentum from it and carries it to the 
other objects. In other words, the interaction, if it is not fully described, may 
violate the existing for any extended object connection between its differential 
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(local) and integral properties. And a fully described interaction means to 
say, and to take into account mathematically, where the energy-momentum 
losses go to. These losses are just the force, applied onto the other object(s). 

The curvature of a given differential system V measures mathematically 
(at a given point and in a differential way) something very close to a "flowing 
out" of V, namely, it determines what parts of the Lie brackets of vector 
fields in V belong no more to V. We may say, that these parts of the Lie 
brackets, namely the vertical projections of the Lie brackets, define local 
flows, directed out of V. Therefore, from physical point of view it is natural 
to choose curvatures as measuring quantities of the external exertion on test 
particles in outer fields. The physical measure of this exertion is a change in 
the universally conserved quantities energy and momentum (of the particle). 

As far as test (point-like) particles in outer fields are considered the us- 
age of curvature as a direct participant in describing the energy-momentum 
transfer seems natural, since we are not interested in the dynamics of the ex- 
ternal fields. " No change in the energy- momentum of the particle" is equal 
to "no presence of external fields", i.e. the force is zero. This approach works 
no more in case of local interaction of two (or more) continuous objects, i.e. 
two fields, where together with the local energy-momentum transfer we are 
interested also in the dynamics of the two interacting fields. So, the local 
dynamical changes of any of the two interacting fields should depend on 2 
things: the proper dynamical character of any of the fields and the kind of 
interaction. These two components of the physical system should be consis- 
tent, and when there is no interaction, the mathematical expression for zero 
energy-momentum transfer should also be consistent with the proper dynam- 
ics of any of the fields, considered now as free. Because of the local character 
of the interaction, derivatives of the fields' components participate in the 
corresponding mathematical expression. In this way, putting the interaction 
expression equal to zero, we obtain general enough differential equations for 
the free fields. And if the fields are mathematically described by curvatures 
we obtain differential equations for the curvatures. 

Let's consider now a more complicated situation, when the interaction 
energy-momentum flows along several channels. Every such channel may 
or may not produce dissipation. The produced by some of the channels 
dissipation of energy-momentum could be utilized or not by some of the rest 
channels. Our proposal to interpret mathematically the non-zero dissipation 
by Frobenius non-integrability would imply in such cases, that some of the 
subdistributions of the initial distribution may be nonintegrable, but these 
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nonintegrable subdistributions may participate in integrable subdistributions 
of a higher dimension. The same can be said about the corresponding Pfaff 
systems. This estabhshes some kind of hierarchy among the subsystems of 
the Pfaff systems and leads to a more complete study of the initially assumed 
Pfaff system, defining the energy-momentum interaction channels. 

We shall see in the next section that the considered examples-solutions of 
our equations, representing a (3+l)-intcrprctation of all known (l-l-l)-soliton 
solutions illustrate the above outlined idea: availability of differential and in- 
tegral conservation laws, non-integrabihty of 1-dimensional Pfaff subsystem, 
integrability of all 2-dimensional Pfaff systems. This is in the general spirit 
of our approach, in which the leading ideas are the extended character of the 
real objects and the interrelation between their local and integral properties. 

We are not going to consider here non-utilized dissipations at the highest 
possible level. This would lead us far behind our purposes. 

After these mathematical and interpretational deviations we go back to 
EED. First we note that our base manifold, where all fields and operations 
are defined, is the simple 4-dimensional Minkowski space. According to our 
equations (1-43) the medium reacts to the field fl by means of the two 71^- 
valued l-forms $ and ^. So, we obtain four 7?.- valued 1-forms a^, a^, a^, a^. 
Because of the 4-dimensions of Minkowski space and the kind of the equations 
(3.9) it is easily seen that only 1-dimcnsional and 2-dimcnsional Pfaff systems 
may be of interest from the Frobenius integrability point of view. All Pfaff 
systems of higher dimension are trivially integrable. Note also that closed 1- 
forms and linearly dependent 1-forms define always integrable Pfaff systems 
in an obvious way. 

The integrabihty equations for 1-dimensional Pfaff systems are 



Every of the 4 equations (3.14) is equivalent to 4 scalar nonlinear equations 
for the components of the corresponding 1-form. We note, that the solutions 
of (3.14), as well as the solutions of the general integrability equations for 
a p-dimensional Pfaff system are determined up to a scalar multiplier, i.e. 
if a* are solutions, then /j.a* (no summation over i), where fi are smooth 
functions, are also solutions. 

In case of 2-dimensional Pfaff systems {a\a^), defined by four 1-forms, 
their maximal number is 4.3=12. The Frobenius equations read 



dQ;*Aa' = 0, i = 1,2,3,4. 



(3.14) 



da' A a* A q;-^ = 0, i ^ j. 



(3.15) 
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We have here 12 nonUnear equations for the 16 components of a*. Clearly, 
these equations (3.15) make some interest only if the corresponding «*, the 
exteriour differential da* of which participates in (3.15), does not satisfy 
(3.14). 

After these general remarks we pass to finding exphcit solutions of (1.43) 
with non-zero currents. 

3.3 Explicit Solutions with Nonzero Currents 

3.3.1 Choice of the anzatz and finding the solutions 

From purely formal point of view finding of a solution, whatever it is, le- 
gitimizes the equations (1.43) and (3.15) as a consistent system. Our pur- 
pose, however, is not purely formal, we consider it as physically meaningful, 
namely, we are interested in solutions, which are physically interpretable as 
models of real objects in the above commented sense. That's why we have to 
meet the following. First, the solutions must be physically clear, which means 
that the anzatz assumed should be comparatively simple and its choice should 
be made on the base of a preliminary analysis of the physical situation in 
view of the mathematical model used. Second, it is absolutely obligatory the 
solutions to have well defined integral energy and momentum. Third, to be in 
the spirit of the soliton-like comprehension of the real natural objects when 
the solution is interpreted as a model of such an object. Fourth, to realize 
the above given physical interpretation of the Frobenius integrability equa- 
tions and the mentioned dimensional hierarchy of Pfaff systems as a possible 
model of suitably chosen and intrinsically structured interrelated processes. 
Finally, it would be nice, the solutions found to be comparatively simple and 
interesting as corresponding generalizations, or extensions, of "popular" and 
well known solutions of "well liked" equations. Probably not all solutions 
will satisfy these requirements, but there must be such solutions, since the 
physical significance of our equations depends strongly on this circumstance, 
to admit solutions with the above mentioned features. Let's now get started. 

The first, we take into account, is the necessary time- dependence of the 
solutions, therefore, the "electric" and the "magnetic" components should 
present. The simplest Vt, or {F,*F), meeting this requirement, looks as 
follows (we use the above assumed notations): 

F = —udy Adz — vdy A d^, *F = vdx Adz + udx A d^. (3.16) 
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When choosing the 1 -forms a* we shall obey the requirement, that the 
"medium" does not involve in the F *F exchange, so we put 

q;2 = q;3 ^ 0. (3.17) 

The rest two 1-forms , and must be linearly independent. This would 
be guaranteed if we choose one of them to be time-like., and the other to 
be space-like. The simplest time-like 1-form is of course a — y, 2;, 
moreover, any such 1-form defines integrable 1-dimensional Pfaff system: 

da A a = {A^dx A d^ + Aydy A d^ + A^dz A d^) A Ad^ = 0. 

The choice of the last one, denoted by /3, reads (3 = Pidx + P2dy + Psdz. 
Clearly, < 0, and a. (3 = 0, so they are independent. Now, j3 participates 
in the equations through the expression (3 A *F, and from the explicit form 
of *F is seen, that the coefficient Pi in front of dx does not take part in the 
equations, therefore, we put /3i = 0. So, we get 

= l3 = hdy- Bdz, = a = Ad^, o? = o? = 0. (3.18) 

We note, that the so chosen /3 does not define in general an integrable 1- 
dimensional Pfaff system, so the requirement for hierarchy of the exchange 
processes, i.e. daAa = 0, d/3A/3 7^ 0, d/^A/^Aa = 0, is obeyed at a definite 
level. 

At these conditions our equations 

(5*FAF = aAF, (5FA*F = /?A*F, 5*FA*F-(5FAF = 0, 

dQ;AQ;A/3 = 0, d/3AQ;A/3 = 
take the form: 5*FA*F — 5FAF = 0is reduced to 

—VUy -|- UVy — 0, —UVx + VUx — 0, 

the Probenius equations dQ;AaA/5 = 0, d/3AQ;A/? = reduce to 

{-h^B + BJ)) .A = 0, 
5*FAF — aAF\s reduced to 

U {u^ — Vz) = 0, V {u^ — Vz) = 0, UUx — VVx = Au, 
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SF A*F ^ P A*F reduces to 

V {v^ — Uz) = —bv, u {v^ — Uz) 



-bu, uuy — vvy = Bu. 



In this way we obtain 7 equations for 5 unknown functions u,v, A, B,b. The 



equations —vUy + uVy = 0, —UV:^ 



have the following solution: 



u{x,y,z,C) = f{x,y)U{z,C), v{x,y,z,C) = f{x,y)V{z,C). 
That's why 

AU = U (U' - , BU = fy (U^ - V^) , f. {Uz - V^) = b, U^-Vz = 0. 
Now the equation B^b — Bbx — takes the form 

f fxy fxfyi 

and the general solution of this equation is f{x,y) = g{x)h{y). Besides, the 
equation gh (V^ — Uz) = —b requires b[x, y, z, ^) = g{x)h{y)b"{z, so we get 

V^-Uz^ -b°. 

The relations obtained show how to build a solution of this class. Namely, 
first, we choose the function V{z, ^), then we determine the function U{z, ^) 

by 

C/(^,0= JVzdC + liz), 

where l{z) is an arbitrary function, which may be assumed equal to 0. After 
that we define b" — U^ — U^- The functions g{x) and h{y) are arbitrary, and 
for A and B we find 

A{x, y, z,i)=g {x)h{y) , B{x, y, z, = g{x)h {y)- 



U 



U 



In this way we obtain a family of solutions, which is parametrized by one 
function V of the two variables {z,^) and two functions g{x), h{y), each 
depending on one variable. 

For *{a A F + P A *F) wc obtain 

*{a AF + (3 A *F) = Audx - Budy — budz — bvd^ = 
= \{U'-V')[{ghf] dx-\{U'-V')[{ghf] dy- 
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-{ghf UWdz^dz - {ghf (I Vh^di^di - dx\ 
where the interaction energy-momentum tensor is defined by the matrix 



pu 



■\{ghfZ 

lighfz 

{ghfJUWdz 

{ghfJVWdi 



and the notation Z = U^ — V'^ \s used. For the full energy tensor 



rp3 



{ghf J Ub^dz-^iU^ + V^) 

Ts = -T! = {ghfUV, 
Tt^ighf \^jvh'^di + \{U' + V') 

and all other components are zero. 



3.3.2 Examples 

In this subsection we consider some of the well known and well studied (1+1)- 
dimensional soliton equations as generating procedures for choosing the func- 
tion V{z, and only the 1-soliton solutions will be exphcitly elaborated. Of 
course, there is no anything standing in our way to consider other (e.g. mul- 
tisoliton) solutions. We do not give the corresponding formulas just for the 
sake of simplicity. 

We turn to the soliton equations mainly because of two reasons. First, 
most of the solutions have a clear physical sense in a definite part of physics 
and, according to our opinion, they are sufficiently attractive for models of 
real physical objects with internal structure. Second, the soliton solutions de- 
scribe free and interacting objects with no dissipation of energy and momen- 
tum, which corresponds to our interpretation of the Frobenius integrability 
equations, as we explained in the preceding (sub)sections. 

1. Nonlinear equation Klein- Gordon. In this example we define our functions 
U and V through the derivatives of the function f{z,^) in the following way: 
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U = fz, V = f^. Then the equation — = fz,^ — f(z = ^ is satisfied 
automatically, and the equation Uz — = h° takes the form fzz — = b°. 
Since b° is unknown, we may assume 6° = b"{f), which reduces the whole 
problem to solving the general nonlinear Klein-Gordon equation when b° 
depends nonlinear ly on /. Since in this case V — we have 

I vb^im = / mm = / / ^wf] j b%f)df. 



For the full energy density we get 

1 



Tt^-{9hf[fl + fl + 2jbV)df 



Choosing b°{f ) = m'^sin{f) we get the well known and widely used in physics 
Sine-Gordon equation, and together with this, wc can use all solutions of this 
nonlinear equation. When we consider the (3-1-1) extension of the soliton so- 
lutions of this equation, the functions g{x) and h{y) have to bo localized too. 
The determination of the all 5 functions in our approach is straightforward, 
so we obtain a (3-1-1 )-dimcnsional version of the soliton solution chosen. As 
it is seen from the above given formula, the integral energy of the solution 
differs from the energy of the corresponding (l-l-l)-dimcnsional solution just 
by the (x, y)-locahzing factor {ghy. 

For the 1-soliton solution (kink) we have: 



f(z,^) = Aarctg l^exp 



±2m 



A = g'{x)h{y) 



^ch (z - 
±2m7 



±-{z--0 }, 7 
7 c 



'1- 



±2mw 



b° = U,-V^^ 



ch {z - ^C)] 
-2w?sh [d=^ (z - ^i) 



B = g{x)h\y) 



nch (z - H^) 

±2m7 



ch 



ch {z - He) 
and for the 2-form F wc get 



rp4 



{gKfA.m^ 



{z - fe); 



F = _^^!Mx)!^dy Adz+ ^ ^2m,(x)My) ^ 



^ch [±r^{z-^0 



c^ch ±^{z-^0 
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In its own frame of reference this soliton looks like 



ch{±mz) 

From this last expression and from symmetry considerations, i.e. at homo- 
geneous and isotropic medium, we come to the most natural choice of the 
functions g{x) and h{y): 

^ ch{mx) ' ch{my) 



2.Korteweg-de Vries equation. This nonlinear equation has the following 
general form: 

fi + aiffz + a2fzzz ^ 0, 

where Oi and 02 are 2 constants. The well known 1-soliton solution is 



V OoOl 



w c 

We choose ^(2;,^) = f{z,^) and get 

[4-^^ ^ \Lw cj 



2a, 



L cL 



Ti = {ghf 



ayil + L) 



2w^Lch^ 



- - — 



3. Nonlinear Schroedinger equation. In this case we have an equation for 
a complex-valued function, i.e. for two real valued functions. The equation 
reads 

ih + fzz + \f\'f^O, 
and its 1-soliton solution, having osciUatory character, is 



^exp [-i {2az + 4(^2 _ p^)^ _ 0)] 



ch {2pz + 8ap^ - 5) 
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The natural substitution f{z,^) — y/p.exp{i(f) brings this equation to the 
following two equations 



For the 1-soliton solution we get 



, <^ = - 



ch^ {2(3z + 9>a(3i - 5) 
We put U — p, V — —Ip^pz and obtain 



11 4/?^ 16«/?4 



c/i2 (2/?;^ + %a^i - S) ' "^^^ ch^ {2(3z + 8a/3C - 5) ' 

^ = g^(a:)%) , ^ ' B = g{x)h\y) 



16/53(16a2/?2 - l)sh{2(3z + 8a(3^ - S) 



ch%2Pz + 8aPC - 
n = {ghf 



.4 . ,^2 16/9' 



We note that the solution of our equations obtained has no the oscilla- 
tory character of the original Schroedinger 1-sohton solution, it is a (3+1)- 
localized running wave and moves as a whole with the velocity 4cq;. 

A.Boomerons. The system of differential equations, having soliton solu- 
tions, known as boomerons, is defined by the following functions: K : TZ^ —>■ 
71^, H : 71^ ^ TZ, and besides, two constant 3-dimensional vectors r and s, 
where / is a unit vector: |s| = 1. The equations have the form 

- s.K^ = 0, K^^ = H^^s + r X - 2K^ x (K x s). 

Now we have to define our functions U{z,^), V{z,^) and b°{z,^). The defin- 
ing relations are: 

U = H, = \s\'^H„ V = s.K, = (s.K)^, b° = -s. [r x - 2K^ x (K x s)] . 
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Under these definitions our equations — Vz — 0, V^c — Uz — —h° look as 
follows: 

[Hz - s.-Kz]z = 0, s. [K,^ - HzzS - r X K, + 2K, X (K X s)] = 0. 

It is clear that every solution of the "boomeron" system determines a solution 
of our system of equations according to the above given rules and with the 
multiplicative factors g{x) and h{y). The two our functions A{z,^) and 
B{z,^) are easily then computed. 

Following this procedure we can generate a solution to our system of 
equations by means of every solution to any (l+l)-soliton equation, as well 
to compute the corresponding conserved quantities. It seems senseless to 
give here these easily obtainable results. The richness of this comparatively 
simple family of solutions, as well as the availability of corresponding cor- 
rectly defined integral conserved quantities, are obvious and should not be 
neglected. In particular, it would be interesting to analyze the abilities of the 
&rea^/ier-solutions of some soliton equations as possible models of bounded 
systems of the type of hydrogen atom. It is worth to note that in this ap- 
proach the classical quantity potential is out of need. The proton and the 
electron participate equally in rights as relatively isolated subsystems of a 
more general dynamical system. The discrete character of the energy spectra 
trivially follows from the fact, that the transitions among the various station- 
ary states of the more general system are caused by creation or annihilation 
of photons, taking into the system or out of it the corresponding conserved 
quantities. As for the quantitative description of this spectra, probably, it 
will be achieved by a corresponding choice of the integrability constants. 
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